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Die Beugung elastischer Wellen an der Halbebene 
: Von A.-W. Maue in Freiburg i. Br. 


Das vorliegende Beugungsproblem wird zweidimensional behandelt. Der durch Überlagerung ebener 
Wellen gewonnene Lösungsansatz führt auf ein System linearer Integralgleichungen. Dieses reduziert 
sich nach Olemmow auf die Aufgabe, eine komplewe Funktion so in zwei Faktoren zu spalten, daß 
jeder Faktor einen vorgeschriebenen Anteil ihrer singulären Stellen enthält. Die Aufspaltung läßt sich 
mit Hilfe des Cauchyschen Satzes durchführen. 


The presented diffraction problem is treated two-dimensionally. The proposed solution that is gained 
by superposition of plane waves leads to a system of Linear integral equations. The problem can be reduced, 
as was shown by Olemmow, to the task of splitting a complex function in two factors in such a way, 
that each factor contains a prescribed portion of singular points. The splitting is carried out by aid of 
Cauchys integral formula. 

Le probl&me de diffraction en question est traitE pour deux dimensions. L’arrangement de solution gagne 
par superposition d’ondes plaines conduit & un systeme d’Equations integrales lineaires. Cela se reduit 
selon Ol em m ow au problöme de fendre une fonction complexe en deux facteurs, de maniere que chaque 
facteur contient une part prescrite des points singuliers de la fonction. La fonction peut &ire fendue & l’aide 
du theoreme de Cauchy. 


Ilpennaraemaan AnuhpakrunoHHaa IPoÖ1eMa paccMaTpaBaeTca B AByMepHoä TPakToBke. 
Teopema pemenns, NO1yYeHHAA IOCPeACTBOM CYNEPHOSUNHH ILIOCKHX BOJIH, IIDHBOAUT K CHCTeMe 
JAHHeÄHBIX HHTETPAJIBHEIX YP&BHeHHK. 3amaya cBonuTca, NO KuemmoBy (Clemmow), K TAKOMy 
PaCÖHeHHIO KOMINIEKCHOU DYHEIAH Ha BA MHOMHTEIA, YTO Ka’kıbIa MHOKHTENB CONEP:KAsı Öbl 
SaNaHHYIW MOM ee OCOÖBHIX MectT. PasÖueHne IPOBONKuTcH HPH NmoMomm Teopemzı Komm. 


Die zunehmende Bedeutung des Ultraschalls in Physik und Technik, sowie die Verwendung 
kinematographischer Methoden!) zur Beobachtung elastischer Wellen, lassen es wünschenswert. 
erscheinen, die Dynamik der elastischen Medien genauer zu studieren, als das bisher geschehen ist. 

Das einfachste, seit langem behandelte, dynamische Problem der Elasitzitätstheorie ist 
die Reflexion elastischer Wellen an einer freien Oberfläche?). Im Gegensatz zum Reflexions- 
vorgang der Optik treten dabei im allgemeinen z wei reflektierte Wellen, eine longitudinale 
(Reflexionswinkel &) und eine transversale (Reflexionswinkel ß) auf, wobei & und ß durch das 


Brechungsgesetz a: =7 verknüpft sind (a =longitudinale, b =transversale Schall- 


geschwindigkeit). Nächstes naheliegendes Ziel wird die Behandlung eines Beugungsproblems 
sein. Wir wählen als einfachstes Problem dieser Art die Beugung an der Halbebene: Ein un- 
endlich ausgedehntes elastisches Medium sei längs einer Halbebene aufgeschnitten; zu unter- 
suchen ist die Ausbreitung einer aus dem Unendlichen einfallenden ebenen Welle®). Dabei wollen 
wir uns auf die Betrachtung des zweidimensionalen Problems beschränken und annehmen, daß 
die Welle in einer zur Kante der Halbebene senkrechten Richtung einfällt. 

Die Beugung elektromagnetischer Wellen an der vollkommen leitenden Halbebene ist be- 
kanntlich von Sommerfeld*)in Strenge gelöst worden. Sommerfeld benutzt dabei 
zweideutige, längs der beugenden Kante verzweigte Lösungen der Wellengleichung und erfüllt 
die Randbedingung auf der Halbebene nach dem Spiegelungsprinzip. Es liegt zunächst nahe 
zu versuchen, diese Methode auf den elastischen Fall zu übertragen. Dabei stößt man jedoch 
auf Schwierigkeiten, die damit zusammenhängen, daß am beugenden Schirm nicht wie im 
optischen Falle (bei geeignet gewählter Polarisation der einfallenden Welle) nur eine, sondern 
zwei Randbedingungen (Verschwinden von Normal- und Schubspannung) erfüllt sein müssen. 


2) H.Schardin,D. Elle und W.Struth, Z.techn. Physik 21 (1940), S. 203; H. Schardin, 
Glastechn. Ber. 23, (1950), S. 1, 67 u. 325. 

2)Lord Kelvin, Phil. Mag. 26 (1890), S. 422; P. Drude, Ann. d. Phys. 41 (1890), S. 759; 
A.Schoch, Ergebn.d. ex. Natw. 23 (1950), S. 144; F.Sauter,Z. angew. Math. Mech. 30 (1950), S. 94. 

3) Anm. bei der Korrektur: Das entsprechende Problem für eine einfallende Stoßwelle mit Stufenprofil 
wurde von M. M. Fridman, Vorträge der Akademie der Wissenschaften der VdSSSR 60 (1948), S. 1145, 
behandelt, worauf mich Herr Schoch liebenswürdigerweise hinwies. Die dabei verwendete Randbedingung ist 
die der festgehaltenen Oberfläche. 

4)A.Sommerfeld, Math. Ann. 47 (1896), S.317;inFrank-v.Mises, Differential- u. Integral- 


eichungen der Mechanik und Physik, 2. Aufl., Braunschweig, 1934 Bd. II, Kap. 20; A. Sommerfeld, 
lenken über theor. Physik, Bd. IV, Optik, $ 38, Wiesbaden 1950. 
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Das Spiegelungsprinzip versagt daher im elastischen Falle. Im folgenden soll ein anderer Weg 
eingeschlagen werden. Der dabei gewonnene Ausdruck für die gebeugte Welle entspricht in- 
sofern dem bei Sommerfeld als die Welle wie dort formal durch Überlagerung ebener Wellen 
dargestellt wird. 

In den letzten Jahren sind eine ganze Reihe von Arbeiten erschienen, in denen Beugungs- 
probleme mit Hilfe von Integralgleichungen formuliert und gelöst werden. Es handelt sich dabei 
im wesentlichen um zwei mögliche Wege. Die eine Methode?) geht aus vom Huygensschen 
Prinzip und benutzt als unbekannte Funktion eine Belegung des Beugungsschirms. Im elektro- 
magnetischen Falle beispielsweise stellt die Belegung den elektrischen Flächenstrom im Schirm 
dar. Als Kern der Integralgleichung tritt die Green’sche Funktion des unendlichen Raumes 

ikr 
° _ auf. Die zweite Methode®) ist von Clemmo w und von Hönl entwickelt worden und 


stellt die gesuchte Lösung des Beugungsproblems durch Überlagerung ebener Wellen dar, 
die sich in allen möglichen Richtungen ausbreiten. Unbekannte Funktion ist die spektrale 
Amplitudenverteilung innerhalb des betrachteten Richtungsspektrums. Diese zweite Methode 
hat verschiedene Vorteile gegenüber der ersten, weshalb sie im folgenden angewandt werden 
soll. Einmal ist die Lösung der Integralgleichung leichter durchführbar; zum anderen liefert 
die Amplitudenfunktion ohne weiteres die angulare Intensitätsverteilung in weitem Abstand 
vom Beugungsobjekt. 

Da es uns bei den vorliegenden Untersuchungen in erster Linie darauf ankommt, die all- 
gemeinen Verhältnisse und grundsätzlichen Schwierigkeiten kennenzulernen, mit denen in der 
Beugungstheorie elastischer Wellen zu rechnen ist, begnügen wir uns damit, die Rechnung so 
weit zu führen, bis das Problem im Prinzip gelöst ist. Auf eine numerische Auswertung haben 
wir zunächst verzichtet, da die hierzu notwendige Rechenarbeit erst dann lohnend erschiene, 
wenn die Möglichkeit eines Vergleiches zwischen Beobachtung und Theorie gegeben wäre. 


$ 1. Die Randwertaufgabe 


1. Differentialgleichungen. Wir betrachten Schwingungen des elastischen 
Kontinuums mit der Frequenz w. Es sei «a die longitudinale und 5 die transversale Schallge- 
schwindigkeit. Dann sind die Wellenzahlen der beiden Wellentypen 


(0) © 
as Keen he el 
Die zugehörigen Wellengleichungen für die Verschiebungen 3, und 3, der longitudinalen und 
transversalen Wellen lauten 
43, +53 =0, rot 3, =0 (2) 
43, +K?3,—=0, div, —=0 ® 
oder nach Einführung eines skalaren Potentials ® und eines Vektorpotentials und Zusammen- 
fassung der beiden Bestandteile 


A8+8=0, AE+KRY=0| 
Dabei ist über die bei der Einführung der Potentiale bestehenden Freiheiten schon teilweise 


verfügt. Aus (3) entnimmt man, sofern man noch div #=0 festsetzt, daß ® und %# bis auf 
konstante Faktoren die kubische Dehnung und die lokale Drehung angeben: 


dv3=k!®, Tot 3 = KU NEE RB ELCH 
also einen einfachen physikalischen Sinn haben. 
2. Randbedingungen. Die Verschiebung erfolge in der x-y-Ebene. Dann 


hat # die z-Richtung und kann als skalar betrachtet werden. (3) lautet in Komponenten .ge- 
schrieben 


od ow ad 0w 5) 
dm Do = EEE 


5)G.Jaffe, Phys. Z. 22(1921), S. 578; H. A. Bethe, Phys. Rev. 66(1944), 8.163; E.T.Copson, 
Proc. Roy. Soc. A 186 (1946), S. 100; Quart. J. Math. 17 (1946), S. 19; H. Levine und J. Schwinger, 
Phys. Rev. 74 (1948), S. 958; 75 (1949), 8.1423; A. Sommerfeld, Vorlesungen (l.c.) Bd. IV, $ 39; 
A.-W. Maue, Z. Physik 126 (1949), S. 601; Cl. Müller, Abh. dtsch. Akad. Berlin 45/46, Nr. 3; Math. 
Ann. 123 (1951), S. 345. Weitere Literatur bei Clem mo w®), 

°)P.C.Clemmow, Proc. Roy. Soc. A. 205 (1951), 8. 286; H. Hönl, Z. Physik 131 (1952), S. 290; 
E.Groschwitzu.H.Hönl,Z. Physik 131 (1952), S. 306. 
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Bild 1 zeigt die längs der negativen x-Achse aufgeschnittene x-y-Ebene. Als Randbedingung 
haben wir beiderseits des Schnittes das Verschwinden von Normal- und 
Schubspannung zu fordern: 


en eletun =, EUR TEREIG). 


Die Spannungen ergeben sich aus der Spannungsdehnungs- 
beziehung, die bei der Dichte o und Verwendung von a und 5b 
als elastischen Materialkonstanten liefert: 


dE 7 (& dE ) 
en aurER, ame ER >. u wel 22 
DEE (@ eb) a oe m ed det 0Y (M). Bild 1 
Führen wir hier mittels (5) die Potentiale ein und formen mit Rücksicht auf (1) und (3) um, 
so geht (6) über in 


m (re 

2ob? = dw = mt 9 Sm oy == | N 

Gy BEI) 92 K: UI EISEN (8). 
a Intl t 5)? =0| 


Zu diesen Randbedingungen tritt die sogenannte Kantenbedingung”), die das Verhalten 
der Welle in der Umgebung der Kante betrifft. Sie sorgt dafür, daß an der Kante keine Kräfte 
auf das elastische Kontinuum ausgeübt werden, und läßt für die kubische Dehnung und lokale 
Drehung und wegen (4) auch für die Potentiale ein Unendlichwerden an der Kante (r = Ab- 
stand von der Kante) wie 

DEREN LTer (a nee (9) 
zu®). 

Ferner ist de Ausstrahlungsbedingung im Unendlichen zu erfüllen. 


$ 2. Die Integralgleichungen 
1. Lösungsansatz. Wir setzen die Potentiale der gesuchten Welle aus mehreren 
Bestandteilen zusammen: 
®=9+9+9,,. EEE DER EEE er (10). 


Hierbei beschreibe ®,, 7, die durch den Beugungsschirm unbeeinflußte einfallende Welle. Die 
restlichen Glieder beschreiben den Einfluß des Schirms, also Reflexion, Schattenwirkung und 
Beugungserscheinungen. Die Aufspaltung in die beiden Bestandteile mit den Indizes 1 und 2 
erfolge dabei so, daß der erste bezüglich einer Spiegelung an der &-Achse symmetrisch, der 
zweite antisymmetrisch ist. Für die Potentiale gilt dann, wenn man beachtet, daß die zu 7 pro- 
portionale lokale Drehung im symmetrischen Falle antisymmetrisch ist und umgekehrt, 


DR Y-Oa@y), ABR-N=- FR | Ron 
EN) Y), By) Aay) | 


Erfüllen wir jetzt die Randbedingungen (8) auf der oberen Seite des Schirms, indem wir dort 
verlangen 


(ra (rem © 
2-45) 0 | 
a al) 20 (e) 
elta) nu tan 22 0 @ 


so sind diese Bedingungen auf der unteren Seite des Schirms wegen der Symmetrieverhältnisse 
(11) von selbst erfüllt. Wir dürfen uns daher darauf beschränken, die Funktionen ®,, Y;(W=1, 2) 
in der oberen Halbebene y> 0 zu untersuchen. Es ist dann noch erforderlich, die Stetigkeit 


?) Für den Fall skalarer und elektromagnetischer Wellen vgl. J. Meixner, Z. Naturforschg. 3a (1948), 
S. 506; Ann. d. Phys. (6) 6 (1949), S. 2; A.-W.Maue(l.c.) $3. 
| .®) Vgl. hierzu A.-W. Maue, Z. Naturforschg. 7a (1952), S. 387, Gl. (17). 
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der Funktionen und ihrer ersten Ableitungen längs der positiven z-Achse zu gewährleisten 
und zusätzlich zu (11) zu fordern 


od, w_g 
a 
r DE AU De a NEE (13). 
0%, 
d,=-0, Zel 
°Y 


In der oberen Halbebene setzen wir die Potentiale als Überlagerung ebener Wellen an, die 
die Wellengleichungen (3) erfüllen: 


P12 = N (MeirerVEsrN dA 


R fan VOR UN 
?12= [9.2 aetTE RN) ER 


Das positive Vorzeichen vor der Wurzel sichert da- 
bei für Al<k, K den Ausstrahlungscharakter der 
Wellen (Zeitabhängigkeit e=i®t). Für Al>%,K ist 
die Wurzel positiv imaginär zu wählen, entsprechend 
Au EBADE einem exponentiellen Abfall der Wellen nach pogi- 
Bild 2 tiven y hin (Oberflächenwellen). Diese Vorzeichenwah 

erfordert, daß der Integrationsweg in der in Bilde 

dargestellten Weise um die Verzweigungspunkte A = +k, + K herumgeführt wird. 


Bildung der Ausdrücke (12) für y = 0 ergibt 


© 
(1,2) | 


y = | IE = Pu2()+AYR? —RQ>0)| eirz dA | 
Be . (15). 
a _ A ya Pı2(4)+ e — ) OR) etz di | 


Hiervon müssen I y, 2) und Ba für beliebige Werte von & verschwinden und zwar für <<0 
wegen (12b, c) und für & >0 wegen (13). Daraus folgt das Verschwinden der entsprechenden 
beiden Fourierkoeffizienten, d. h. der geschweiften Klammern, in (15). Dem tragen wir ER SrhUnUR, 
indem wir unter Einführung zweier neuer Funktionen R,, (4) ansetzen 


Pra-(—) Ro,  MQ=-AyFZERA), 


But) = -AYRı— ZERO), 0) (5 —r) Rd) 


2. Integralgleichungen. Aus (12a, 2 (15) und (16) folgt die erste Integral- 
gleichung für R,, (A) 


oo 


x: Me co 
ie: Ha VRR) RR) Rısd) ers di = — 1° für 2<0 (7) 
— oo 823; y=0 
Als zweite Integralgleichung fordern wir 
[ve je R,(A) etz di} = 0 
für. > Ü.Swee Rue (18) 


VRR: R,() rad —0 


und erfüllen damit die durch Einsetzen von (14) und (16) in (13) entstehenden Bedingungen. 
Um das einzusehen, hat man neben (18) deren erste und zweite Abteilung nach x zu verwenden. 


et D 
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Die Lösung des Beugungsproblems läßt sich nach (10), (14) und (16) mit Rücksicht auf die 
Symmetrieverhältnisse (11) für die ganze x-y-Ebene in der Form schreiben 


es ) 
Barker I — 2) RO) FAyR? a By) et VE=RD 2, | 
den (19) 


v_w (| RER: K° 2 | (ax +VR?—3®|y| 
-%+ EA —®eR@A)+\5 —#) BR eit“t A 


wobei sich das doppelte Vorzeichen auf die beiden Halbebenen y>0 und y< 0 bezieht. 


Wir spezialisieren nun die einfallende Welle. Der Einfall erfolge schräg von unten. Wir 
unterscheiden dabei zwei Fälle: 


I. Einfallende Druckwelle. 


D, = eeHe Zen, pP —=0 eh ER (20 a). 
II. Einfallende Schubwelle. 
D=0, P = eiwatVR' ey) —Ka@n<Ks AR 22221(20b); 
Hiermit ergibt sich für die rechte Seite von (17) nach (12a, d) 
5 — I) | = Aeir:, a Ares.) Diary, (21), 
y=0 y=0 
wobei 
E: 
I A, = Se). 4, — ayk? —n: 2 Se (22a) 
a K: ; 
II, A=—ayYE?’— x, A=—(—r) EA RR (22b). 


Die Integralgleichung (17) lautet dann 


oo 


| IE: a) Hay RCRT A| Rs Merrdr = Ayere füra<0 (2). 


$ 3. Lösung der Integralgleichungen 


1. Lösungsverfahren. Die beiden Integralgleichungen für R,,(4) sind (18) und 
(23). Wir lösen sie nach dem Vorbilde von Clemm o w®) mit Hilfe des Cauch yschen Satzes 
und betrachten zunächst R,. Es seien U, (A) und U,;, (A) zwei Funktionen mit algebraischem 
Verhalten im Unendlichen der komplexen /-Ebene (vgl. hierzu Ziffer 2 dieses $), von denen die 
erste in der positiv-imaginären, die zweite in der negativ-imaginären A-Ebene regulär sei. Die 


5 Begrenzung der beiden komplexen Halbebenen verlaufe dabei längs des Integrationsweges 
von Bild 2, wobei neu festgesetzt werde, daß der Punkt A = x oben umlaufen wird. Setzen wir 

nunmehr | 
EEE 2 A) ES A EEE (24a) 


K® ) ; FRE) Bel ch Une) 
( 2 ES) ur 2ni A—x U) 
und verlangen U,,(x) #0, so löst der Ansatz (24) beide Integralgleichungen. Der Integrations- 
weg in (18) läßt sich nämlich wegen der vorausgesetzten Eigenschaften von U, (A) ins Unendliche 
der positiv-imaginären A-Ebene verschieben, wo das Integral wegen des Verhaltens der Expo- 
nentialfunktion ei#* und 2>0 verschwindet. Entsprechend kann der Weg in (23) ins Un- 
endliche der negativ-imaginären A-Ebene verlegt werden, wobei der Pol A=x von (24b) die 
rechte Seite von (23) liefert. 

Es ergibt sich nun die Aufgabe, U, und D,;, so zu bestimmen, daß (24a) und (24b) mit- 
einander verträglich sind. Elimination von R,(A) führt zu 


1 a ee N) 
ar a) ne a A) Hare ze 

K:—k: 

| | Er 


A, . . (4b) 


(25) 
(A —a)YRER8 : F (9) 


angew. Math. Mech. 


6 Maue, Die Beugung elastischer Wellen an der Halbebene „4.33 A 1/2 Jan./Febr. 1953 
mit 
En 
2 2 
ee Li hr . (26). 
K:—k: Y(k® — 22) (K? — 2°) 


Die Abspaltung eines Faktors F (A)erfolgte dabei so, daß in der gemäß Bild 3 von —K bis —k 
sowie von k bis K aufgeschnittenen A-Ebene gilt 


conste 5, 
FiA)>1+ 7: für A200. - en al 2005 
wie jetzt zu zeigen ist. Hierzu stellen wir zunächst fest, daß 
Authene F(A) in der aufgeschnittenen Ebene eindeutig ist. Weiter gilt, 
wie wir erkennen, wenn wir A beispielsweise längs der reellen 

Bild 3 }-Achse anwachsen lassen, 


VER’) = VRR HEHE N für Aw (8) 


woraus durch Entwicklung für große A direkt (27) folgt. 

Um U,(A) und U,,(A) zu finden, muß (25) in zwei Faktoren zerlegt werden, von denen 
der eine in der oberen, der andere in der unteren A-Ebene regulär ist. Für F(A) sei eine solche 
Zerlegung 


PA PlA) ER IA) e RR RER (29) 
vorgenommen, wobei entsprechend (27) gelte 
FF, >11: für AS near (30). 


Dann folgt mit Rücksicht auf den Verlauf der Trennlinie zwischen oberer und unterer /-Ebene 
gemäß Bild 2 und wegen U;; (x) #0 unter Absehung von konstanten Faktoren 


1 TER 
—— m(A—r): . ET IE rs aa 
Ur REF (N ee ee ee Wa (31b). 

Aus (24b) und (31b) können wir jetzt R,(A) direkt entnehmen. Ferner läßt sich R,(A) angeben, 


indem man A, durch 4, ersetzt und in den Wurzeln k und K miteinander vertauscht. Man 
erhält 


1 == (), = 
RR.) = Fu : a —— . (32). 
 2ri 3 SR F 1) > 2) +42 Y(k?— 22) (K?— 22) 
—H 


Setzen wir diese Ausdrücke in (19) ein und unterscheiden wieder die beiden Polarisations- 
fälle der einfallenden Welle, so ergibt sich für die vollständige Lösung des Beugungsproblems 
mit (20) und (22): 

I. Beugung der Druckwelle. 


D— ei ası+ Ve a EEE he: FA), ei(au+Vk®—32|y|) 


2n%V K—x Fu() 
5 ER 33 
EM Seren | 
BE Fı(A) eiAa+/YR?— | 
Ze Ve N 
(33b) 


FERNE 


(5-2) + VRR) 


. 
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II. Beugung der Schubwelle. 


oo 


UNTEREN, ) 
dd — KR | K—A ‚Fır(A) . eiaa+ Vr®—a: |y|) 
and, k—x Fı() 
RE. KR: | 33 
& 2) REF + — re) Ay KH 5 (330) 
Y—_ us Een. I Er . APR |y)) | 
el Pez K: 33d). 
tray SEHE He 5, | 
x ® 
A—ı 


2 — ») +22 Y(k® — 22) (K® — 28) 


2. Konvergenz und Kantenbedingung. Zur Konvergenz der Integrale ge- 
nügt es, daß der Faktor der Exponentialfunktion im Integranden für A— © gegen O geht, da 
dann die Exponentialfunktion für Konvergenz sorgt. Die Bedingung ist erfüllt, wie sogleich 
gezeigt wird. 

Eine Ausnahmestellung nimmt in dieser Hinsicht der Punkt = y =0 ein, und es bleibt 
zu untersuchen, wie sich die Welle in der Umgebung dieser Stelle verhält. Zu diesem Zweck 
bemerken wir, daß wegen (28) gilt 


(5) HET SR für arm 2... 


und weiter wegen (30) ein Verhalten des Faktors der Exponentialfunktion in den Integranden 
von (33) wie A-12. Der Integrationsweg (vgl. Bild 2) wurde so um die Verzweigungspunkte 
herumgeführt, daß gilt 

Ve—2>ilil, YE®-2oi1| fü AoLt® ...... (35). 
Für den gesamten Integranden ergibt sich somit das Verhalten A-!/? e#=—llvl. Hieraus folgt 
für die Integrale ein Verhalten wie r—V/? in der Umgebung der Kante r = 0 in Übereinstimmung 
mit der Forderung (9) für die Potentiale. Die Kantenbedingung ist also erfüllt). 


$ 4. Das Aufspaltungsverfahren 


Die Produktaufspaltung (29) führen wir durch, indem wir den Logarithmus von F(A) in 
zwei Summanden zerlegen. Dementsprechend setzen wir 


Ben a Wera. 3...) 
Wir schreiben zunächst f(A) mit Hilfe des Cauch y’schen Satzes als Integral 
1 dz 
= = — EN IE TE (DL), 
J@) =1In F(A) 5m $inFß) RR) (37) 


Der geschlossene Weg umschließt den Punkt z=4. Der Integrand ist wegen (27) im Unend- 
lichen regulär und kann daher auf die beiden Wege C; und Cr: (Bild 4) zusammengezogen 
werden, die gemeinsam sämtliche singulären Stellen von 
In F(z) umschließen sollen. Die gesuchte Zerlegung er- .a,-x 
gibt sich nunmehr von selbst: 


1 d = LDEen, 
5 Au) = 57; [nr@- = 142(88), z-Ebane 


Bild 4 


SAL 
(A) und F, (A) sind überall außerhalb des Weges C,, insbesondere in den Punkten A=k,K, 
regulär und eindeutig. Entsprechendes gilt für fr(4) und F,(A) bezüglich der Punkte 


9) Tatsächlich sorgen die Integralgleichungen (17) und (18) von selbst für Erfüllung der Kantenbedingung, 
vgl. hierzu H. Hönlu. A.-W.Maue, Z. Physik 132 (1952), S. 569. 


. . 5 . th. h. 
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A=—k,—K. Die Forderung der Regularität oberhalb bzw. unterhalb des Weges von Bild 2 
ist daher erfüllt, wenn wir den Weg noch so deformieren, daß er links oberhalb C, und 
rechts unterhalb C\, verläuft (Bild 5). 

Um die Wege C; und (\, im einzelnen festzulegen, suchen wir die singulären Stellen von 
In F(z) auf. Das sind einmal die Verzweigungen von F (z), d.i. die Punkte z= tk, +K, 
zum andern die Null- und Unendlichkeitsstellen von 
F(z). Aus (26) erkennt man, daß die Unendlichkeits- 
stellen mit den Verzweigungen zusammenfallen und 
n den Nullstellen gelten muß 


Bild 5 


2 2 
= — + 22 /(k® —z2) (K?— 22) =0 (39). 
Der linksstehende Ausdruck ist dabei in der gemäß Bild 3 aufgeschnittenen z-Ebene eindeutig. 
Ordnen wir, entsprechend GI]. (1), der Wellenzahl z mittels 


= — 40 
=! (40) 
eine Phasengeschwindigkeit c zu, so lautet die Bedingung (39), in den drei Geschwindigkeiten a, 
b und c geschrieben, sofern man c reell und kleiner als b annimmt: 


(5) pe (ao (a1). 


Das ist aber gerade die Gleichung für die Geschwindigkeit c =czr der Rayleigh’schen Ober- 
flächenwellen !°). Gehen wir wieder zur Wellenzahl z über, so haben wir damit als gesuchte Null- 
stellen von F(z) 


PER. 2 RS RE ee (42). 
CR 
Die Integrationswege CO, und C,; umschließen diese Punkte (Bild 4). Der Beweis dafür, daß 
(39) neben (42) keine weiteren Lösungen besitzt, sei nur angedeutet. Durch Multiplikation von 
(39) mit 


KE: an > 
E —) Say oe . (43) 
gewinnt man eine Gleichung dritten Grades für 2°. Durch Einsetzen eines speziellen Zahlen- 
wertes für das Verhältnis X:k läßt sich zeigen, daß von den drei Lösungen dieser Gleichung 
nur eine (39), die beiden anderen aber (43) befriedigen. Hieran kann sich durch kontinuierliche 
Variation des Zahlenverhältnisses X:% nichts ändern. 

Zur Vorbereitung für eine numerische Auswertung von (38) ist es zweckmäßig, das Integral 
umzuformen. Der Integralwert ändert sich nicht, wenn im Integranden ein innerhalb der Schleife 
O,, 77 reguläres additives Glied weggelassen wird. Wir können daher schreiben 


TE A ee. de ER lady 
Lu 2mi 2? Y(k® — 2?) (K?— 2) gan | 


0,11 


Aus demselben Grunde ist auch die Vieldeutigkeit des Logarithmus im Integranden belanglos. 
Das Schleifenintegral läßt sich noch um die Punkte Fk, =K, FA, zusammenziehen 
und als reelles Integral schreiben. Die Rechnung ergibt AT 


K: N 2 
in a Em dz | 
A ma ar engen ae (AR 


rk 


$ 5. Die Beugungserscheinungen 


Es soll nunmehr das asymptotische Verhalten der gebeugten Wellen (33) in weitem Ab- 
stand von der beugenden Kante untersucht werden. Hierbei müssen sich neben der einfallenden 
Welle fünf Bestandteile ergeben, drei ebene Wellen und zwei Zylinderwellen. Die ebenen sind 


10) Vgl.etwa A. Sommerfeld, Vorlesungen (l. c.), Bd. II, 1947, $ 45.26. 
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die beiden reflektierten Wellen auf der Vorderseite des Beugungsschirms und eine dritte Welle 
hinter dem Schirm, die hier die einfallende Welle kompensiert und so den Schatten erzeugt. Die 
beiden Zylinderwellen, eine longitudinale und eine transversale, liefern die’ eigentlichen 
Beugungserscheinungen. 

In Bild 6 sind die Strahlenrichtungen der einfallenden und reflektierten Wellen eingezeichnet 
und durch die Winkel a, ß festgelegt. Von den einfallenden Wellen existiert jeweils nur eine, 
während die reflektierten Wellen beide gleichzeitig vorhan- 
den sind. Für die Winkel gilt 


RER DIR Sin DB .: 21.0.8507 (A6), 


Die Beugungswellen in (33) haben die allgemeine Form 


D = | Pi) east lr ln aa, | 


£: a (SEN 
v- [onen | 
Durch die Substitutionen BildER 
A oe Sin On ne ee (48) 
und ; 
t=rsinp, yl =rcosp 2. <p<? er (Ay) 
(Bild 6) entsteht hieraus 
a/2 —i oo 
D, = f  P(ksiny) et" cs 6 To cosydy, 
eh | ir . (80). 
MR j‘ Q(K sin ö)eikr os —pcosödd | 
— z/2 +io J 


Zur Durchführung des Grenzüberganges r— co haben wir den Sattelpunkt des Integranden 
aufzusuchen, der ersichtlich bei 


BU DEN a NER er 451) 
liegt. Das Integral über den Sattelpunkt liefert bis auf unwesentliche Zahlenfaktoren 
®,—k:P(ksin op) cos ae ı pP, SR :Q(Ksin p) cos en ID2); 
Ykr YKr 


d.i. die auslaufenden Zylinderwellen. Um die beiden Intensitäten miteinander vergleichen zu 
können, haben wir an Stelle der Potentiale die Verschiebungen 3 zu betrachten, und zwar ist 
die (zeitlich gemittelte) Energiedichte proportional zu |$|”. Beim Übergang zur Intensität 
(Energiestrom) ist noch mit der für beide Wellen verschiedenen Schallgeschwindigkeit zu 
multiplizieren. So erhält man unter Weglassung von beiden Wellen gemeinsamen Faktoren 
für die relativen Intensitäten mit Rücksicht auf (3) und (1) 


k? 2 } 2 

ko ER I Pksing), Kan rg Keine) (8). 

Das ist das gesuchte Beugungsbild. . | 
Weiter bleibt zu beachten, daß der Integrations- 


weg in Bild 2 oberhalb des Poles A =x herumläuft. Es x A, 
sind die beiden Fälle zu unterscheiden, daß der Sattel- 
punkt in der A-Ebene A, = ksin g bzw. Ksin o links oder 
rechts von x liegt: N 
n p<abzw.ß, h<x h PA (54) } Bild 7 
2. o>abzw.ß, h>r i 


Da, wie man sich leicht überlegt, der Weg in der A-Ebene durch den Sattelpunkt unter 45° 
von links oben nach rechts unten verläuft, bleibt er imersten Falle (Bild 7a) am Pol A= » hängen, 
im zweiten Falle nicht (Bild 7b). Im ersten Falle liefert also das Residuum in A = x einen Bei- 
trag zum Integral. Aus Bild 6 erkennt man, daß g <abzw. ß gerade der Bereich ist, in dem 
Reflexion und Schatten auftreten. 


Te 
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Wir betrachten zunächst den Fall y>0, d.h. oberes Vorzeichen in den Formeln (33). 
Der in negativem Sinne umlaufende Pol A = xliefert in (33a, d) einen Beitrag, der die einfallende 
Welle genau kompensiert, während er in (33b, e) nichts beiträgt. Für y<0 unteres Vorzeichen, 
liefern alle Integrale einen Beitrag. Die Formeln lassen sich in einfacher Form schreiben, wenn 
man die Winkel «, ß, 9 gemäß (46) und (49) einführt. Es wird dann 


2 2 ee Tg 
EZ os, Re = sin2ß, 
LIU K: . (55) 
x yYk? — x: =7sin2ß .tgß-ctga 
) 
und aus (33) entsteht 
I. Reflexion der Druckwelle. 
ikrcos(p+% tg? 2ßtgß —tga ikr cos (p —& 
EUR (pP +0) k (p ) r 
o=e METTPTITTERTTE . (56a), 
I 2tg2PBtgß lan (P—B) ö (56b). 
tg’ 2ßtgß-+iga 
II. Reflexion der Schubwelle. 
De 2 tg 2 ß ig & eikr cos (P—&) Y (56c), 


tg?2ßtgß +iga 


y _ gKree(eTP) tg 2PIgßP—tga irrcos wm EEE RE E (56.d) 
tg’2PßtgPp +tga 


in Übereinstimmung mit den bekannten Verhältnissen bei der Reflexion elastischer Wellen 
an einer freien ebenen Oberfläche). 


Liegt ® innerhalb des Winkelbereiches der Totalreflexion für Schubwellen, d.h. 
Ksin |p|>k, so bleibt der Integrationsweg in (33b, d) auch an einem der beiden Verzweigungs- 
punkte + k hängen. Der Beitrag des entstehenden Schleifenintegrals um diesen Punkt ver- 
schwindet jedoch wie 1/r für r— co, ist daher für das asymptotische Verhalten der Wellen 
belanglos. 


Eingegangen am 14. Februar 1952. 


Theorie der Druckstabilität der Sandwichplatte. II. 


Von H. Neuber in Dresden 


Zweiter Teil der in dieser Zeitschrift, Bd. 32, 8. 325— 337 erschienenen Arbeit. Auswertung der Theorie 
für die technische Anwendung, Diagramme für das optimale Wandstärkenwerhällnis. 


Second part of the paper, publicated in this jorurnal, Bd. 32. 8.325—337. Development of the results 
for technical applications, diagramms for the optimum ratio ofthickness. 


Deuzxieme partie de la theorie (la premiere partie &tait apparue dans Bd. 32, $. 325—337. du journal). 
Explanation des resultats pour l’application technique avec des diagrammes pour les dimensiones optimales. 


Bropas yacTb B 9TOM »kypHase, Tom 32, cTp. 325—337 onyÖsukoBaHHON pa6oTsI. Mcnon3o0- 
BAHHe TeOpHÄ AA TeYHHYeCKOH IIPHMEHHMOCTH C AHATPAMMAMH AA Iy4LuX usMepenuf. 


6. Linearisierung 


Da die Zusammendrückung der Platte praktisch nur in jenem Bereich untersucht zu werden 
braucht, der unterhalb der Bruchfestigkeit liegt, kommen für die Auswertung nur Werte von & 
und damit von Sin Betracht, die sehr klein gegenüber 1 sind. Mithin ist eine Entwicklung sämt- 
licher Ausdrücke nach Potenzen von $ zweckmäßig. Wie nachstehend gezeigt wird, ergibt sich auf 
diese Weise für die Stabilitätsbedingung eine Potenzreihe, die bei der dreischichtigen Platte mit 
der dritten, bei der zweischichtigen Platte mit der vierten Potenz beginnt. Wird jeweils soweit 
entwickelt, daß noch der Faktor der nächst höheren Potenz exakt herauskommt, und werden alle 
weiteren Potenzen vernachlässigt, so folgt für $ jeweils eine lineare Gleichung, welche für die 
Auswertung wesentliche Vorteile aufweist. Zunächst seien die Größen _?, y und u, sowie die in der 


weiteren Rechnung auftretenden Kombinationen dieser Größen nach Potenzen von $ entwickelt. 


Mit Hilfe der Binomialreihe erhält man; 


© 
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1+S 
2 FR 14254284284 --, 


e-I+(EYE-IHS+ SE + LS +, 


y=l-SUHL2SETI=IH+242)EH+U-2IR 2) +, 

B=Ü4-91—(A-2)8]1 +1 28 
=1+(—-2+20)8+ (1— 20) (2 — 2c) 8?+ 

B-1)Wt+R—D-4U- NEUN, (102) 


?-1-84(-5+2) 845 -2)8°+4. 
(2=1-254- 2441984 


ne=1+(-1+2)5+{-g+2)84., 


my+u—2=(—-4+40)5+(—-2+208°+ 
» Für die Hyperbelfunktionen ergibt sich mittels der Taylor-Reihe: 
Sin go = Sin [d + (e — 1) PB] = Sinf+(e—1)BCof + -- 


— Einß+ SPC + 82 (BCof + A Sin) + 
22 


1 £; . (103). 
+ 288008 + 3 Sinß+ BCfB)+:--, 


1 
Col oA= Co ß + SP Sin +z(B Sinß+BCojP)-+:-: 
Bei Anwendung dieser Entwicklung steht der Auflösung der Basisfunktionen nichts im Wege. 
Bei ®, zeigt sich, daß die Reihe mit S? beginnt, so daß an sich bei ©in e £ und Co] 0 £ eine Ent- 


wicklung bis zur dritten Potenz von S erforderlich wäre. Dies läßt sich jedoch vermeiden, wenn 
eine Aufspaltung von ®, vorgenommen wird, wie z. B. die folgende: 


9,(d)= (up ta — 2) (CE Pf — E23) + (u —1)u+y—D+LP,(ß)D,(P) (104). 


Für die kürzere Darstellung des Ergebnisses seien folgende Bezeichnungen eingeführt: 


M-5ITNd+SHN, MM -ZBM+SDD, ) 
A-HSTE UD-TndH SB | - (109) 


1 
D; () = 4: 5 [Ddo(P)+ 5 PD; (#)] 
Bei Unterdrückung der Zwischenrechnungen ergibt sich: 


Du(h)= Ein2B—2ß, Bu(d)= (3-20) Sin20+B6|2%. 


Du (B) = Ein2ß+2%, Du (B)= (5 — 2) Ein 24 + C0f20. 

Pu (P)= (3 —4c) Sin2P+P, 

Du (9 =(3 —20)(Sin2f + 2BCof2P)+ 2%, | ... (106), 
Du(B)= 3 —4o) Sin2B—2P, 

ud) = (5 —20)(in 26 + 24602928, 


Du (B) = (6 — 80) Co AB + 10 — 2c-+ 16024162, 
Bs1 (B) = (3 — 4) [(4d — 4c) Cof 4 ß+4L Sin4ß]+4—4c-+ 16. £? 


angew. Math. Mech, 
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Nunmehr läßt sich die Stabilitätsbedingung in eine für S lineare Beziehung überführen. Die Auf- 
lösung nach S führt auf folgende Formeln: 


Erste Linearisierung (S <]) 


A. Dreischichtige Platte mit gegensymmetrischer Verformung. 


aM On Bu 3 )+ SC HD HH)+ 
BR Je 


Linda) Du )o.(*) Da (R) 


+4H(H 1) |0.(4)\0.(2)02(0+ 2u(2) 019) + 21 (2) old) 0.0] 
(107) 


B. Dreischichtige Platte mit symmetrischer Verformung. 


MB) Bol) + HR + 


RER: Pa Balz) + Po mBulz) + 180 -or Hp ++ DaB HAN 


=E AH HN.) Di 9 Pa (%) 


Rp ERENPN SERIEN, 
(108), 


C. Zweischichtige Platte. 
a —H)D, (8) Du (f) Da) + 


en BE 9. (B) Du (B) Ds, Sr 2 (R) Ba: (P) 245 an De] 


+16 0? 0. +) ot) + 
+10 [o,(e+2)0u(6+5)+ 0u(e+ "ou (e+2)+ 0u(# +2) on(a+2)]+ 


LEÄH— Due) ®, (3) Ds(P) 


+ Hu —1)0,(2) Io, on + Du )®, (8) + ls ) Due Du (8) 
(109). 


Da die Anforderungen des Leichtbaues — wie erwähnt — ausschließlich auf Kombinationen von 
dünnen, festen und damit tragenden Schichten mit dicken, weniger festen und besonders leichten 
‚ Füllschichten führen, kommen für die praktische Anwendung nur Werte von A, in Betracht, die 
sowohl gegenüber der Wandstärke A, als auch der Knickhalbwelle 1 als klein angesehen werden 
können. Der Parameter £, bleibt daher auf Werte beschränkt, die sowohl gegenüber £, als auch 
gegenüber 1 klein sind. Hieraus ergibt sich ein zweiter Linearisierungsprozeß, indem die Stabili- 
tätsbedingung außer nach S auch nach ß, entwickelt wird. 


Mit Bezug auf die Gl. (107), (108) und (109) sind hierzu die Potenzreihen von ®,, ©, und do, 


erforderlich. Bei Anwendung der Potenzreihen der Hyperbelfunktionen ne aus den n Gleichungen 


(106) für das Argument £,/2: 


a- Mo. (9 0.(E) +o.004(@)]+ 100-0 0+ m) +onC+ + 
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ul) 4R+ u P-+:-- 0) (2)-@-20)8 ı18-_20@ a 

p) 6’ı Tpporı » u\9 RSIE TE c)Ar > 

bı 1 ß 1 
at Rt le) @-29h4+4B- 294, 


(110). 


ar twa-nnt, 


Du (2) - 151 —c? +81 —-)5—-4)f-+'-: 
Für die Gleichungen (107) und (108) werden außerdem Funktionen mit dem Argument (B + %ı) 
benötigt. Die zugehörigen Entwicklungen nach £,; sind: 


Do +B)= Duo (d)+2hı (Ef 2% —1)+., 

Du +A)= Du (f)+ 2,12 —2c)00f 2 +PSin2ß]+. a ae en eh) 
Da (+ B)= Dad) +2 (E28 +N +, | 
Da P+h)= Pad) + 2A [Ra — 20) Co 2f+PESin2ß] + ) 


Schließlich treten in Gl. (109) auch die entsprechenden Funktionen mit dem Argument (# + £/2) 
auf. Die zugehörigen Entwicklungen ergeben sich unmittelbar aus vorstehenden Gleichungen 
wenn ß,/2 statt £, gesetzt wird. 
. „Werden die so gewonnenen Ausdrücke in die Gl. (107), (108) und (109) eingesetzt. soemp- 
fiehlt sich die Ordnung von Zähler und Nenner nach Potenzen von H, wobei innerhalb des Koeffi- 
zienten jeder Potenz von H im Rahmen der beabsichtigten Linearisierung nur die niedrigste auf- 
tretende Potenz von £, zu berücksichtigen ist. Wird noch in Zähler und Nenner durch 16 (1—c) 
dividiert, so ergeben sich die nachstehenden Formeln: 


Zweite Linearisierung (S<1, Pf, <1, Pı < PB). 


HR? BA 
Isar Hbf tn 


—s 
H? 2 
gt 2Hßı; el 
Hierin ist zu setzen: | 
A.Für die dreischichtige Platte mit gegensymmetrischer Verformung: 


h= Em2B— 28, A=Cf2B—1, = @-4)Gin2ß+2L, 


5 
w=(2 —20) Ein2ß+BCo28—28, = @-20)C128—-1+REn2e | MM 
B. Für die dreischichtige Platte mit symmetrischer Verformung: 
f= Sin2ß+2%, f=&f2f+1, hA=8-4d)6in26—2%, 
5 
9=(5— 20) Ein24 + 6028428. H=ß—2)Cf2A +1 +BSm2R| 9 
C.Für die zweischichtige Platte: 
1 1 
h=zU 1 -5—AR, h=z &inap—2P, 
3 5 
n=(2-2)C418 +5 —-8c+1e+aP, BERN: le 


= (2 —2c) (Co AB —1)+B Sin 4B—4P?, 
= (7-20) Sin1a+ HC] aR 24 


Auf diese Formeln wird insbesondere im neunten Abschnitt Bezug genommen. 

Für die Diskussion der Ergebnisse ist außer den bisherigen Linearisierungen noch eine 
weitere sehr aufschlußreich, welche sich für den Bereich der langen Knickwellen durchführen 
läßt, d.h. für den eigentlichen ‚„‚Knickbereich‘“ (im Gegensatz zum „Beulbereich‘“, der das Gebiet 
der kurzen Wellen umfaßt). Im Knickbereich, der in den anfangserwähnten zweiten Euler bereich. 


a 
is 
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überführt, sind # und £, sehr klein gegenüber 1, da die Knickhalbwelle l als groß gegenüber A und 
h, angesehen werden kann. Daher können für diesen Bereich durch Entwicklung nach £ und £, 
vereinfachte Formeln erzielt werden. Da hierbei S<1, £<1, A, <1, jedoch nicht I, <£ 
vorauszusetzen ist, können die Beziehungen der zweiten Linearisierung nicht als Ausgangspunkt 
der Rechnung dienen, sondern es muß auf die erste Linearisierung zurückgegriffen werden. Zur 
Aufstellung der gewünschten Entwicklung ist es erforderlich, die Basisfunktionen (Gl. (105) und 
(106)) nach dem Argument zu entwickeln, wobei die Potenzreihen der Hyperbelfunktionen zur 
Anwendung kommen. 

Die Entwicklungen von ©; Diı Dzo Dar, Ds, und ©,, können aus den Gl. (110) entnommen 
werden. 

Für die übrigen Funktionen folgt: 


Du) R+ (1 ec) + De Du (= BB) + Fe) + en 
(116). 
Du) AB) + [a +, dA —EdR+B—F +: 


Durch Einsetzen in die Gl. (107), (108) und (109) ergeben sich nach kurzer Rechnung nach- 
stehende Formeln: 


Dritte Linearisierung (S<1, Pf, <1, P <]). 
A. Dreischichtige Platte mitgegensymmetrischer Verformung. 
h „ FRh 
1—H)®+HR+B)+ 
st) EEHA) FOR PER Er 


B. Dreischichtige Platte mit symmetrischer Verformung: 


1—2c 
HR. 0 ren ga 
SE En ee (118). 
2 )(+ HB)+T Heß 


C.Zweischichtige Platte: 


ihre ib 
3(1 —9(6+ 45) 


Eine weitere Vereinfachung erfahren die Formeln, wenn außer S<1, 9, <1,ß <1lauch I, <Pß 
vorausgesetzt wird. Diese Approximation sei als vierte Linearisierung bezeichnet; die zugehörigen 
Endformeln folgen entweder direkt aus den Gleichungen (117) bis (119), wenn innerhalb des 
Koeffizienten jeder Potenz von H nur die niedrigste Potenz von £, berücksichtigt wird, oder aus 
den Gleichungen (112) bis (115) durch Entwicklung nach £: 


Vierte Linearisierung (S<1, B<1, Ph,<ß). 
A.Dreischichtige Platte mit gegensymmetrischer Verformung: 
HH: ß ß4 
3 2 4 
31 —)P+HR)+6HPE 
B.Dreischichtige Platte mit symmetrischer Verformung: 


B+HB+ I HIER 
_ ga 2 a FRE ee (121) 
120 , 
@-)+HR)+ TE HM 


C. Zweischichtige Platte: 


—sS 


C 


BH2HRBH HB 


I gute . Se Ne ee 


Ö 


2 


a 


Sn 
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7%. Nachweis der E uler bereiche 


Wie schon anfangs erwähnt wurde, muß sich für Jange Knickhalbwellen, d. h. bei Knickung 
über die volle Plattenlänge (beiderseits frei drehbare Ränder vorausgesetzt) ein Verhalten der 
Platte herausstellen, daß dem eines schlanken Stabes ähnlich ist; d.h. der Bereich der langen 
Knickwellen muß sich als ‚„„E ule r‘‘-Bereich nachweisen lassen. 

Zur Unterscheidung vom selbständigen Knicken des Bleches bei sehr kurzen Knickwellen 
(„erster Euler -Bereich‘‘) wird jener Bereich als ‚zweiter Euler -Bereich‘ gekennzeichnet. 
Der Nachweis ergibt sich unmittelbar bei Anwendung der üblichen Näherungsvorstellung: Es 
wirken im wesentlichen nur Spannungen i,,, und die Querschnitte bleiben eben. Bei Anwendung 
der Formänderungsgleichungen (20) folgt mit t,, = 0: 


oryR.. 23 LER 15% 
days. 124200) 2): 
Definitionsgemäß war entsprechend Gl. (7) und (2) 
Ve Tor, 
ER NE ER ent 
x 0% 0y GEL 
U N oV, 
so daß nach Elimination en 
1—2coV, 
* 
e Tee 7 20.241259) 
zu setzen ist. 
Mit Bezug auf die erste der Gl: (20) gilt daher: 
Oz: E AR, 
Aehsek make ee EP eu (126). 


— n Be S2 
1+—\(1— 
| + „)‘ d) 
Andererseits folgt aus t,„=0 bei Integration nach y, wenn nunmehr die Verschiebung P, 
„Durchbiegung‘‘) als eine Funktion von z allein angesehen wird, 


dV. 
GV se ren ru dar) 


Demnach ist 

d:V 

bus = W (Y — Y) Fey . 
k 5: „) Ne) 


die Formänderungsgleichung, die für Biege- und Knickvorgänge bei Platten anzuwenden ist, 
deren Dicke als genügend klein gegenüber dem Krümmungsradius ihrer Biegelinie angesehen 
werden kann. Wird hierbei der Elastizitätsmodul E als Funktion von y aufgefaßt, so ist 
Gl. (128) für beliebig geschichtete Platten (für kontinuierliche und diskontinuierliche Schichtung) 
anwendbar. Bei Integration über die Plattendicke treten die Hilfsgrößen 


Yı Yı Yı 
J E(y)dy=K,, VE una 109) 


auf. Da beim Knickvorgang keine zusätzliche Axialkraft auftreten kann, muß das Integral der 
Spannungen über die Dicke verschwinden. Hieraus folgt eine Bedingung für die Integrations- 
konstante %9: 


8128) 


BR (130). 


Wird die auf die Längeneinheit der Plattenbreite bezogene Knicklast mit P bezeichnet, so hat 
das Moment der Knicklast, wenn ihre Wirkungslinie vor Eintreten des Knickvorganges den Ab- 
stand (\, von der Stelle y = 0 hatte, während der Knickung an jeder Stelle der Platte den Betrag 
P(C,—V,). Andererseits hat das Moment der Spannungen t,, den Betrag 


rn. Bahr EV .........caı) 


(1+)@-9 er” 
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Als Gleichgewichtsaussage ergibt sich demnach die Differentialgleichung: 
d?V, 
Bone + RB (V, per C;) —— 0 BY he, ae EN (132). 


Hierbei stellt die Größe 


die „‚Biegesteifigkeitdergeschichteten Platte‘ dar. Aus der Lösung der 
Differentialgleichung 
V=QGF+ Gsn (or) re (134) 
folgt die Knicklast 
DEOEBN LT se ee (135). 


Herrschte während des Aufbringens der Druckbelastung H o o kesches Gesetz mit derselben 
Verteilung der Elastizitätsmoduln, so läßt sich die Knicklast aus der Zusammendrückung — e 
nach der Gleichung 


Peam—maRK 2 A N (136) 
errechnen. Demnach folgt 
er 2 
eg a (137) 
0 Er e 2 
(1 — -) (1—c)K: 
und 
Bu K, R, Zu Ki 2 
RS l-oJR: DE N ae RE (138) 


Bei der dreischichtigen Platte wird 
Kh,=Eh-+2Eh, Kb; RK, ,[(E—E)b°+ E,(b+21,)] el‘)? 

_E-E)k®*+E,(h+2h) 
2([1+ 2) 0) 


_A—H)ß+H(ß+ 6) 
DASLS 3 (1540) (fe Day AR (141). 


Bei der zweischichtigen Platte wird 


B 


regen BORN. (140), 


RZ BRLRN R=-,(-BR+ Em), R,— (EB, M-+ EN). . (142), 
_BE(BE-B)M+ER(b+h+BH(H— HM 
= i 

12(1 lee (Eh+ Eh) 


_s-A-DAFALRH IM HHeH —1) (Br) 
31 —o)(+ HA): 
Bei Vergleich mit den Beziehungen (117) und (119) erkennt man völlige Übereinstimmung für 
die zweischichtige Platte. Bei der dreischichtigen Platte ist zu beachten, daß die Glieder mit H®? 
in Gl. (117) mit #4 (d.h. 9°) im Zähler bzw. 8? (d.h. 9%) im Nenner behaftet sind und daher bei 
langen Knickwellen, also kleinem Betrag p klein sind gegenüber £3 und (+ £,)? (d.h. 9°) im 
Zähler, bzw. £ und £, (d.h. 9)im Nenner. Beim Grenzübergang zum Euler -Verhalten stellen 
daher die Glieder mit 4? in Gl. (117) Größen dar, die mit abnehmenden g-Werten beliebig klein 
werden, so daß in der Tat Gl. (117) gegen Gl. (141) konvergiert. Damit ist der Grenzübergang 
zum zweiten Euler-Bereich nachgewiesen. 
Andererseits leisten die Gl. (117) und (119) für H—= & den Übergang zum ersten 
ae lerbereich (selbständiges Knicken also ‚„‚Beulen‘‘ der Außenschicht); in diesem Falle 
gilt: i ’ 


B 


3 


se 118 Sanıfrebr. 1953 Neuber, Theorie der Druckstabilität der Sandwichplatte II 17 


‘G. (141) des zweiten E ule r-Bereiches ist dagegen an das Ebenbleiben des Plattenquerschnittes 


gebunden und liefert mit H= o: 


ss PHBA+ AB 


1—c 


2 ,(146), 


also einen Ausdruck, welcher dem Knickvorgang zweier in festem Abstand unverschieblich 
gegeneinander festgehaltenen Bleche entspricht. 


8. Erweiterung auf nicht-lineares Formänderungsgesetz 


Bei guter Ausnutzung der Tragfähigkeit eines Bauteiles muß die Beanspruchung möglichst 
dem’ sicherheitsmäßig zulässigem Spannungswert genähert werden. Dieser liegt bei den meisten 
technischen Werkstoffen in einem Gebiet, in welchem die Formänderungskurve nicht mehr 
linear verläuft, so daß bei Stabilitätsbetrachtungen die bisherigen Ableitungen nicht mehr direkt 
anwendbar sind. Es läßt sich jedoch zeigen, daß die dargelegte Theorie bei Einführung gewisser 
vereinfachender Annahmen auch für nicht-lineares Formänderungsgesetz gültig bleibt. 

Wird zunächst ein Stab oder eine Platte von gleichem Werkstoff betrachtet, so führt die 
Frage nach der Stabilität eines Druckzustandes (Druckspannung — 6, Zusammendrückung — &) 
auf die Betrachtung der näheren Umgebung des zugehörigen Punktes A der Formänderungskurve 
(Bild 3). Wird die Formänderungskurve 
als stetig, differenzierbar und umkehrbar 
vorausgesetzt (von Verfestigung wird 
also abgesehen), so kann ihr Verlauf in 
der Umgebung von A durch die Gerade 


o=ö+(ce—) (2). (147) 


ersetzt werden. 

Die Spannung ö entspricht hierbei 
der- Vorspannung T',., während die be- 
liebig kleinen Spannungszuwüchse — 0 
mit den Spannungen t,,, die Dehnungs- 
zuwüchse e—£ mit den Dehnungen 


E R E > Bild 3. Zur Stabilitätsbetrachtung bei nicht-linearem 
oVx/ox identisch sind. An die Stelle Formänderungsgesetz 


des Elastizitätsmoduls tritt demnach der 

„fiktive Elastizitätsmodul‘ oder „Engessermodul‘”) während in Gl.(35) an die Stelle der 
Zusammendrückung — e eine fiktive Zusammendrückung zu setzen ist, welche mit — &* be- 
zeichnet werden möge; sie ist durch die Bedingung festgelegt, daß bei Multiplikation mit dem 
fiktiven Elastizitätsmodul die tatsächlich wirkende Druckvorspannung — 6 erhalten wird: 


ß | 
un = wi —— . . . . . . . . . . . . . . 
(02 & (2 2: (148). 
Hieraus folgt bei Division durch (do/de).=:: 
& a 
€ 57 0 Dee)! 


Eh 


Liegt — e* als Rechnungswert auf Grund der Theorie vor, so kann demnach sowohl — £, als 

auch —ö an Hand der Formänderungskurve mit Hilfe der Gl. (149) bestimmt werden. Die 

geometrische Betrachtung läßt erkennen (Bild 3), daß der zugehörige Punkt der Formänderungs- 

kurve der Bedingung entspricht, daß die Subtangente gleich — e* wird. Ein anderer Weg wäre 
z als Ordinate und — e als Ab- 

dino 

szisse. Der zugehörige Kurvenpunkt hätte hierbei die Ordinate — e*. 


die Konstruktion einer Kurve mit dem Funktionswert — 


Bei der mehrschichtigen Platte wurde zur einfacheren Durchführung der Rechnung an- 
genommen, daß alle Schichten bei Knickbeginn die gleiche Zusammendrückung — & besitzen. 
Für nicht-lineares Formänderungsgesetz bedeutet diese Annahme die Übereinstimmung von — &. 


?) Die Einführung der Größe do/de an Stelle des Elastizitätsmoduls bei Knickung von Stäben mit 
nicht-linearem Formänderungsgesetz geht auf F. Engesser zurück; vgl. Z. VDI 42 (1898), S. 927. 
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Wird ferner zur Vereinfachung angenommen, daß auch — e* übereinstimmt, so schneiden im 
Augenblick des Knickbeginnes die Tangenten der Formänderungskurven für die Werkstoffe der 
verschiedenen Schichten die e-Achse an der gleichen Stelle (Bild 4). Als Bedingungsgleichung 
folgt demnach: 


de de ) 
rer WeralR: ERNEST er a: 


Soll das Stabilitätsverhalten für alle Punkte der Formänderungskurve untersucht werden, 
so wird diese Bedingung zweckmäßig für alle Punkte vorgeschrieben, d.h. die Formänderungs- 
kurven unterscheiden sich nur durch 
einen Faktor: 

8 (E)= K onle)e na lan 
Demnach wird das Verhältnis Z der 
Elastizitätszahlen von e unabhängig. 
Ferner wurde angenommen, daß die 
Poissonschen Konstanten der 
Schichten untereinander übereinstim- 
men und von der Zusammendrückung 
unabhängig sind. 

Es ist zwar ohne weiteres möglich, 
die Theorie von diesen Einschränkun- 
gen frei zu machen; nur muß dann ein 
entsprechender Rechenaufwand in 
Kauf genommen werden. Es bleibt in- 
des fraglich, ob derartige Verfeinerun- 
gen noch praktisch interessieren; 
Bild 4. Zur Stabilitätsuntersuchung der mehrschichtigen Platte daher wurde in der vorliegenden Arbeit 

bei nicht-linearem Formänderungsgesetz hiervon abgesehen. 


9. Auswertung 


Für die praktische Anwendung kommen nur kleine Werte von S und £, in Betracht, also 
die Gl. (112) bis (115) der zweiten Linearisierung. Um das Stabilitätsverhalten einer bestimmten 
Ausführungsform mit festen Werten kA und A, beurteilen zu können, ist die Einführung des 
Wandstärkenverhältnisses «& zweckmäßig. Es wurde 


bei der zweischichtigen Platte: 


a=hjlh = ß,/2ß | 
a 2h/h=Bß | 


gesetzt. Die Gl. (112) bis (115) stellen daher bei festen Werten H und a die Größe — S und damit 
=— e* als Funktion von £ dar. Damit die Kurven näher zusammenrücken, wurde für die Auf- 
tragung die Größe ß, (= aß bzw. 2 «ß), als Abszisse gewählt, welche als ‚‚reduzierte Wandstärke“ 
bezeichnet werden kann. Die Auswertung wurde für breite Platten (c= 1/m) mit der 
Poissonschen Konstanten 1/m = 0,3 und einem Verhältnis der Elastizitätszahlen E,/E = 
H = 700 durchgeführt. Bild 5 und 6 zeigen die kritische Zusammendrückung der dreischichtigen 
Platte (gegensymmetrischer Fall) und der zweischichtigen Platte, wobei aus Gründen der Über- 
sichtlichkeit doppelt-logarithmischer Maßstab verwendet wurde. Von der Auftragung des sym- 
metrischen Knickfalles der dreischichtigen Platte konnte abgesehen werden, da die Rechnung 
im praktischen Bereich stets höhere Knicklasten liefert als im gegensymmetrischen Fall. 


bei der dreischichtigen Platte: + (182), 


Infolge der doppelt-logarithmischen Auftragung sind die beiden E u le r-Bereiche unmittel- 
bar am gradlinigen Verlauf der Kurven & = konst zu erkennen. Für & = 0,01; 0,15 und 0,10 sind 
in Bild 5 auch die der Arbeit von Gough-Elam-de Bruyne entsprechenden Kurven 
eingezeichnet. Sie entsprechen der Vernachlässigung von fo, 9, und g,. Für große £,-Werte ver- 
laufen sie noch in der Nähe der exakten Kurven, verlassen diese jedoch für abnehmendes £, mehr 
und mehr und streben für kleines 6, gegen den Festwert 1/2 Ha (1-+ 1/m), statt gegen Null, so daß 
sie den Übergang zum zweiten E u le r-Bereich nicht wiedergeben. 


Für die exakten Kurven & = konst ist außer dem klar hervortretenden Übergang zum 
zweiten Euler-Bereich die Ausbildung eines Minimums an der Stelle #,= 0,167, — e*= 0,0078 


O 
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= — :; von Bedeutung. Dies gilt bei der dreischichtigen Platte (Bild 5) für alle Kurven mit 
& S 0,06; bis etwa & = 0,08 findet eine kleine Verlagerung des Minimums statt; für «x > 0,08 
tritt schließlich kein Minimum mehr auf, und die Kurven verlaufen vom zweiten zum 
ersten Euler-Bereich in durchweg positiver Steigung. Bei der zweischichtigen Platte 
(Bild 6) zeigt sich das ausgeprägte Minimum bis etwa & = 0,04, während es für & > 0,05 
verschwunden ist. 


Da für die Ausbildung der Knickverformung der Platte als Knickhalbwellen / alle ganz- 
zahligen Bruchteile der Plattenlänge !* in Betracht kommen und andererseits nur der niedrigste 
— e*-Wert für das Eintreten der Verformung maßgeblich ist, so kommen die gestrichelt ge- 
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> Bild5. Kritische Stauchung der dreischichtigen Plattein Abhängigkeit von der reduzierten Wandstärke 


zeichneten Kurvenäste, welche oberhalb des Minimums — e} verlaufen, für die praktische An- 
wendung nicht in Betracht und sind durch die Stelle — e*= — ef, die kurz als „‚Beulgrenze‘ 
gekennzeichnet sei, ersetzt zu denken. 


Zur Erläuterung dieses Sachverhaltes wurde in Bild 7 die Zusammendrückung der drei- 
schichtigen Platte in gewöhnlichem Maßstab über ff = rh,/l* aufgetragen (l*— Plattenlänge). 
In dieser Darstellung kann die Beulgrenze an unendlich vielen, einzelnen Punkten der Geraden 
— e*= — e} liegen, je nachdem, wie groß die ganze Zahl der Knickhalbwellen ist, die sich über 
die Plattenlänge erstrecken. Die Einzelpunkte liegen im praktischen Bereich so nahe beieinander, 
daß die Gerade — e*= — e} direkt als Beulgrenze aufgefaßt werden kann (bei genauer Dar- 
stellung wären die Einzelpunkte girlandenartig durch Kurvenstücke verbunden; der so ent- 
stehende vielfach geknickte Kurvenzug fällt jedoch mit ausreichender Genauigkeit mit der 
Geraden zusammen). Die Linien « = konst steigen in diesem Diagramm zunächst bis zur Beul- 
grenze nahezu parabelförmig an. Unterhalb der Beulgrenze tritt stets Knicken über die gesamte 
Plattenlänge ein (n=/!*/i=1). Im Schnittpunkt A mit der Beulgrenze ist sowohl Knicken 
(n = 1), als auch Beulen (z. B. n = 348 für & = 0,005) möglich. Darüber hinaus tritt nur noch 
Beulen ein. 

2* 
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Bild 6. Kritische Stauchung der zweischichtigen Plattein Abhängigkeit von derreduzierten Wandstärke 
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Bild 7. Kritische Zusammendrückung der dreischichtigen Plattein gewöhnlichem Maßstab 
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10. Das optimale Wandstärkenverhältnis 


Zur Bestimmung der Tragfähigkeit der mehrschichtigen Platten war die Einführung einer 
dimensionslosen Hilfsgröße zweckmäßig, welche mit der Knicklast P in direktem Zusammen- 
hang steht. 

Ein solcher Kennwert ist die Größe 


p=—e(H+") N . (153). 


Wird die Plattenbreite mit b bezeichnet, so gilt für die Knicklast der dreischichtigen Platte 
(Index 3): 


P,=—e*(Eh+2E,M)b=—2eEh(H+L)b= bir En, er (154) 
bzw. — der zweischichtigen Platte (Index 2): 
P,=—e*(Eh+Eh)b=—e*Eh, (n4! )b=1 dl Ep er (155). 


Bei linearem Formänderungsgesetz sind hierbei alle Größen, außer P, p und e* Konstante, so 
daß P unmittelbar dem Kennwert p proportional ist. Bei nicht-linearem Formänderungsgesetz 
ist unter dem vereinfachenden Annahmen des Abschnittes 8 für E der Engessermodul do/de zu 
setzen, der seinerseits eine Funktion von e und damit von p darstellt. Daher wird P in diesem 
Falle eine nicht-lineare Funktion von p, und zur Beurteilung der Tragfähigkeit muß statt p das 
Produkt p do/de benutzt werden. In den Diagrammen 8 und 9 ist für y,= 16 y der Kennwert p 
über /#7 aufgetragen. Die zugehörigen Kurven & = konst erfahren an der Beulgrenze einen 
Knick und verlaufen darüber hinaus gradlinig. In die so gewonnenen Diagramme wurden zur 
Entscheidung der Frage des Gewichtoptimums Kurven eingezeichnet, längs denen ein für das 
- Baugewicht Q der Platte maßgeblicher dimensionsloser Kennwert 


elle ne) 
konstante Werte besitzt. Es folgt für das Baugewicht der dreischichtigen Platte (Index 3): 
EAN ENSE N AED DIRZAN- +) 2a A N 
bzw. der zweischichtigen Platte (Index 2): 
br (yhtyıl)= ya +) I deal 


Demnach ist Q direkt proportional g, und alle Punkte einer Kurve qg = konst entsprechen den 
Abmessungen gewichtsgleicher Platten. Die günstigste Ausführung ist mithin jene, zu welcher 
auf der zugehörigen Kurve q = konst der höchste Wert von P, d.h. beilinearem Formänderungs- 
gesetz der höchste Wert von p gehört. Die Maxima der Kurven q = konst stellen daher die 
optimalen Ausführungsarten geschichteter Platten mit linearem Formänderungsgesetz dar. Die 
Verbindungslinie der Maxima wurde mit der Bezeichnung ‚‚Optimum‘‘ versehen. Sie verläuft 
bei der dreischichtigen Platte (Bild 8) im zweiten Euler-Bereich gradlinig. Für größere ff- 
Werte erreichen die Kurven g = konst die Beulgrenze und erfahren dort einen deutlich erkenn- 
baren Knick, derart, daß die so entstehende Spitze zugleich dem Höchstwert und damit der je- 
wejligen Lage des Optimums entspricht. Die Optimumlinie fällt daher nach Erreichen der Beul- 
grenze mit dieser selbst zusammen. 

Im zweiten Euler-Bereich läßt sich das optimale Wandstärkenverhältnis an Hand der 
Ausführungen des siebenten Abschnittes leicht herleiten. 

Im zweiten Euler-Bereich der dreischichtigen Platte (Index 3) gilt nämlich wegen 
ßı= £T mit Bezug auf Gl. (141): 


year) 
TE er nee ana nu. .5 (159) 

(1 le —c)(1+ Ha)o? 
Entsprechend Gl. (153) wird 
= . 3 
Krmlen+s 


ea 23. a (160). 


ps; = 
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Andererseits gilt für den Gewichtskennwert mit Bezug auf Gl. (156) 


Mithin gilt 
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Bild 8. Reduzierte Knicklast der Gewichtskurven der dreischichtigen Platte 
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Wird diese Beziehung zur Elimination von #7 in Gl. (160) benutzt, so folgt für die dreischichtige 
Platte 


3 ser SıE: ea RR = (163). 


rannte) 


Bei konstanten g, ist demnach p, eine Funktion von «a allein. Das Maximum entspricht dem Ver- 
schwinden der Ableitung nach a: 


(1+% oJ) +24) |; +Hla++° 3) = ER: (164). 
Die Auflösung führt auf die quadratische Gleichung 
vi 
EN A ee ET fi De Pe (165) 
ai a 
Y 
mit der positiven Wurzel 
a 1) e £ # 
= —1+ ve N, bzw. für 2 >|, yı>2y: %=,— (166). 
Ay — y) 2 yı 


Wie in Bild 8 ersichtlich ist, entspricht die für — 26 eingezeichnete Optimumlinie im zweiten 


Euler-Bereich sehr genau dieser Formel. 
Der Kennwert p, erreicht für das optimale & den Beträz 


H—-Dyg 
© TREE Tr 16 
P3 3(1+ Nr ntn dr! (167), 
HA yı —y 
bzw. für H>1 und y,>y: 
3 
Ps = NEE (168). 


1 
7(1+)u-9n 
Für die zweischichtige Platte (Index 2) gilt im zweiten E ule r-Bereich entsprechende Gl. (144) 

8 4 ER 
Be en ee Fe) HE Ne VER ER (169). 
12[1 + Ja —c)(1+ Ho)? a? 
Mit Bezug auf Gl. (153) wird 
MID HH Ve! 
121 + )a BIN. 


Da für den Gewichtskennwert wieder Gl. (156) gilt, kann zur Elimination von #7 Gl. (159) be- 
nutzt werden. Dann ergibt sich 


I-HHFAUFe LH N) 5 
1214) -90+Hajlı +2.) ; 
m Y 


Durch Nullsetzen der Ableitung nach « folgt unter Weglassung der Zwischenrechnung eine sym- 
metrische Gleichung vierten Grades für aVH: 


(2%) [ea +1] + aVa(ı 7) [(erI°+ av] 


% +2|2(1-=2)+ 2-2. (m =0 


n= re (170) 


Pr = 


7) 
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Da H sehr groß gegenüber 1 ist, wird die gesuchte Wurzel zweckmäßig in Form einer Entwicklung 
nach H dargestellt. Der Aufbau der Gl. (172) führt auf folgende Entwicklung: 


1 1 1 
= l —— DD Vier i L  n (173) 
% VE He 
Durch Einsetzen in Gl. (172) folgt für a? die quadratische Gleichung: 
«+2(1-2%)e+1=0 nr 2172) 
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mit den positiv-reellen Wurzeln 


ER - r-ı +y (2-1) -1 KT (175). 
y y 


Der maximalen Knicklast entspricht das untere Vorzeichen vor der zweiten Wurzel. Für H>1 
und y,>y folgt mithin als optimales Wandstärkenverhältnis der zweischichtigen Platte im zweiten 
Euler bereich 


a ee 

Er Em nn mu mn nn m a e 

IFFSEEHHS RS ER 

a FE EI IE I EZ Era BREI EEE I FR 

ee eh 

22 ER ERDE RE RE 

ee SE] 

„11290 SCÄRE Ze Bez a IA BI DE Er BZ EEE RA 

Pe || 

Er ar 12 u er ee Anz FRE: THE ERRE BREET Dr GI HEN 

Peer 

: ee 

NE EEK ea TS 

a RE RL HL 
IN 


e 
+ 
n 


SA NaL. 


SOSE HT HEN 
AN FNA- 


BE EmER271071177////ANNENN 
BEE ER RR ER ALT TB ZZIN ARE 


BEEsa 


BE 


VEN FRHE 


P Ta a ANAND @ 
: Io orNdneE = 
y a a By IP AU UN AN 70 Var 

3 En Erna, ze 


IH—— -— ‘ 


P =Knıcklast p® Zr 


«Ah, EEE, 


Diagramm für: u. =700, c= 4=03 f 
Bert N le TLoiT 
ARRRA 7371 2A RR RR PER IN 


RAS li2 3 56 Bar: 2 3 4.567890" 


Bild 10. Optimumlinien der dreischichtigen Platte für verschiedene Gewichtsverhältnisse 
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Der zugehörige p-Wert ergibt sich aus Gl. (171). Bei Beschränkung auf den Sonderfall > 1, 
yı>y folgt nach kurzer Rechnung 


CN 9 . (177). 


Die Optimumlinie für die zweischichtige Platte wurde im Bild 9für den zweiten Eulerbereich nicht 
mehr eingezeichnet, da sich herausstellte, daß die Beziehungen der zweiten Linearisierung für ihre 
exakte Ermittlung nicht ausreichen. Die im Bereich höherer #7-Werte eingezeichneten p-Kurven 
zeigen wieder ausgeprägte Spitzen auf der Beulgrenze, die also hier wieder mit der Optimumlinie 
zusammenfällt. Bei abnehmenden #7 werden die Spitzen flacher, sodaß eine Abzweigung der 
Optimumlinie im zweiten Eulerbereich im Sinne der Gl. (176) und (177) als wahrscheinlich 
anzusehen ist. 

Die Berechnungen der Optimalausführungen im zweiten Eulerbereich erlauben un- 
mittelbar einen Vergleich der drei- und zweischichtigen Platte. Mit Hilfe der Beziehungen 
(154), (155), (157), (158), (168) und (177) folgt für das Verhältnis der Knicklasten bei gleichem 
Baugewicht (9, = Q;) unter Umgehung der Er 

BB TER ee (178). 
'P; are salz, ae BI 


Für 4=700 und y},=16y ist demnach die HL ENT. der dreischichtigen Platte im zweiten 
Euler bereich etwa 9mal größer als die der zweischichtigen Platte. Im ersten Eulerbereich 
sind infolge des selbständigen Knickens des Bleches beide Ausführungen gleichwertig. Dieser 
Grenzfall wird jedoch in praktischen Fällen nicht erreicht. Im Zwischengebiet zeigen die Dia- 
gramme 8 und 9 stets die Überlegenheit der dreischichtigen Platte. Bild 10 zeigt die Optimum- 
linien für verschiedene Gewichtsverhältnisse. 

Für nicht-lineares Formveränderungsgesetz sind die Diagramme 8, 9und 10 für das Produkt 
p do/de unter Berücksichtigung des mit e veränderlichen En gesser moduls mit Bezug auf die 
genaue Formveränderungskurve des Werkstoffes aus Bild 5 und 6 neu zu entwickeln, wobei sich 
eine Verlagerung der Kurven g=Konst und damit der Optimumlinien ergeben kann. 


Eingegangen: 24. November 1951. 


Über die Fehlerabschätzung beim Extrapolationsverfahren 
von Adams 


II. Gleichungen zweiter und höherer Ordnung 
Von P. Matthieu in Zürich/Berlin 


Die vom Verfasser kürzlich entwickelte neue Methode der Fehlerabschätzung wird auf Gleichungen höherer 
Ordnung ausgedehnt. 


Thenew method of estimating the error given recently by the author is extended to equations of higher orders,. 


La nouvelle methode de l’estimation des erreurs d&veloppee r&cemment par l’auteur est &tendue sur des 
öquations d’ordre superieur. 


Passat aBTOPOM HeNaBHO HOBbIÜ METOA OUeHKH OMIHÖOK PacıpocTpaHneTcH Ha YPpaBHeHHA 
BbIcHHIeTO NHOPAAKA. 


1. Einleitung 

In einer kürzlich erschienenen ersten Arbeit!), die hier als bekannt vorausgesetzt und im 
folgenden mit „I“ bezeichnet wird, wurde eine neue Methode der Fehlerabschätzung für das Ver- 
fahren von Adams zur Integration von Gleichungen erster, Ordnung entwickelt. In der vorlie- 
genden zweiten Arbeit (mit „II“ bezeichnet) wird gezeigt, daß sich diese Methode ohne wesent- 
liche Änderungen auch auf die Integration sämtlicher Gleichungen höherer Ordnung ausdehnen 
läßt, ein Problem, das bis jetzt nie systematisch in Angriff genommen wurde. Es ergibt sich auf 
diese Weise eine sehr allgemeine und geschlossene Näherungstheorie für das Adamssche Verfah- 
ren zur Integration beliebiger gewöhnlicher Differentialgleichungen. 

Der in Ieingeschlagene Weg beruht im wesentlichen auf zwei wichtigen Prinzipien. Erstens 
wird zu einer gegebenen Gleichung erster Ordnung die zugehörige Fehlergleichung (1,21) auf- 


1) Matthieu, P.: Über die Fehlerabschätzung beim Extrapolationsverfahren von Adams, Z. angew. 
Math. Mech., Bd. 31 (1951), S. 356—370. 
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gestellt, deren Lösung mit größter Näherung den Fehler darstellt. Zweitens kann die wesentliche 
Funktion g(x)in der Fehlergleichung (sie wird im folgenden als die infinitesimale Funktion der 
Fehlergleichung bezeichnet, da sie die Größenordnung der Differenzen Art! f# mitn 2zm-+1 
hat, wenn zur Interpolation Polynome vom Grade m verwendet werden) stets mit größter Leichtig- 
keit aus dem Differenzenschema abgeschätzt werden, wodurch sich dann sofort eine Abschätzung 
für den Fehler ergibt. Als weiteres wichtiges Resultat folgt aus diesen Betrachtungen ferner, daß 
sich der zunächst viel komplizierter scheinende Fall der mittelbaren Extrapolation ganz leicht auf 
den der unmittelbaren zurückführen läßt. 


Alle diese Gedankengänge werden jetzt unmittelbar auf Gleichungen höherer Ordnung 
übertragen. Zunächst wird auf Grund der Methoden der Störungsrechnung wieder die zugehörige 
Fehlergleichung aufgestellt. In dieser tritt auch jetzt eine infinitesimale Funktion auf, welche mit 
größter Leichtigkeit aus dem Differenzenschema abgeschätzt werden kann, wodurch sich dann im 
Prinzip wieder eine Abschätzung für den Fehler ergibt. Weiterhin ergeben sich auch bei Gleichun- 
gen höherer Ordnung ganz einfache Beziehungen zwischen unmittelbarer und mittelbarer Extra- 
polation. Der Hauptunterschied gegenüber dem in I behandelten Fall liegt darin, daß jetzt die 
Fehlergleichung nicht mehr eine lineare Gleichung erster Ordnung ist, welche in geschlossener 
Form integriert werden kann, daß sie vielmehr linear von höherer Ordnung ist. 


Obwohl die vorliegenden Betrachtungen ohne jegliche Schwierigkeit für den Fall durch- 
geführt werden können, daß Extrapolationspolynome beliebigen Grades verwendet werden, 
beschränken wir uns der Einfachheit halber wie in I auf Polynome dritten Grades. Wir gehen nun 
daran, die verschiedenen hierher gehörenden Fälle einzeln zu besprechen. 


2. Gleichungen zweiter Ordnung bei unmittelbarer Extrapolation 
Es sei eine reelle Differentialgleichung zweiter Ordnung in der Gestalt 
ur TIL, Yy,Y') NE N Rn er (1) 
vorgelegt mit den Randbedingungen 
Yy(%) = Yo» Ylo)=Y m: rennen l2) 


Bezeichnet man wieder mit y* (x)irgendeine Näherungslösung, welche nur den Randbedingungen 
(2) genügen soll und mit 


den Fehler, so kann man genau nach dem gleichen Verfahren wie in I die Fehlergleichung auf- 
stellen, die jetzt 


0 
"=, y", v3 (a YES U U Ye yeT 3, 2 (4) 
lautet oder bei Einführung geeigneter Abkürzungen die (1,21) entsprechende Gestalt 
ö" BO ee kernel een (0) 
annimmt. Hier hat die infinitesimale Funktion die Form 
U EN N Ah) a. (6) 
Diejenige Lösung von (5), die den Randbedingungen 
ö(2,) = 0 Et: (7) 


genügt, stellt auf Grund von (2) und (3) wieder mit großer Annäherung den Fehler dar, und es ist 
auch leicht, aus der Gestalt der gegebenen Gl. (1) die Zulässigkeit dieser Näherung abzuschätzen. 
Praktisch wird sie immer erlaubt sein. 


Bei der Anwendung dieser allgemeinen Prinzipien auf das Verfahren von Adams ist wieder 
zu berücksichtigen, daß sich die Näherung %* (x) aus einzelnen Teilen zusammensetzt, die fort- 
laufend für die verschiedenen Intervalle gebildet werden. Die Bildung der einzelnen Teile erfolgt 
ganz entsprechend wie in I durch Verwendung des Differenzenschemas?) 


2) Vgl.die folgenden Werke: Willers, Fr. A.: Methoden der praktischen Analysis, 2. Aufl. Berlin 
1950, S. 372ff. — Sanden, H. v.: Praxis der Differentialgleichungen, 3. Aufl. Berlin 1945, S.48ff. — 
Collatz, L.: Numerische Behandlung von Differentialgleichungen, Berlin 1951, S. 68—86. — Collatz, L.: 
Differenzenverfahren zur numerischen Lösung von gewöhnlichen Differentialgleichungen n-ter Ordnung. 
Z.angew. Math. Mech., Bd. 29. 


a a ca a 
a I { x 
f K. s 


28 Matthieu, Fehlerabschätzung beim Extrapolationsverfahren von Adams „4 33 a 


*/ x * 
Top Yo Yo f 


Aft 
nun v je Af$ 
Aff Asft 
mn Won PR af} A 
Aft Af, 8), 
zu je At a3 
Aft Asfs 
DB LE. Aeft 
Äfz | 3 


N RT TEN 
von dem wieder die ersten Werte vorausberechnet werden müssen. Wie in I machen wir die Vor- 
aussetzung (von der man sich nötigenfalls auch sofort befreien kann), daß diese Anfangswerte 
nicht merklich von der exakten Lösung abweichen. An Stelle von (1,25) bis (1,28) treten jetzt 
unter Verwendung entsprechender Bezeichnungen für n>3 die Gleichungen ?) 


h ‚ Afe- Aeft_ 
ade) = Baıla)esye st % et Le EN 


ne ey | 
a Ö 
Ye) = at hf am at A 3 Bf + 2A) 
DE Sea. ie (10), 
FR A) St + Aue a] 
Vila) = ya that Ri 


5 
Un | 


+ (Aa + A) 55 + Aa Ion, 


’ 7 Y 1 FR = FE 3 2. 
ws alfta tt At + At girl - - 2) 


Y=Vr(m)= yrı + hy-ı 
* 188 a 19 — 13), 
a r BR ER 6 Afn-ı tr g Afn-ıHt 180 Adfı-ı (13) 
von denen zur Fortsetzung des Schemas (8) nur die beiden letzten gebraucht werden, während die 
drei ersten für die Fehlerabschätzung von Bedeutung sind. : 
Aus (9) bis (11) ergibt sich weiterhin gemäß (6) sofort die infinitesimale Funktion g(x) im 
n-ten Intervall für n>3 (für n<3 verschwindet g(x) nach den gemachten Voraussetzungen 
identisch) in der Form 


a(@)=f(e, Yil), Yi(2)) — Fr (a) ei: 
=f(x, Yi(a), Y# (2)) —f (8, y(a), y()) — [YR" (8) — y’(@)l 


für die man leicht entsprechend (1,32) und (1,33) im n-ten Intervall die Beziehungen 


gm )=0 gla,) = AS Re N ee (15) 
herleitet. Aber auch die Schlüsse, die in I gezogen wurden betreffend den monoton ansteigenden 
und nach oben konkaven Verlauf von g(x) in jedem Intervall können unverändert übernommen 
werden, denn der wesentliche Term in (14) ist wiederum der eckige Klammerausdruck, der von 
vierter Ordnung in hist, während der übrig bleibende Ausdruck die 5. Ordnung hat und sich sogar 
von dem entsprechenden früheren Ausdruck nur um Glieder 6. Ordnung unterscheidet. Zusam- 
menfassend kann also gesagt werden, daß g(x) auch für Gleichungen höherer Ordnung genau 
dieselben Eigenschaften hat wie sie in I für Gleichungen erster Ordnung hergeleitet wurden. Es 
verläuft zackenförmig vom Werte Nullam Anfang des n-ten Intervalls (n>3)bis zum Werte A@f, 
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an dessen Ende (vgl. Bild 1), und der Verlauf wird im wesentlichen, d.h. bis auf Glieder höherer 
Ordnung gegeben durch 


ENT Ei 1 BR | 
g(&) 27, ee mt 94 nt 4 Un-ı + 34 Un—1 ) miascesm Iswmisl (16). 


Auch die Möglichkeiten der strengen Kontrolle für den Verlauf von (16) übertragen sich 
unmittelbar auf den jetzigen Fall. Entsprechend den Gleichungen (1,47) bis (1,49) leitet man jetzt 
aus (9) bis (11) z.B. für die Mitte des n-ter Intervalles (n > 3) sofort die Formeln ab 


S EEE 10% IN ;e 5’ TE RB = 
Yy_ı2 = mit > Y—ı-+ RR ER: r 18 Afn-ı + 384 Afnı + 11550 A, | (17), 


r 1 1 - 10 DR 
rat ats tat... .....08, 
pt At Aa A a9) 
n—1/2 FeN. 5 nn 8 n-1 16 ls: 4 A Mer ArFnE Kine Nr s) re & 
An—1/2 2. f (n—1/2 A Yn—ıj2 , Yn—ıj2) —— Yn-1/2 NE er a RE LAT < (20). 


Die Formeln (18) und (19) stimmen im wesentlichen mit (1,47) und (1,48) überein. Andere 
Zwischenwerte ergeben sich ganz entsprechend aus (9) bis (11). 
Es handelt sich nun weiterhin darum, auf Grund der eben hergeleiteten Eigenschaften von 
7 q(x) gemäß der Fehlergleichung (5) eine Abschätzung für den Fehler zu finden. Hierbei bieten 
sich vor allem zwei Möglichkeiten dar. Man kann erstens von der allgemeinen Lösung der zu (5) 
gehörenden homogenen Gleichung ausgehen, die man in manchen Fällen in geschlossener Form, 
stets jedoch approximativ bestimmen kann, wobei im letzteren Fall sehr rohe und einfache 
Methoden, z.B. graphische®), genügen. Sind dann ö, (x) und ö,(z) irgend zwei linear unabhängige 
Lösungen dieser homogenen Gleichung, so stellt sich der Fehler ö(x) als Lösung von (5) und (7) 
vom vierten Intervall an (in den drei ersten verschwindet ö(x) nach Voraussetzung identisch) 
bekanntlich in der Form dar*) 


a aelae) 20: 2 Tkate)are) - 
ö(x) == öı (x) | 8,6, en; ö, de sa | Be dE N - (21). 
Für den Nenner in den Integranden gilt dabei die Identität 
& 
n Sr(ea 
Are nn _ AD) _ 


wobei C’eine Konstante ist, die von ö, (x) und ö,(x) abhängt und aus (22) leicht bestimmt werden 
kann, z.B. indem man £ —=z, setzt. Führt man noch die Abkürzung ein 


61 (X) d,(2) 


VE SE ES ns CE RE RE RENENER 23), 
| EU EROENG es 
br) so läßt sich (21) auf Grund von (22) und (23) auch in der übersichtlichen Gestalt 

a 
1 a de 
ö(z) = 6 Ir ; RE Rp I (24) 


da 
schreiben. Diese Form hat überdies den Vorteil, daß sie sich unverändert auch auf Gleichungen 
höherer Ordnung übertragen läßt (vgl. dazu Abschnitt 5). 
Aus (24) folgt nun sofort die (1,51) entsprechende Abschätzung 


= 


5 EYE 1 [ze RE ; | ze) 
[ö(@)) < =, 1(8) e | (2, )| Di 5,5] 


D(z, E)|dE (25). 


3) Vgl. die folgenden Werke: Runge, C. und König, H.: Vorlesungen über numerisches Rechnen, 
1. Aufl. Berlin 1924, 10. Kapitel. — Sanden, H. v.: Praxis der Differentialgleichungen, 3.,Aufl. Berlin 1945, 
S.41—48. — Runge, C.: Graphische Methoden, Leipzig und Berlin 1915, $16. — Hütte: 27. Aufl., 
8. 203—205. £ 

4) Vgl.z.B. Frank, P. und Mises, R. v.: Die Differential- und Integralgleichungen der Mechanik 
und Physik, 2. Aufl. 1930, S.298—300. Ferner Forsyth, A. R.: Lehrbuch der Differentialgleichungen, 
2. Aufl. Braunschweig 1912, S. 119ff. 
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Wählt man die mittlere Majorante 9(x) wie in I, S. 364 so, daß sie ungefähr in halber Höhe der 
Zackenspitzen von q(x) verläuft (Bild 1)und wechseln D(z, £) und g(x) im ganzen Integrations- 
intervall das Vorzeichen nicht, so folgt wieder aus (16), daß die Fehlerschranke stets etwa das 
Anderthalbfache des Betrages des wirklichen Fehlers ist, während sich das Vorzeichen aus (24) 
und (22) ergibt. 

Im folgenden erweist es sich als bequem, neben der mittleren Majorante 9(z) für Ia(z)| noch 
in genau analoger Weise eine mittlere Majorante 9,, (2) für q(x) einzuführen, ebenso eine mittlere 
Minorante q„(x). Falls g(x) Z0, kann z.B. qg,(2)=9(x) und q,(x)=0 gesetzt werden, bei nega- 
tivem g(x) dagegen qy,(x2)=0 und 9,(2)= —q(R). Es sind aber auch engere derartige Schranken 
denkbar. Setzt man 9y(x) an Stelle von qg(x) in die Fehlergleichung (5) ein, so erhält man das 
Gleichungssystem 


E Or = rösr+ Po + Ur, da (m) du lm) 0 ee eh 
dessen Lösung ö,, (2) gemäß (24) in der Gestalt 


T 


f LT 
due) = 1 [zo “N Dim Das 


G 
dargestellt werden kann. In genau gleicher Weise kann man auch 9, (x) in (5) einsetzen, wodurch 
eine Lösung ö,„(x) definiert wird. Falls die Funktion 


Da 3) bzw. ODE) NR 


bei beliebigen Werten von x und £ innerhalb des Integrationsintervalls stets positiv bleibt, gilt 
dann, wie man aus (24) unmittelbar erkennt, die Fehlerabschätzung 


dla) älr)< du le) Re er a 


Wenn umgekehrt (28) stets negativ ist, sind in (29) die beiden Schranken zu vertauschen. 

In dem eben genannten Falle, daß D(x, £) im Integrationsintervall stets das gleiche Vor- 
zeichen hat, können die Formeln (25) und (29) unmittelbar zur Fehlerabschätzung verwendet 
werden. Bei Kenntnis von ö,(x) und ö,(x) kann das Vorzeichen von D(z, &) direkt festgestellt 
werden. Weiterhin ist (28) nach einem nächstens (S. 30—31) auszusprechenden Satz stets positiv, 
wenn p(z) im Intergrationsintervall überall =0 ist. Bei der Abschätzung wird man sich natür- 
lich wieder geeigneter, z. B. graphischer Methoden bedienen?). Natürlich kann man auch den 
Integranden oder Teile von ihm durch irgendwelche Majoranten ersetzen, wodurch man zu sehr 
einfachen Abschätzungen gelangen kann. Oft wird man auch die rechte Seite von (25) ent- 
sprechend (21) wieder in zwei Summanden zerlegen und mit diesen einzeln rechnen. 

Wenn dagegen D(z, €) das Vorzeichen wechselt, braucht (29) bzw. die entsprechende Formel 
bei negativem Wert (28) nicht mehr richtig zu sein, während (25) zwar noch richtig bleibt, aber 
nicht mehr einfach zu berechnen ist. Man greift dann am besten auf (21) bzw. (24) zurück, aber 
in diesen Gleichungen führt jetzt die Ersetzung von q(z) durch eine Majorante oder Minorante 
nicht mehr notwendigerweise zu einer Majorante oder Minorante für ö(x). Man kann aber in 
diesem Fall dadurch zum Ziel kommen, daß man in (21) und (24) gleichzeitig mit Majoranten 
oder Minoranten arbeitet. Die leichten Überlegungen führen in Verbindung mit den bisherigen 
Resultaten zu folgender Regel: 

Wenn D(z,£) im ganzen Integrationsintervall stets das gleiche Vorzeichen hat, benutzt man für 
die Fehlerabschätzung am besten die Formeln (25) und (29). Andernfalls erhält man dadurch eine 
Majorante für ö(x), daß man in (21) bzw. (24) überall dort, wo die Ausdrücke O' d,(z) d,(E)g(E) und 
—Cö,(a) d1(E) q(E) positiv bzw. negativ sind, g(&) durch Iu(&) bzw. Qm($) ersetzt. In ganz 
entsprechender Weise kann man auch eine Minorante für ö(x) herleiten. 

Die zweite Methode der Fehlerabschätzung beruht darauf, daß man die Fehlergleichung 
direkt löst oder ihre Lösung abschätzt. Falls p(x) im Integrationsintervall positiv ist (der Fall 
eines negativen p(x) wird im 6. Abschnitt behandelt), stützt man sich dabei am besten auf das 
folgende Theorem: 

Es sei eine reelle Differentialgleichung zweiter Ordnung mit Anfangsbedingungen in der Gestalt 


vv thR@ytol@), ya)=c, Yi)=d ...... (30) 


gegeben, und es sei ferner im Intervall asz<b fs(z)20. Wenn man dann zu (30) ein Glei- 
chungssystem der Form 


= f(J"+fhle)Y+C(R), Y.a)=(, Yıla)=3D ar 2, (31) 


5) Außerdem kann man ö (x) auch wieder in tabellarischer Form berechnen wie in I, S. 366. Es soll dies 
hier jedoch nicht ausgeführt werden. 


C 


Me 
Pr 
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bildet mit G(x)Z g(e), ic, DZ ERS. Baar (32), 
so ist die Lösung Y (x) eine Majorante der früheren Lösung y(x). Es gilt also 
Y(x)>2 y(x), NEAR 1) (33). 


In ganz entsprechender Weise läßt sich auch eine Minorante für die Lösung herleiten, indem man g(x) 
durch eine Minorante ersetzt und die Anfangswerte verkleinert. 

Die Richtigkeit dieses Satzes erkennt man am besten, wenn man (30) durch das lineare 
Gleichungssystem 

"= fıla)z+Rla)y+ ge), Era. ala) det (3A) 
ersetzt, die entsprechende Umformung für (31) vornimmt und den Verlauf der Funktionen %(z) 
und 2 (x) mit demjenigen der entsprechenden Funktionen Y (x) und Z (x) vergleicht. Man erkennt 
dann, daß es nicht möglich ist, daß Z(x) irgendwo kleiner wird als z(z). Denn wenn die Kurve Z (x) 
an einer Stelle die Kurve z(&) erreicht, so folgt aus (34) und der entsprechenden Umformung von 
(31) sofort, daß sie sich von da an wieder über 2(z) erheben muß oder im Grenzfall mit ihr zu- 
sammenfällt. Für den Fall einer Minorante verläuft der Beweis genau analog. Dieses wichtige 
Theorem läßt sich ohne Schwierigkeiten auch auf Gleichungen höherer Ordnung übertragen. Es 
muß dann ‚vorausgesetzt werden, daß alle Koeffizienten mit Ausnahme desjenigen des zweit- 
höchsten Gliedes > 0 sind. 

Die Anwendung auf unseren Fall ist unmittelbar. Ersetzen wir die Fehlergleichung (5) 
durch (26), so stellt nach dem obigen Satz die Lösung öy, (x) eine Majorante für den Fehler dar. 
In gleicher Weise bestimmt sich auch die Minorante Öm(z), so daß man den Verlauf des Fehlers 
gemäß (29) beherrscht. Bei dieser zweiten Methode wird jetzt also. öyr(z) nicht mehr gemäß (27) 
berechnet, sondern durch direkte Lösung von (26) bestimmt, wobei wieder die einfachsten und 
primitivsten Methoden Verwendung finden können?). Entsprechendes gilt natürlich für d,,(8). 
Zusammenfassend kann man also sagen: 


Ist p(x) im Integrationsintervall überall positiv, so stellen die Lösungen dur () und ön( x) ge- 
mäß (29) stets Schranken für den Fehler dar. 


3. Beispiel zur Fehlerabschätzung bei unmittelbarer Extrapolation 
Die Gleichung y"+2y’—+ y=0 mit den Anfangsbedingungen y(0)=0 und y’(0)=1 soll 
im Intervall O<Sxz<s1 nach Adams integriert werden, und es soll der dabei begangene Fehler 
abgeschätzt werden. 
Die Lösung kann für die Anfangswerte wieder sehr rasch und genau durch Potenzreihen 
berechnet werden. Das Schema (8) nimmt dann, in analoger Weise wie bei dem in I, Abschnitt 5 
berechneten Beispiel, die Gestalt (35) an. Ferner wurde, obwohl das für die Fehlerabschätzung 


z, ys yr fi 107 Afs  107A2f% 10rASft 107 AfE 1 

0,1 —0,1105171 1,2156880 —2,3208589 
3208589 

0,0  0,0000000 1,0000000 —2,0000000 —400500 
2808089 47168 

0,1  0,0904837 0,8143537 —1,7191911 — 353332 5729 
2454757 41439 

0,2 0,1637462 0,6549846 —1,4737154 —311893 —3801 
2142864 37638 

0,3  0,2222460 0,5185915 —1,2594290 — 274255 —5086 
1868609 32552 

0,4  0,2681305 0,4022188 —1,0725681 —241703 BASE 1..(85) 
1626906 .. 29309 

0,5  0,3032707 0,3033034 —-0,9098775 — 212394 3791 
1414512 25518 

0,6  0,3292960 0,2195652 —0,7684263 —186876 —2634 
1227636 29884 

0,7. 0,3476227 0,1490200 —0,6456627 —163992 —2890 
1063644 19994 

0,5 0,3594803 0,0899090 —0,5392983 —143998 — 2146 
919646 17848 

0,9 0,3659339 0,0406999 —0,4473337 —126150 
793496 | 


1,0  0,3679046 0,0000397 —0,3679841 


TE ES RE N Me 
x > x as 2: 
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nicht notwendig ist, die Funktion |q(x)| berechnet. Sie ist, gleichzeitig mit (x) und einer ge- 
eigneten Funktion 9(z), in Bild 1 dargestellt. Die Berechnung ist dort sogar bis & = 2,0 weiter- 
geführt. Man erkennt daraus durch 
Vergleich mit I, Bild1 den völlig 
analogen Verlauf der beiden Funk- 
tionen g(x), wie er im vorigen 
-IHLX)I=-g0%) Abschnitt theoretisch erschlossen 
ee) wurde. 
—— 1 — = 50)+032:10°? e°* Aus dem Schema (35) kann 
nun weiterhin. der Fehler abge- 
schätzt werden. Es soll dies, da 
einerseits die Fehlergleichung 


6" = —20'—Ö i 
BEN 
6(0,2) = 6’ (0,2) = 0 
in geschlossener Form lösbar ist 
und andererseits p(x) negativ ist, 
nach der ersten der beiden an- 


gegebenen Methoden geschehen. Die zu (36) gehörende homogene Gleichung hat die beiden 
partikulären Integrale 


Bild 1 


lal=ie*, Öle) >= ee REN ee DR el; 
Daraus folgt gemäß (23) die Beziehung 
Dix, &$)= —(2 —!) a EN IE A (38). 


Demgemäß hat D(z, &) im ganzen Integrationsintervall das negative Vorzeichen, und es kann also 
die Abschätzung (25) verwendet werden. Man erhält gemäß (25) 


Si) = [UE—Herds BE (39). 


Gemäß (24) und (22) erkennt man weiter, daß der Fehler das gleiche Vorzeichen hat wie q(z), 


d.h. wie die Differenzen A 4f7, also negativ ist. Um seinen Beitrag zu berechnen, muß man in 
(39) für 9(x) eine geeignete Majorante einsetzen. Wir wählen 


712),.0,32- 1029 ee Ser see (40), 


was nach Bild 1 eine sehr geeignete Majorante darstellt und erhalten damit die Abschätzung (39) 
in der Form 


1ö(2)| < ö(x) — 0,32 - 10-9 e* f@ — E)dE 

Fr . (41). 
— - 103 68 © 0,224 0,02) 
— 0,32 - 103e © 


Damit ist die gesuchte Fehlerabschätzung bereits gefunden. Für =1 erhält man ö(1)= 
0,377. 10%, Demgegenüber hat der wirkliche Fehler den Betrag ö(1)= —0,252- 10”. Die Schranke 
beträgt also tatsächlich etwa das Anderthalbfache des Fehlerbetrages. Statt auf dem angegebenen 
analytischen Wege kann man natürlich auch auf andere Weise, z.B. graphisch zum Ziel kommen. 


4. Gleichungen zweiter Ordnung bei mittelbarer Extrapolation; Beispiel 


Hier können die Überlegungen von I, Abschnitt 6 fast unverändert übernommen werden, 
nachdem auf S. 23—29erkannt wurde, daß die infinitesimale Funktion («) bei Gleichungen zweiter 
Ordnung in Bezug auf das Differenzenschema die gleiche Gestalt hat wie bei Gleichungen erster 
Ordnung. Bei unmittelbarer Extrapolation äußert sich diese Übereinstimmung in der Gleichheit 
der Formeln (16) und (1,45). Aber auch bei unmittelbarer müssen die entsprechenden Formeln 
gleich lauten, und die genannten a von I Re sich also unmittelbar auf den 
jetzigen Fall. 


Entsprechend der iterativen Verbesserung sind jetzt die aus (12) und (13) hervorgehenden 
Werte nur als erste Näherungswerte 7. und %),o zu betrachten. Aus diesen werden unter 
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Beibehaltung der gleichen Bezeichnungsweise wie in I die verbesserten Werte 


, r In Ianr 1 - 
Yaı= yiathlfno 5 Afno — 75 A’fn0 37 4° 3.) RR (49) 


und 
Yn,ı = Yn Hay tag Ar Aefa Zieh 

n— DERZ Se ED TEE 77) 

abgeleitet, welche in derselben 
Weise weiter verbessert werden bis 


As.) (43) 


sich die neu entstehenden Werte R dm 
nicht mehr ändern®). Auch hier R x "Ipeia6 
Binnenzuie Kormelnalä2Wund (Aa) ur) ame a Br 

g — — — — :15ß).025:0° e°* 


für den praktischen Gebrauch 
noch wesentlich vereinfacht wer- 
den®). Entsprechend dieser Itera- 
tion-vesteht jetzt dieFunktion g(x) 
wieder aus einzelnen Kurvenbögen, 
die sich in den Intervallendpunkten 
stetig aneinander schließen (vgl. 
Bild 2) und denen analytisch ge- 
wisse Funktionen Q,(z) entspre- 
chen (vgl. I, Bild 3). Dabei be- 
stehen, immer unter Verwendung 


derselben Schreibweise wie in I, die Bild 2 
Beziehungen 
On (&n-ı) FE 0, Q(%) eo: 0, Qu (&n+1) = Atfa+1,0 ge ‚O7 Or FE Re (44) 
und 
g(2) = Atfz (u al) x Sr —1sus0 (45) 
n+1,0 4 n 94 N 4 N 94 ’ n—1 n» 


Die Übereinstimmung von (44) und (45) mit (1,66) und (1,68) stellt den mathematischen Ausdruck 
dar für die eben erwähnte gleiche Gestalt der Funktion g(x) bei Gleichungen erster und zweiter 
Ordnung. 

Mit der gemäß (45) bestimmten Funktion g(x) kann man nun in derselben Weise weiter 
rechnen wie mit der früheren Funktion (16). Speziell können die Abschätzungen (25) und (29) 
unmittelbar übernommen werden, wenn darin nur für q (x) die Funktion (45) eingesetzt wird. 
Dabei geht man am besten wie in I, S. 12so vor, daß man eine Funktion g (x) definiert, die in den 
dreiersten Intervallen verschwindet, während sje für n> 3 vom Werte Null am Anfang des n-ten 


: til 19  - Y x { 
Intervalles in gerader Linie bis zum Wert —— |A*fa;ı,o| an seinem Ende ansteigen soll. Die un- 


241 
stetige Funktion 7(z) ersetzen wir dann gemäß (1,50) wieder durch eine mittlere Majorante 9 (2). 
Dann können wir unsere Resultate folgendermaßen zusammenfassen: 


Die in Abschnitt? ausgesprochenen Abschätzungsregeln können ohne weiteres auch auf den Fall 
der mittelbaren Extrapolation übertragen werden. Dabei sind einzig die Differenzen At fa durch die 


Werte — a Asftyıo (angenähert — ilzıne) zu ersetzen, bzw. es ist die Definition von 


T(z) und G(x) in der eben genannten Weise abzuändern. 


Die Formeln für die strenge Kontrolle können entsprechend (17) bis (20) sofort aus (9) bis (11) 
abgeleitet werden, was jedoch hier nicht explizite durchgeführt werden soll. Ferner ergibt sich wie 
in Idurch Vergleich von (16) und (45) die folgende interessante und wichtige Aussage: 

Auch bei Gleichungen zweiter Ordnung ist das Adamssche Verfahren bei mittelbarer Extra- 
polation stets rund 13 Mal genauer als bei unmittelbarer. 

Schließlich soll zur Erläuterung dieser Methoden das Beispiel von Abschnitt 3 durch mittel- 
bare Extrapolation gelöst werden. Das (1,70) entsprechende Differenzenschema nimmt dann die 
Gestalt (46) an, welche gemäß Fußnote®) auf vereinfachte Weise berechnet wurde’). Maßgebend 

6) Vgl.dazu Willers, Fr- A.: Methoden der praktischen Analysis, 2. Aufl. Berlin 1950, 8. 372—374; 
oder Collatz, L.: Numerische Behandlung von Differentialgleichungen, Berlin, 1951, S. 77 ff. 

?) In der ursprünglichen Rechnung wurden die Werte auf 8 Stellen nach dem Komma genau berechnet, 
woraus dann (46) durch Rundung der letzten Stelle hervorgeht. Auch in Wirklichkeit wurde nicht nur zweimal 
sondern viermal iteriert, was der Einfachheit halber in (46) weggelassenist. Die beiden in (46) angegebenen 
Iterationen stellen die durch einmalige und durch viermalige Iteration hervorgehenden Werte dar. Aus 
diesen Vereinfachungen erklären sich gelegentliche Abweichungen in den letzten Stellen von (46). An unseren 

rlegungen kann dies selbstverständlich nichts ändern. 


3 


ni. 
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für die Fehlerabschätzung sind die Differenzen 4,0. Sie sind in Bild 2 als Punkte eingetragen 
Damit können nun die Funktionen 9(z) und %(z) bestimmt werden (Bild 2). Auch die (für die 
Fehlerabschätzung nicht nötige) Funktion q(x) wurde interessehalber berechnet und in dieselbe 
Abbildung eingetragen. Für 9(x) wurde die Gestalt 


(z) = 0,26-10-2e. . 


en DET a N Er ER N arte > ot 


on Ym,i Ym.i In,i 10° APfn,: 107 A%fn,: 107A® fm, 107 Afn,o 
—0,1 —0,1105171 1,2156880 —2,3208589 
3208589 
0,0 0,0000000 1,0000000 —2,0000000 — 400500 
2808089 47168 
0,1  0,0904837 0,8143537 —1,7191911 — 353332 —5729 
2454757 41439 
0,2 0,1637462 0,6549846 —1,4737154 — 311893 
2142864 
0,3  0,2222460 0,5185915 —1,2594290 
41844 
0,2222454 0,5185700 —1,2593854 —311488 —5185 
2143269 36659 
0,2222454 0,5185716 —1,2593885 — 274829 
1868440 
0,4  0,2681283 0,4022081 —1,0725445 
37026 
0,2681278 0,4021886 —1,0725051 — 274462 —4517 
1868807 32509 - 
0,2681278 0,4021900 —1,0725078 —241953 
1626854 
0,5  0,3032653 0,3032785 —0,9098224 
32828 
0,3032649 0,3032616 —0,9097881 — 241634 —4047 
1627173 28781 
0,3032649 0,3032628 —0,9097905 — 212853 
1414320 
0,6 0,3297867 0,2195359 —0,7683585 
29067 
0,3292862 0,2195207 —0,7683277 —212567 —3600 
1414606 25667 
0,3292863 0,2195218 —0,7683299 — 187100 
1227506 
0,7  0,3476091 0,1489652 —0,6455794 
25721 
0,3476087 0,1489717 —0,6455520 —186846 —3197 
1227760 22524 
0,3476087 0,1489726 —0,6455539 —164322 
1063438 
0,8  0,3594622 0,0898740 —0,5392101 
22750 
0,3594818 0,0898620 —0,5391858 —164096 —2836 
1063664 19914 
0,3594619 0,0898628 —0,5391875 — 144182 
919482 
0,9  0,3659114 0,0406640 —0,4472393 
20114 
0,3659111 0,0406533 —0,4472177 — 143982 — 2521 
919682 17593 
0,3659111 0,0406541 —0,4472192 — 126389 
793293 
1,0  0,3678778 0,0000061 —0,3678900 
0,3678775 —0,0000034 —0,3678707 


0,36787756 —-0,00000273 


Genaue Werte: y (1, 0) = 0,36787944 und y’ (1, 0) = 0,00000000 


/ | 
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angenommen, welche sehr angenähert in halber Höhe der Spitzen von 9(x) verläuft. Damit kann 
nun genau die entsprechende Rechnung durchgeführt werden wie im Beispiel von Abschnitt 3. 
Der einzige Unterschied gegenüber diesem letzteren besteht darin, daß an Stelle von (40) jetzt die 
dazu proportionale Funktion (47) tritt. Infolgedessen verändert sich auch das Resultat (41) 
in demselben Verhältnis, und wir erhalten für den Fehler die Abschätzung 


ol <a) =0,26- 10-10 ( —0,22+0,02) EN HL EN 


Für ©=1 hat man z.B. ö(1)=0,306-10=. Dieser Schranke steht der wirkliche Fehler 
ö(1)=0,188.10° gegenüber. Das Verhältnis ist also wieder ungefähr 3/2. Das Vorzeichen des 
Fehlers bestimmt sich auch hier aus (21), wobei jetzt g(x), wie immer bei mittelbarer Extra- 
polation, das entgegengesetzte Vorzeichen hat wie die Differenzen A? fo. Schließlich findet 
man auch die oben ausgesprochene Regel bestätigt, daß das Resultat bei mittelbarer Extra- 
polation stets rund 13 Mal genauer ist als bei unmittelbarer. 


5. Gleichungen höherer Ordnung 


Hier treten keine neuen Gesichtspunkte mehr auf, und es genügt wohl, die sich ergebenden 
Resultate am Beispiel der Gleichung dritten Grades kurz zu skizzieren. Bei einer Gleichung 
dritter Ordnung 


EEE EN. ie (49) 
mit den Anfangsbedingungen 
y(m) = Yo» y(m) = Yo» Yr)=Ww 0.0.00 (60) 
lautet die zu einer Näherungslösung 4*(x) gehörende Fehlergleichung entsprechend (4) 
6 — nu Er (2 yr, N y*") — Ö nr (z, 4", ur, ee) 
of r n ti " [ZA ; 3 ; (51) 
art, VER UN) a EYE | 
J 
oder in abgekürzter Schreibweise analog zu (5) 
BR 0 Nr. SER RENT 52) 
wozu noch die Anfangsbedingungen 
d(%) = 0 (H) > 0 6".(%) = 0 TE EEE « (53) 


hinzukommen. Bei der Integration nach Adams ist wieder zu berücksichtigen, daß sich y*(x) 
aus den einzelnen Teilen Y7 (x) zusammensetzt, die für n > 3 in Verbindung mit dem zugehörigen 
Differenzenschema fortlaufend für die verschiedenen Intervalle gebildet werden ®). Dabei wird 
Y% (x) in genau gleicher Weise aus dem Differenzenschema abgeleitet wie (9) aus (8). Die rechte 
Seite von (9) stellt demnach bei Gleichungen dritter Ordnung die dritte Ableitung des Extra- 
polationspolynoms Yy(x) dar, aus der sich dieses letztere und seine beiden ersten Ableitungen 
durch einfache Integrationen ergeben, genau wie (10) und (11) aus (9). Die rechten Seiten von 
(10) und (11) stellen jetzt also die Ableitungen Y,”, (x) und Y7’(z) dar, während zur Berechnung 
von Y,(x) eine nochmalige Integration nötig ist. 

Mit Y, (x) und seinen Ableitungen kann nun weiterhin vor allem die Funktion g(x) gebildet 
werden. Die Berechnung erfolgt in ganz analoger Weise wie dies auf S.28—29 für Gleichungen zwei- 
ter Ordnung durchgeführt wurde. Speziell ergibtsich auch bei Gleichungen dritter Ordnung für (x) 
wieder die Gestalt (16) bzw. bei mittelbarer Extrapolation die Form (45). Dieselben Schlüsse 
lassen sich auch für Gleichungen höherer Ordnung ziehen, und wir können also das folgende all- 
gemeine Prinzip für unsere Fehlerabschätzung aussprechen: 


Eine beliebige gewöhnliche Differentialgleichung sei nach Adam s integriert. Dann bestimmt 
sich der dabei begangene Fehler aus der Fehlergleichung, welche entsprechend (4), (5), (51) und (52) ın 
elementarer Weise aus der gegebenen Gleichung abgeleitet werden kann. Die in der Fehlergleichung 
als Störungsglied auftretende infinitesimale Funktion q(x) hat bei allen Gleichungen in erster Nähe- 
rung die Form (16) bzw. bei mittelbarer Extrapolation die Form (45) und kann im übrigen mit großer 
Leichtigkeit genau berechnet oder abgeschätzt werden. Die übrigen in der Fehlergleichung auftretenden 
Funktionen ergeben sich sofort aus der Lösung der Differentialgleichung. Man erhält dann eine Ab- 


 schätzung für den Fehler, indem man die Fehlergleichung integriert oder ihre Lösung abschätzt, was je 


4 
r 
” 


8) Vgl. die unter *) genannten Werke. 
3* 


\V = R ” : C Z. angew. Math. Mech. 
36 Matthieu, Fehlerabschätzung beim Extrapolationsverfahren von Adams „4 33 Nr. 1/2 Jan./Febr. 1953 


nach ihrer Gestalt in geschlossener Form oder approximativ erfolgen kann. Dabei stehen sehr ver- 
schiedene Möglichkeiten zur Verfügung, deren wichtigste vorstehend und in I besprochen wurden. Auch 
die bei diesen Rechnungen gemachten Vernachlässigungen können stets kontrolliert werden, so daß 
diese Methoden als vollständig streng betrachtet werden können. 

Weiterhin folgt ganz allgemein: 

Bei sämtlichen Differentialgleichungen führt die Integration nach Adam s bei mittelbarer 
Extrapolation zu rund 13 Mal genaueren Resultaten als bei unmittelbarer. 


Auch die beiden Methoden zur Lösung bzw. Abschätzung der Fehlergleichung, wie sie für 
Gleichungen zweiter Ordnung angegeben wurden, können unmittelbar übernommen werden. 
Bei der ersten (S.29—30) muß wiederum zuerst das allgemeine Integral der homogenen Gleichung 
bestimmt werden. Sind z.B. ö,(x), ö,(x) und ö,(x) drei linear unabhängige Integrale der zu (52) 
gehörenden homogenen Gleichung, so stellt sich die Lösung der Fehlergleichung (52) unter Be- 
rücksichtigung der Anfangsbedingungen (53) entsprechend (24) in der Form dar 


z 


1 — Ss). 
ld) ae D.(x, ae Re 
wobei D(x, &) jetzt analog zu (23) eine Abkürzung ist für 
dl) del) 5 (e) 
Dd=- 0) 0 5%) 
CEO NG 
und wobei Ü wiederum eine Konstante ist, die von der Wahl von 6, (z), 6,(x) und ö, (x) abhängt °). 
Die Identität, gemäß welcher dies je der Fall ist, lautet jetzt entsprechend (22) 


dl)  4(E) B RE Ä oa 
(ld) (8) ()=(e” BE re 
8 8 


Sie kann wiederum leicht zur Berechnung von C dienen. 
Aus (54) folgt dann die (25) entsprechende Abschätzung 


. (55) 


BL 
lö(z)| <ö(e) = co ro ID, EildEr ae rel): 


Für sie gelten die gleichen Bemerkungen wie für (25), wobei wiederum das Vorzeichen von D(z, &) 
eine ausschlaggebende Rolle spielt. Ganz entsprechende Verhältnisse bestehen für Gleichungen 
höherer Ordnung. All dies soll jedoch hier nicht mehr näher ausgeführt werden, da dieser Lösungs- 
weg bereits von der dritten Ordnung an praktisch recht kompliziert wird. 
Sehr einfach ist dagegen die Übertragung der zweiten Methode der Fehlerabschätzung, wie 
„sie S.30—31 entwickelt wurde. Bei Gleichungen dritter Ordnung erhält man an Stelle von (26) 
das Gleichungssystem 


u = söhr + röhr + Po + Au: Sr (2) = du (2) = du (m) = 0 fe. (88). 


Macht man hier zunächst die Annahme, daß p(x) und r(x) im Integrationsintervall positiv sind 
(die übrigen Fälle werden im nächsten Abschnitt behandelt), so können wir aus dem auf S. 30 bis 
31 ausgesprochenem Theorem, das sich in ganz analoger Weise auch für Gleichungen höherer 
Ordnung formulieren läßt, sofort schließen, daß die Lösung öy,(&) wieder eine Majorante ist für den 
Fehler ö(z). Genau entsprechend liegen die Verhältnisse bei Gleichungen n. Ordnung. 


6. Die zweite Abschätzungsmethode im allgemeinen Fall 
Die zweite der beiden Methoden, die in einer direkten Integration der Fehlergleichung 
besteht und die im Prinzip die einfachere und ansprechendere ist, ist nach den bisherigen Aus- 
führungen bei Gleichungen zweiter Ordnung auf den Fall p(z2)Z0 beschränkt, während bei 
Gleichungen höherer Ordnung noch mehr solche Ungleichungen hinzutreten. Wenn diese Voraus- 
setzungen nicht mehr erfüllt sind, ergeben sich wesentlich neue Verhältnisse. Das auf S. 30—31 
ausgesprochene Theorem gilt dann nicht mehr. Tatsächlich erkennt man aus (24) und (54) sofort, 


?) Vgl.z.B. Forsyth, A. R.: Lehrbuch der Differentialgleichungen, 2. Aufl. Braunschweig 1912, $ 75, 
insbes. Beispiel 1. 


a I 
- % « “ R - 
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daß im allgemeinen Fall (nämlich dann, wenn 'D(z, £) das Vorzeichen wechselt) eine Majorisierung 
von g(x)in der Fehlergleichung nicht mehr eine solche der Lösung zur Folge zu haben braucht, 
wie das bei Bestehen der genannten Ungleichungen gefolgert werden konnte. Es sind demnach 
in diesem Falle neue Gesichtspunkte heranzuziehen. Als wichtigste Methoden der Fehlerab- 
schätzung ergeben sich dabei die folgenden: 

a) Oft ist es möglich, die Fehlergleichung (5) durch eine einfache Umformung in eine lineare 
Gleichung mit p(z) = 0 überzuführen. Z. B. kann das Gleichungssystem 


"= —ö-tgle), LE 1 VE ME 
im IntervallO< x< 1 durch die Substitution 
= usin(e+%) EEE A 9 00) 
übergeführt werden in 
u" = —2 u cotg (.+ 2) Fr ] N BUN OFFEN: (61): 
# sin (e —- 2 i 
Für große Intervalle ist (60) jedoch nicht mehr brauchbar, da dann die Koeffizienten unstetig 


werden. 
Auch das Beispiel von Abschnitt 3 kann auf diese Weise berechnet werden. Die Fehler- 
7 gleichung (36) geht durch die Substitution | 


DENE ae ra er(B2 
über in 
er Yayeu (Val etrerlars (03). 
Ersetzt man hier q wieder durch die Majorante (40), so kommt man genau auf die frühere Lösung 
(41) zurück. Im Gegensatz zu (61) ist (63) für beliebige Intervalle brauchbar. 


b) Eine andere Möglichkeit besteht darin, daß man die Fehlergleichung (5) zunächst durch 
ö’= ein ein lineares Gleichungssystem erster Ordnung überführt und dann in diesem Linear- 
kombinationen von ö und e als neue Veränderliche einführt. Z.B. kann man in dem aus (59) 
hervorgehenden linearen Gleichungssystem 


O4 8(0)=0 | 


NE 
= —ö+gla)  2()=0 | 1) 
die neuen Veränderlichen v und w einführen durch 
5= p(a)v + yla)w, e=rl@aoto@w 2.2... (65) 
mit zunächst unbekannten Funktionen $(z), y(x), x(2), (x). Es folgt dann für v und w das 
Gleichungssystem 
(we — yX)vV = (we - op +Yy+xXy)o+ (oe — yo+y?+o'y)w—qgy 
) (we  y)wW = (— —- XP — 2 +tpr)o + pPy — 90’ —go+y'x)w-tgq 2 (66). 


v(0) = 0 (0) = 0 


Hier ist es offenbar auf mannigfache Weise möglich, die obigen Funktionen so zu wählen, daß der 
Koeffizient von win der oberen Gleichung, derjenige von v in der unteren sowie der bei den Ab- 
leitungen stehende Koeffizient alle > 0 werden. Z. B. kann man 


P=cosT, y=sinz, = -sinz, DIS OE Dr art OL) 
setzen, wodurch (66) die Gestalt 
v”=—gsns v()=0, BO =O N. 2% ...2(68) 


annimmt. Ersetzt man hier die rechten Seiten — qsin z und gcos x durch Majoranten bzw. 
Minoranten, so ergeben sich durch Integration der entstehenden Gleichungen solche für v und w 
sowie schließlich durch Vermittlung von (65) Majoranten und Minoranten für ö und e, wie das der 
Leser selbst im einzelnen ausführen möge. Damit ist die Fehlerabschätzung geleistet. 

Die vorstehend am Beispiel von (59) ausgeführte Methode kann in gleicher Weise auch auf 
den Fall der allgemeinen Gleichung zweiter und höherer Ordnung ausgedehnt werden. Dabei 
läßt sich vor allem an Hand von (24) und (54) unschwer einsehen, daß man stets auf mannig- 
fache Weise solche linearen Kombinationen von ö und seinen Ableitungen als neue Veränderliche 
einführen kann, daß bei dem entstehenden linearen Gleichungssystem die notwendigen Unglei- 
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chungen erfüllt sind. Jedoch soll dieser zunächst außerordentlich naheliegend und vielverspre- 
chend erscheinende Weg hier trotzdem nicht weiter verfolgt werden, da es dem Verfasser bisher 
nicht gelungen ist, ihn trotz seiner theoretischen Einfachheit zu einer wirklich brauchbaren 
Methode auszubauen. 

c) Ein praktisch sehr geeigneter Weg besteht dagegen darin, daß man die Fehlergleichung 
ein oder mehrere Male ableitet, die so entstehenden Gleichungen mit geeigneten (konstanten oder 
veränderlichen) Faktoren multipliziert und sie schließlich addiert. Man erhält dann eine neue 
Gleichung für den Fehler. Diese ist zwar von höherer Ordnung, jedoch kann man versuchen, die 
genannten Faktoren so zu wählen, daß die in Frage stehenden Ungleichungen erfüllt sind. Dabei 
sind natürlich auch die Anfangsbedingungen auf die Gleichung höherer Ordnung zu übertragen. 

Eine kleine Schwierigkeit ergibt sich dabei noch dadurch, daß g(x) als unstetige Funktion 
an den Sprungstellen nicht differenzierbar ist. Diese kann aber leicht dadurch überwunden 
werden, daß man g(x) durch eine stetig verlaufende Funktion K(z) ersetzt, welche g(x) unter 
Vermeidung der Zacken im Mittel möglichst gut approximiert, genau wie g(x) durch 9(x) approxi- 
miert wird, wenn in (1,50) das Gleichheitszeichen steht. X(x) kann dabei als sehr rohe Annähe- 
rung gewählt werden, wie das nächstens noch deutlicher ersichtlich wird. Eine solche kann stets 
leicht gefunden werden. 

Ist nun wieder z.B. die Gl. (59) vorgelegt, so ergeben sich durch zweimalige Ableitung die 
Differentialgleichung e 
0" + KEN 9 KHK Er (69) 


sowie die Anfangswerte 
ö"(0) = K(0) 2ö (O). == RE (0), Re 

Damit kann also (59) ersetzt werden durch das Gleichungssystem 
"=6—K-+K", 5(0,=5.0)=,0, OO) FERKO): 6” (0)= K’(0) (71), 


bei welchem die fraglichen Vorzeichenbedingungen erfüllt sind. Ersetzt man nun in (71) die Funk- 
tion — K(x) + K"(«) und die Anfangswerte 0,0, K(0), X’(0) durch irgendwelche Majoranten 
(wozu sehr rohe Abschätzungen genügen),so ist die Lösung öx(x) des so entstehenden Gleichungs- 
systems eine Majorante für den Fehler. In gleicher Weise kann auch eine Minorante ö,(x) her- 
geleitet werden. 

In genau gleicher Weise folgt für den Fehler im Beispiel von Abschnitt 3 das Gleichungs- 
system 

6" =36+25—2K+K, 5.(0,2) = 8’(0,2)= 0, 6" (0,2) = K (0,2) (72). 
Ersetzt man hier gemäß Bild 1 die Funktion 
—2K-+ K’ durch die Majorante 0,8-103 e® 
bzw. durch die Minorante 0,5-103e? und 
K (0,2) durch die Werte — 0,15-103 bzw. 
— 0,2-103, so ergeben sich als Lösungen der 
entstehenden Gleichungssysteme die in Bild 3 
dargestellten Funktionen öx(x2) und ö.(x), 
welche leicht approximativ (und im obigen 
Fall auch leicht exakt) bestimmt werden kön- 
nen. Auf Grund von diesen ergibt sich also die 
Abschätzung 
(JS) SHlT . » --. (2). 

In ähnlicher Weise kann man auch in 
anderen Fällen vorgehen. Die weitläufige 
Frage, ob es stets möglich ist, für einen Fehler 

Bild 3 nach der obigen Methode eine Gleichung mit 

den fraglichen Vorzeichenbedingungen abzulei- 

ten, soll hier nicht berührt werden. Ein Nachteil dieser Methode besteht natürlich darin, daß 
die Ordnung der zu lösenden Gleichung erhöht wird. 

d) Die letzte und wichtigste Methode beruht schließlich darauf, daß die Integration der Feh- 
lergleichung bei Ersetzung von g(z)durch eine geeignete Mittelfunktion auch im jetzigen Fall weit- 
gehend eine Näherung für den Verlauf des Fehlers liefert, welche überdies in der überwiegenden 
Mehrzahl der Fälle bereits zur Fehlerabschätzung brauchbar ist. Es erscheint daher naheliegend, 
die Fehlergleichung auch jetzt direkt zu integrieren und dabei solche Korrekturen anzubringen, 
die der Unsicherheit des jetzigen Falles Rechnung tragen. Am besten geht man dabei folgender- 
maßen vor, wobei wir uns wieder auf Gleichungen zweiter PrunUng) beschränken: 
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Zunächst sei nochmals daran erinnert, daß wir im Falle p(x)> 0 auf Grund des am Schluß 
von Abschnitt 2 ausgesprochenen Satzes zum Ziel kommen konnten unter Verwendung der auf 


S.29/30 eingeführten Funktionen I (%), Im (X), dr (8) und na). Im jetzigen Falle ist die Be- 
dingung p(x)> 0 jedoch gerade nicht erfüllt. Um hier weiter zu kommen, teilen wir das Inte- 
dnsttercall in einige Teilintervalle der gleichen Länge /, deren Mittelpunkte mit m;, 
ET M„ bezeichnet seien. Ferner bezeichne s; bei positivem q(m;) den Wert Tu (m;), bei 
negativem dagegen den Wert 9,(m;). Integriert man dann die zu (5) bezw. zu (26) gehörende 
homogene Gleichung 


a ee (74) 
unter den n verschiedenen Anfangsbedingungen 
6, (m.) ==(); ö,(m;) SM AoH ea) nen (79) 


so ergeben sich » Lösungen Ö,,;(z), die wir uns jeweils links von ihrem Anfangspunkt durch die 
identisch verschwindende Funktion ersetzt denken wollen. Die Summe 


nr 
OT N 2 kat u a ah ra (70) 
i=1 
der so definierten » Lösungen stellt dann eine gewisse Näherung für den Fehler ö(x) dar, die 
jedoch in keiner Weise eine Majorante zu sein braucht. 

Indessen erhält man sicher eine solche, wenn man diejenigen Teile der Lösungen Ös, ;(x), die 
negativ sind, einfach wegläßt, die betreffende Lösung dort durch Null ersetzt und dann wieder die 
Summe der so abgeänderten Lösungen bildet, die jetzt mit Ayr(x) bezeichnet sei. Ebenso erhält 
man eine Minorante Anz), wenn man die positiven Teile wegläßt und die Lösungen dort durch 
Null ersetzt. Es läuft dies darauf hinaus, daß man q (x) in jedem Intervall entweder durch eine 
Majorante oder eine Minorante ersetzt (wobei eine von ihnen stets identisch verschwindet), und 
zwar So, wie es für jeden Wert von x am ungünstigsten ist. Daß ein solches Vorgehen zu einer 
Majorante bzw. zu einer Minorante für den Fehler führen muß, ergibt sich unmittelbar aus der 
Darstellung (24). Die Überlegungen sind um so strenger, je gr ößer die Anzahl n der Intervalle 
gewählt wurde. Damit hat man dann also eine Abschätzung i in der Form 


BRASS PTENE TE IE DU NE TORE MINE IR DREI 


Dieses zunächst umständlich erscheinende Verfahren ist in der Praxis sehr einfach durch- 
zuführen. Es genügen dabei natürlich wieder sehr rohe Integrationsmethoden. Ferner braucht 


man durchaus nicht alle Lösungen ö,,;(x) zu berechnen. Vielmehr genügt die Berechnung einiger 
weniger von ihnen, aus welchen dann die übrigen durch Interpolation ; approximativ festgestellt 
werden können. Schließlich kann man auch die sämtlichen Lösungen ö,,;(x) auf Grund von (75) 
linear kombinieren aus zwei Grundlösungen 6, (2) und ö, (x) der zu (5) und (26) gehörenden homo- 
genen Gleichung. Besonders einfach ist dieses Verfahren auch bei maschineller Integration durch- 
zuführen. 

Als Beispiel soll wieder die Aufgabe von 
Abschnitt 3 durchgerechnet werden. Wir teilen 
dazu das Intervall 02 Sxzs 1,0in vier gleiche 
Teile der Länge != 0,2. Es bestehen dann 


2, (x)s0 
FF rd (8) 


die Beziehungen “2 
m=0,3, m= 0,5, | (78) 
Mg = 0,7, My = 0,9 J 2 


Ferner können auf Grund von Bild 1 die Werte 
kı = — 0,24-10%, = —0,19-.10, 
(79) ro 
k; = — 0,16-107, = — 0,13:107 
angenommen werden. Weiterhin nimmt die 
Gleichung (74) gemäß (36) die Gestalt an Bus 4 
nt... er (80). 
Es handelt sich nun darum, (80) unter den vier Anfangsbedingungen (75) zu integrieren, wenn in 
(75) die Werte (78) und (79) eingesetzt werden. Das Resultat ist in Bild 4 dargestellt. Von diesen 
vier Kurven wurden Ö,, ı(2) und Ö,, 3() direkt durch approximative Lösung s der obigen Gleichun- 
gen berechnet, während d,,2(8) und Ö,,1(%) nach Augenmaß durch Inter- bzw. Extrapolation 
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bestimmt wurden. Da hier die Lösungen 6,,,(8) überall negativ sind, müssen sie zur Bildung von 
Ay() alle durch die identisch verschwindende Funktion ersetzt werden. Es ist also Ay(z) = 0, 
während sich andererseits A„(x) direkt als Summe der vier Lösungen darstellt. 


Schließlich sei noch bemerkt, daß man in absolut analoger Weise auch bei Gleichungen 
höherer Ordnungen vorgehen kann, wie das hier jedoch nicht im einzelnen durchgeführt 
werden soll. 


7. Schlußbemerkungen und Zusammenfassung 


In der vorliegenden Arbeit wird die in I für Gleichungen erster Ordnung entwickelte 
Methode der Fehlerabschätzung auf Gleichungen höherer Ordnung ausgedehnt. Als Hauptunter- 
schied ergibt sich dabei, daß auch die Fehlergleichung eine höhere Ordnung annimmt. Für die 
Fehlerabschätzung werden dann die fünf Funktionen öy(2), 6m (x), Ay (&), An (x) und ö(x) 
eingeführt, die in manchen Fällen übereinstimmen, und von denen sich die beiden ersten gemäß 
(26) als Lösungen der Fehlergleichung ergeben, wenn man in dieser q (z) durch Qu (x) bzw. durch 
Am(x) ersetzt. Die Bildung von Ay(z) und (x) ist weiterhin auf S. 39 beschrieben, während 
diejenige von (x) gemäß (25) erfolgt. 

Für die Fehlerabschätzung wichtig sind in erster Linie die vier a der genannten Funk- 


tionen. Sie können stets in ganz elementarer Weise bestimmt werden. Dabei sind verschiedene 


Fälle zu unterscheiden. Wenn erstens gewisse Vorzeichenbedingungen erfüllt sind (bei Gleichungen 
zweiter Ordnung muß z.B. p(x)> 0 sein), können die Funktionen öu,(x) und d,„(x) direkt 
gemäß (29) zur Fehlerabschätzung verwendet werden, wobei meist die Kenntnis einer von ihnen 
genügt. Die Bestimmung kann dabei entweder nach der ersten der beiden genannten Methoden 
gemäß (27) (bzw. einer entsprechenden Formel bei Gleichungen höherer Ordnung) erfolgen oder 
dann nach der S. 30—31 und S. 36—38 beschriebenen zweiten Methode durch direkte, im all- 
gemeinen approximative Integration der Fehlergleichung. Den ersten Weg wird man vor allem 
dann einschlagen, wenn die zur Fehlergleichung gehörende homogene Gleichung in geschlossener 
Form oder sonst einfach integrierbar ist. Ebenso ist diese Methode z.B. bequem in dem häufig 
vorkommenden Fall, daß die zu lösende Differentialgleichung linear und homogen ist. Denn in 
diesem Fall stimmt die zur Fehlergleichung gehörende homogene Gleichung mit der ursprünglich 
gegebenen Gleichung überein. Die nach Adams gefundene Lösung dieser letzteren ist dann also 
gleichzeitig auch eine Lösung der zur Fehlergleichung gehörenden homogenen Gleichung. 

Wenn dagegen die genannten Vorzeichenbedingungen nicht erfüllt sind, kann man zu- 
nächst versuchen, dies durch Umformung der Fehlergleichung zu erreichen, wozu vor allem die in 
Abschnitt 6 unter a) und c) genannten Methoden zur Verfügung stehen. Andernfalls liefert das 
unter d) beschriebene Verfahren stets die Funktionen A,,(2) und A,„(z), aus welchen sich gemäß 
(77) eine Abschätzung ergibt. Auch hier braucht man meistens nur eine der beiden Funktionen zu 


berechnen. 
Schließlich kann oft, und zwar wieder vor allem dann, wenn die homogene Gleichung in 


geschlossener Form lösbar ist, auch die Verwendung der Majorante ö(z) für den Fehlerbetrag von 


Nutzen sein, welche übrigens in sehr vielen Fällen dem Betrag nach mit u (&) oder mit m (2) 
übereinstimmt. Man vergleiche dazu die in den Abschnitten 3 und 4 durchgerechneten Beispiele. 


Diese Majorante ist außerdem von prinzipiellem Wert; vor allem hat sie den Vorteil, daß sie ohne 


Einschränkung in allen Fällen gültigist. Eine weitere Möglichkeit der Fehlerabschätzung, die mit 


ö(&) in engem Zusammenhang steht und die unter Umständen praktisch sehr wertvoll sein kann, 
ist auf S. 38 und 39 ausgeführt. 


Die Majorante ö(x) ist weiterhin die direkte Verallgemeinerung der in I gleich bezeichneten 
Majorante. Jedoch liegen die Verhältnisse bei Gleichungen erster Ordnung insofern wesentlich 
einfacher, als D (x, &) dort automatisch stets das gleiche Vorzeichen hat, so daß es am einfachsten 


ist, mit ö(x) zu operieren. Ergänzend sei jedoch an dieser Stelle noch zu I hinzugefügt, daß es 
selbstverständlich auch bei Gleichungen erster Ordnung möglich ist, die Funktionen öy.(x) und 
ö,, (x) einzuführen. Gemäß (1,22) nimmt z.B. öy,(&) die Gestalt an 


— fpün | pda | 
öu (2) = e gi Gu(x) e®.. da. u 0 er. a 


% 


Die Verwendung von du (x) und On (x) ist sogar bei Gleichungen erster Ordnung besonders sinn- 
gemäß, weil die bei Gleichungen höherer Ordnung so beträchtliche Schwierigkeiten bereitenden 


Ö 
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Vorzeichenbedingungen dort automatisch erfüllt sind, so daß (29) immer gilt. Im übrigen ist der 


Unterschied zwischen der Verwendung von ö(x) einerseits und u (x) und ö,,(x) anderseits bei 
Gleichungen erster Ordnung völlig unwesentlich. Die Verwendung der letzten Funktionen ist vor 
allem dann gegeben, wenn g (x) das Vorzeichen öfters wechselt. 

Bei all den angegebenen Fehlerabschätzungsmethoden ist schließlich noch im Auge zu 
behalten, daß es sich in erster Linie um prinzipielle Möglichkeiten handelt, nach denen eine Ab- 
schätzung vorgenommen werden kann. Die angegebenen Wege können, je nach der gestellten 
Aufgabe, noch auf die verschiedenartigste Weise modifiziert, verbessert und vereinfacht werden. 
“ Ähnlich wie beim Verfahren von Adams selbst ist es ein besonderer Vorteil dieser Methoden, daß 
sie sich zwanglos allen speziellen Forderungen anpassen lassen. Ferner kann man auch versuchen, 
bereits die gegebene Gleichung auf eine zur numerischen Integration und zur Fehlerabschätzung 
günstige Gestalt umzuformen. Beispielsweise wird es in dem häufig vorkommenden Fall einer 
linearen homogenen Gleichung zweiter Ordnung in der Regel vorteilhaft sein, diese in eine 
Riccatische Gleichung erster Ordnung überzuführen. 


Eingegangen am 20. Dezember 1951. 


Taylor=-Entwicklungen und finite Ausdrücke für zu und 4au 
Von J. Albrecht in Haste bei Hannover 


Die zur Herleitung finiter Ausdrücke für Au und AAu erforderliche Taylor - Entwicklung eines all- 
gemeinen Ansatzes wird durch geeignete Zusammenfassung partieller Ableitungen in eine übersichtliche Form 
gebracht, aus der leicht zu ersehen ist, welche Möglichkeiten hinsichtlich des Taylor - Abgleiches bestehen. 
Es werden regelmäßige Gitter aus Quadraten, Dreiecken, Sechsecken und Würfeln zugrunde gelegt. 


The Taylor series expansion of a general statement that leads to finite-difference approximations of 
Au and AAu is brought, by a suitable composition of the partial derivatives, into a clear form. Thhus it is easy 
to survey the possibilities of the procedure. T’he author discusses regular grids of squares, triangles, hewagons, 
and cubes. 


Le developpement de Taylor d’un arrangement general, lequel est necessaire pour la derivation d’ew- 
pressions finies pour Au et AAu, est reduit par une composition convenable de derivations partielles a une 
forme claire, d’oü il est facile de voir, quelles possibilites existent & l’egard de l’ajustement de Taylor. 
Des grilles regulieres de carr£&s, triangles, hexagones et cubes sont supposees. 


TatnopoBo pasıokHue HEKOTOPOTO BbIpasKeHHA B OÖIIEeM BuNe, TpeöyeMmoe Alf BBIBONA 
KOHEYHBIX BbIpaskeHnnä ua Au am AAu, IPHBONUTCA, IIYTeM HAaNIesKaıeroO COYeTAHHA YACTHEIX 
IpPOHSBOAHBIX, K HATJIAAHOMY BHAY, H3 KOTOPOTO JIeETKO VCMATPHBAIWTCH BO3MOKHOCTH TOHIOPOBA 
COTJIACOBAHHAN. B OCHOBY IHOJIOSKEHBI PETYAIAPHbIeE PeImeTKH MW3 KBAAPaATOB, TPEYTOJIBHHKOB, 
INeCTHYTOAIBHHKOB HU KYÖOB. 

Unter Verwendung regelmäßiger Netze aus Quadraten, Dreiecken und Sechsecken bzw. 
aus Würfeln werden im folgenden für Funktionen u(z, y) bzw. u(x, 9,2) Taylor - Entwick- 
lungen von A’u(n =0,1,2,...; A=Laplacescher Operator) aufgestellt und einfache finite 
Ausdrücke für Au und AAu angegeben ). Dabei soll vor allem gezeigt werden, wie sich durch 
Zusammenfassung partieller Ableitungen zu geeigneten Operatoren eine übersichtliche Darstel- 
lung der Taylor -Entwicklung erzielen läßt, die zwei Vorteile bietet: 


1. die Zahl der beim Taylor - Abgleich zu lösenden Gleichungen wird verringert. 
2. Es ist leicht feststellbar, 
a) von welcher Ordnung das Restglied eines finiten. Ausdruckes bei gegebenen 
Gitterpunkten höchstens werden kann, und 
b) wieviele Punkte zur Aufstellung eines finiten Ausdruckes mit einem Restglied von ge- 
wünschter Ordnung mindestens erforderlich sind. 


Besonders günstig sind finite Ausdrücke für Au, die nur die Entwicklungsstelle und die 
Punkte des ersten umliegenden Ringes verwenden. Dann läßt sich gewöhnlich vermeiden, daß 
beim Aufschreiben finiter Gleichungen in randnahen Gitterpunkten auch Werte an Stellen auf- 
treten, die außerhalb des gegebenen Grundgebietes liegen. 


A. Quadratische Netze 


%, bezeichne den Funktionswert an der Stelle des Netzes, an der ein finiter Ausdruck auf- 
gestellt werden soll; u, sei die Summe der Funktionswerte in den unmittelbar benachbarten 


1) Eine allgemeine und ausführliche Darstellung des Differenzenverfahrens findet man in: L.Collatz: 
1. Numerische Behandlung von Differentialgleichungen, Berlin 1951; 2. Eigenwertaufgaben mit technischen 
Anwendungen, Leipzig 1949. 3.Einige neuere Forschungen über numerische Behandlung von Differentialglei- 
chungen, Abschn. I 1d. Z.angew. Math. Mech. Bd. 31 (1951), S. 230—236; 4. Das Mehrstellenverfahren bei 
Plattenaufgaben, Z. angew. Math. Mech. Bd. 30 (1950), S. 385—388. [Im Text und in den Tabellen eingefügte 
Hinweisziffern beziehen sich auf diese Literaturangaben.] 
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Gitterpunkten, I%,, Nu;, ... die Summe der Werte in weiteren, jeweils im gleichen Abstande 
von der Entwicklungsstelle gelegenen Punkten (s. Bild A). Entsprechendes gelte für die Diffe- 
rentialausdrücke Au,, YAw, YAu,, ...; Adu,, YAAw, NAAU,,... usw. hist die Seiten- 
länge der Maschen des Netzes. 

Nun wird (mit a,b, o,...; A, B, C,.... als noch zu 
bestimmenden Koeffizienten ?)) der Ausdruck 


law +b3wte Bu tdiwmt-::} 
+51 {Au,+ B 2m + C YAm,+D Au, ---} 
4m) 


gebildet und an der erwähnten Stelle nach dem Ta y- 
lorschen Satz entwickelt (s. Schema A), wobei durch % 


Zusammenfassung partieller Ableitungen Operatoren 
Ar De (p,gq=0,1,2,...)eingeführt werden mit 


0° 0° o* 
Zur näheren Erläuterung des unter 2a und 2b Bild a 


Gesagten stellen wir uns die Aufgabe, für die Dif- 

ferentialgleichung Au = r(&, y) finite Gleichungen zu ermitteln, wollen dazu aber Funktions- 
werte u sowie Ausdrücke Au{=r(z, y)}, Adu {= Ar(z, y)},... nur in den Punkten heran- 
ziehen, die in Bild A 1 angegeben sind. — Bei der Durchführung des Ta ylor- Abgleiches (d.h. 
Nullsetzen der Faktoren der in der Taylor- Entwicklung auftretenden Operatoren bis zu 
Gliedern möglichst hoher Ordnung hin) sind dann u.a. folgende Gleichungen zu lösen: 


Faktor von x U : a +4(b+ c+ J=0 Erfüllt für: 
24 a — 300 

Faktor von 7i Du, : 4(b+ 16c— 4e)=0 b= —64 
18 Bl 

Faktor von 5; D’u, : 4(b+ 256c +16e) = 0 e=.—12 
h!? 

Faktor von 151 D’u,: 4(b-+ 4096 ce — 64e) — 0 


Die eingeklammerte Gleichung wird von den Werten, die durch 

die ersten drei Gleichungen bis auf einen gemeinsamen kon- 

stanten Faktor eindeutig festgelegt sind, nicht mehr befrie- 
mM : dist, im günstigsten Falle kann also ein Restglied 12. Ordnung 
e ; erreicht werden. (Die Koeffizienten der 


übrigen Glieder der Taylor-Ent- 8. 
% © wicklung bis zur 8. oder 10. Ordnung 
oO e einschließlich lassen sich durch geeignete en. 
Wahl der Größen A, B, C, E; a, ß,y, 
&;... zum Verschwinden bringen.) — ONE 
ARE Ebenso zeigt man: Sollen nur die Ent- BildA2 


wicklungsstelle und die Punkte des um- 
liegenden Ringes herangezogen werden (Bild A2), so lassen sich für Au= r(z, y) finite Ausdrücke 
mit einem Restglied höchstens 8. Ordnung bilden, und zwar wegen: 


0 .. .. 
Faktor von w u :a+4b+ o)=0 Erfüllt für: 


hi b — —4 
Faktor von ar Dun: 4(b— 4e) =0 | ER 
78 
Faktor von 25 DW: 4(b +16) | 


2) Verwechselungen mit a, b, c,... als Punktbezeichnungen sind wohl nicht zu befürchten. 


ige ia ar, a 
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Abschließend werden in Tabelle A einfache Formeln zusammengestellt, die man aus der 
obigen allgemeinen Taylor - Entwicklung herleiten kann. Formel (3) ersetzt die Differential- 
gleichung Au +Au = 0 durch finite Gleichungen der Form 


Su, +4 3w=v:u mit Ki 12 u »[2] 
Tabelle A 
Nr. Formel Restglied | 
ik Zu — Lu — Aug [1] 4. Ordnung 
2. | {Zu.+ 4 2up — 20 u — -; Z2Aup +8 Aud | [1] 6. Ordnung 
3. ? (Fu, " va % >* = N? {3 Aus +42 Au + 52 Ana} [2] I 6.Ordanie 
AN Nu A — Wu — R?{& Au. + DAup + 82 Auas — a ht? AAua In 8. Ordnung 


> 1 I 
5. | Ze. + 42m —-20u + TER h?{2Aup— 34 Auar — 30 ht{2 AAup + 17 AAugN [3] 8. Ordnung 


| — 0 


6. Due — 12 Zu. — 64 Dun + 300 ua} — ee h? {11 2 Au. + 16 2 Aun — 3468 Au} 10. Ordnung 
un 335 h?{& AAue + 164 2 AAup + 2040 AAduu} 

7. | {£u0 +2 Eu, — 8 Fun + 20 ua} — hr Aug a 

8. Zy + Lu — 22u — 10 LZuy + 36u — = he {2 AAun + 2AAug} [4] | 8. Ordnung 


9. | (Zu, +8 Zu + 18 Zu. — 8 Zu, — 144 Zup + 468 u) + Se ht {2 AAu. — 26 2AAu. 
— 72 2AAup — 332 AAu.} [4] | 10. Ordnung 
10. | {- Zus —3 Fu +12 Zu. + 24 Zu. — 57 Zu + 112 u0} 


1 
mern {2 AAu. +2 2 AAup — 12 AAug} 10. Ordnung 


B. Dreieeksnetze 


Die zu Anfang des Abschnitts A ge- 
gebenen Erklärungen für die Bezeichnungen 
Us Du; - . . usw. gelten sinngemäß auch für 
Dreiecksnetze, wobei jedoch die Lage der ver- 
wendeten Punkte jetzt aus Bild B hervorgeht. 
Schema B enthält die allgemeine Taylor- 
Entwicklung; die darin auftretenden Aus- 
drücke Au und ARu(n=0,1,2,...) sind 
wie folgt erklärt 


i+k 
Z bedeutet an) 
[OR 


a) 


oyk 
Mar == Ugo — 15 Wut 15 Ua — Uns 
\ Au = Ugg — 14 Up + 14 Uyg — Ugs 
j- Au = ron — 13 Ug — 14 ug + 14 U + 13 Uns — Ugıo 
; Rt er 
w bzw. 


Au — Ungo— 66 Ua + 495 Ugg — 924 Uge + 495 Ugg — 66 Ugio + one - 
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Schema B 


Au +b Fun +e2u | +dug 
2 
+ 5r {AAus +BEAw +0 FEAu +DE Aug 


- h* 
ai Fr. {& AAua + BZ AAun + y Z AAu. +6 % AAug} 


Re ee ee ARE EEE 
= Ua 6 + c+ d-a 

h? 
+-rAua | 3 B+ 4c+ 394 + Er am 

BR ETNB 4\ 
+ Au | b+ 16c+ 99+(5) HB+ 4c+ 3D+ BR rer 


+ EA Mr ne 640 + 74+(9 2) 8 + 160 + 9m +(4 SB, gar ar 
+ Are 9 ®+ 256c + sıd-+(> Br 640 + 27 +(4 IE 169 + 9 + 
Een E B + 10246 + 2430) 43 + 2560 + aD +[% )Ee+ ++] 
+ At = + 40966 + 7299 + (% )tB + 10220 +213D)+(% EB +20 +1 +] 


+ At 86+ 64c — 27) | 


en 


m 64 2560 — sıa+(3 RE: 640 — 27D}| 
Alte [320 + 10240 — 239 + (5) SB+ »00- um+(Y) + 6y- 20] 
+2 Apus 38 + 4096: - 1294 + (2) EB +10240-243D) +4 1 +2586y— 819 +--] 
+ = Alla 5:8 + 4096 c + 729 | 

© L 


+... Glieder 14. und höherer Ordnung. 


Zur einns finiter Ausdrücke für Au = r(x, y) bestehen folgende Möglichkeiten (siehe 


12 
inder Taylor- Entwicklung die Koeffizienten von N Au, 2 Su, a 29,U,): 
a=6 b=—1 Restglied 6. Ordnung erreichbar 
0168 b= —27 d=—1 ran 12. ” % 
a = 3844 b=—624 d=—16 et 19; dr f 


Von den in Tabelle B aufgeführten finiten Ausdrücken ersetzen die unter Nr. 1—4 genann- 
ten die Differentialgleichung Au+ Au =0 durch Gleichungen der Form Nu, =» u, [2] (‚Spezielles 
Matrizeneigenwertproblem‘‘ [2]); dabei wird: 


Formel Nr. | 


48 +v— 684 + 366» +»? 
30 


46 Albrecht, Taylor-Entwicklungen und finite Ausdrücke für Au u. AAu ga. 33 Re le Se seeBsKinen 


Tabelle B 
Nr. Formel Restglied 
= 3 Mm 
1. | {Eup — bu} — =y h? Aug [1] a 2a Un +: 
1 ni 75 ER 
en eh Te he, ie BE 
2.| {Zu — 6uu} 16 h?{& Aup + 18Aua} 1] ei 354 Uq +yg 23 ta) + 
s 6 
3.| Dun — bu — I 1 Aua— 8 eh = 18 Ara + 1 Aa) En 
4.\ır ; r WEITER 
. |) {Fun — 6Uat — — gs" {2 Aup + 48 Auu} — 6 h2AAug Eon TR Au + 


5. | {Zua + 27 Zup — 168uu} + Li h? {X Aug — 21 Z Aup — 240 Au} | 8. Ordnung 


6 } 9 h® 45 

E Fu — 3 Fun + 12u a} — Tg hAAUa [1] = 145.40 Ab un 4... 

” 1 h® 45 

T. {Zu — 3 2up + 1200} — ZH AAu, + 15 AAua} = 9754... © 


8.| {Zu +2Fug — 10 2un +42 u — ER HZ AAup+6A Au.) [1] | 8. Ordnung 
16 | 


C. Sechsecksnetze 


Wegen der Bezeichnungen verweisen wir wieder auf die Erläuterungen zu Anfang von Ab- 
schnitt A, wobei jetzt Bild C zugrunde zu legen ist. Derinder Ta ylor- Entwicklung (Schema C) 
auftretende Operator A% wurde bereits in Abschnitt B erklärt; wir führen noch ein: 

AzU = Ugg — Auge — I Una. 


Schema C 
NR UtrbLu +cLu +dZua 


ara Aug + BE Au +0 Z Aue + DZ Aus} 


en u|36+20+ tal 


! 
5 Aal u: a a HS 2C+D)+4| 
Ana | + 180 + 169 + (5) 18 + 6C+ 1) C 
n 1 
Adna p+51c+ 09 +(2) 5 1B+18C+ 162} | 
s 3 
I Aa IH ” sa| 
5 15 53 
Aa Ei 32 +(5)—- {B R% 82} | 
ns 3 
Bild 0 + Abt gb -Me+6g| 


-+... Glieder 7 und höherer Ordnung 


Einige einfache finite Ausdrücke sind in Tabelle C zusammengefaßt. — Die Sechsecksnetze 
besitzen gegenüber den Dreiecksnetzen den Vorteil, bei gleichem Grundbereich und gleicher 
Maschenweite 3 weniger Punkte aufzuweisen, so daß zur Lösung eines Problems auch weniger 
Unbekannte einzuführen sind; beim Aufschreiben von finiten Gleichungen derselben Ordnung in 
randnahen Stellen hingegen erweisen sich die Sechsecksnetze als ungünstiger, weil die heran- 
gezogenen Gitterpunkte weiter von der Entwicklungsstelle entfernt sind als in Dreiecksnetzen 
(vgl. z. B. die Formeln 6. Ordnung). 
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Tabelle € 
Nr. Formel | Restglied 
18 {Zup — 3Uua — Ip Ava [1] 3. Ordnung 
2. {Zua +8 up — 27 ua} — 9 h?Aua | 4. Ordnung 
82 ) {Zu +8 Eu — 27 ua} — er h?{& Augs +22 Aup + 27 Aua} | 6. Ordnung 
37% {Zus —22u. +82 — lu} + - h* Adua | 5. Ordnung 


D. Würfelnetze 


Bezeichnungen:w,seider Funktionswert in der Mitte eines Würfels von der Kanten- 
länge 2h; Nu,, Nu, und u, jeweils die Summe der Werte in der Mitte der Seiten, in der Mitte 
der Kanten und in den Ecken des Würfels. Entsprechendes gelte für Au, N Au,, .. . usw. 


e) EEE 


a — — DT er 
= +3 tra 02? (n ’ E) > ) Ujkı = 92 Oyk 92! 9! 
Us,2,0 = He Unge + Yoge ; Ug,a,o — Unso 4 Waoa + os - 
Schema _D 


M* Um 48 2us + k Zur teLu. ® 
+ {MM - Aum +SEAu +KZAur +EE Auf} 


{—8s+ 80k+ 192 + (5) 68 + 36K + 96.2 + (4)t- 40+ 8% + 328} 


{ 68—156k+ 3840)| 


+ {m Adum +0 EAAus + # 2 AAup + © EAAu) 
NE RE, SA BL EEE NR ARE 
a E 6s+ 12k+ Be +m}| 
ST Am E 2s+ 8k+ 83+ { 68$+ 2K+ 8E + | 
+ Arm & I 89+())t 298+ 8K+ 8H+ {o+12x+ Bc+p| 
€) + um er Beach sy +(2)t 25+ 8K-+ sm +(42e+ 8xt 8-+- 2 
+ 5 Arm E 2:4 8k+ 89 +(2)t 39. SEK sm +(4)20+ SCH 8 +: 
Hamm [| 4e+. 84 229] 
FE Auya,om [#804 36K + 960 + (9) 48 + 8K-+ 32E} 
| 
| 


8 
+ gr Aumm { 166 — 496k + 2304 9 + (2) 68 — 156 K + 3842) | 


h® 
+ gr laaom l& 4s+ 72k+ 128 9 | 


+... Glieder 10. und höherer Ordnung 
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Aus der in Schema D gegebenen allgemeinen Taylor - Entwicklung lassen sich die in 
Tabelle D genannten finiten Ausdrücke herleiten; wir erwähnen auch hier die zur bequemen 
Lösung von Eigenwertproblemen Au +4 u = 0 besonders geeigneten Formeln: 


Formel Nr. Finite Gleichung h? 3} 
3 | 57 re 12 24 —y 
| Zus +2 2U=V* Um 184» 
p | „56 —v 
5 | Zue +8 2u= 7: um 12 505 
Tabelle D 
Nr. Formel Restglied 
1.| {Zus — 6um} — h? Aum [1] | 4. Ordnung 
2.| {Zur +2 Zus — dumm — nn h>i& Aus + 6 Aum} [1] | 6. Ordnung 
3. (Zur + 2 2u — 24 um} — 4 {2 Au +2 2% Aus + 48 Aum} 6. Ordnung _ C 
4.| {Zus 48 Zus — 56um} — hr {Z Aus + 6 Aum} 6. Ordnung 
5.| {Zu +8 Zus — 56 um} — n {2 Aue + 8% Aus + 88 Aums 6. Ordnung 
6.) Zu. +3 Fur +14 Zus — 128 um} + B{Z Aus — 36.Aum} — h{Z AAus + 15 AAum}| 8. Ordnung | 
TI ES Dir MA Zu 18 . h? {Aug — ZAus 84 Aumt — = h* Adum 8. Ordnung 
8| Zu. +3 Fur +14 Fu — 128 um} — — h? {2% Aug +2 2% Aus + 160 Aum} — = h% en 8. Ordnung 


Die Göttinger elektronischen Rechenmaschinen 


Die Göttinger Arbeiten über elektronische programmgesteuerte Rechenmaschinen werden im I. Teil 
in ihrem Verhältnis zur neueren Entwicklung mathematischer Maschinen charakterisiert. Im II. Teil 
wird die bandgesteuerte Rechenmaschine im einzelnen beschrieben und ihre Arbeitsweise erläutert. Sie ist 
mit einem Aufwand von. knapp 500 Röhren und 100 Relais erstellt. 

x - In part I, the work that was done in Göttingen in the field of electronic sequence controlled computers, : 
is characterized in its relation to the recent development of mathematical machines. In part II, the tape- C 
controlled machine is described in detail, and its working principles are explained. T'he computer con- 

y tains less than 500 vacuum tubes and 100 relays. 

Dans la Iere part les travaux de Göttingen sur des machines & calcul £lectroniques dirigees suivant 

un programme fiwe sont caracterises dans leur relation au developpement present des machines math&matiques. 
Dans la IIe part la machine & calcul dirigee par ruban est decrite en d&tail, et su maniere de fonction est 
expliquee. Elle est construite avec un appareil d’un peu noins de 500 tubes et d’un peu moins de 100 relais. 

B nepsoä yacrı xaua XapakTepmcTuka HPOBOAHMBIX B Terrunreue pa0oT B OÖNacTH Blnek- 
TPOHHEIX IPOTPaMMHO-YIPABIAeMbIX CYeTHBIX MAMIHH IIO OTHOIMEHHIO K HOBeÄNIeMy PA3BHTHIO 
MATEMATUYeCKUX MalmuH. BoBTOpoü YacTH NeTAIBHO OUUCAHA NEHTOYHO-YIPABIAEMAA CYETHAA 
MAIIUHa H OÖBSICHAeTcH ee cIocoO pa6ors. Mamma aTa NocTpoeHa e 3aTpaToi MeHee 
500 „ıamı m 100 pen. 


I. Allgemeines } 


Von L. Biermann und H. Billing 
(Max-Planck-Institut für Physik und Arbeitsgruppe numerische Rechenmaschinen im Institut 
für Instrumentenkunde in der Max-Planck-Gesellschaft.) 
Unter den mathematischen Maschinen lassen sich gegenwärtig zwei Haupttypen unter- “ 
scheiden, nämlich die nach dem Analogieprinzip arbeitenden Geräte, deren wichtigstes Beispiel e 
die Integrieranlagen darstellen!), und die logisch-mathematischen Maschinen, von denen be- a 


; !) In Deutschland besitzt das Institut für praktische Mathematik der TH Darmstadt (Prof. A. Walther) Er 
eine ausgebaute Anlage dieser Art. Te) L 
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sonders die elektronischen Rechenmaschinen bekannt geworden sind. Beide Typen gehören in 
den weiteren Bereich der Kybernetik im Sinne von N. Wiener). 

Innerhalb der zweitgenannten Gruppe kommt neben den Ziffernrechenmaschinen den 
logistischen Maschinen (im engeren Sinn) in folgender Hinsicht eine besondere Bedeutung zu. 
Die bisherigen Konstruktionen?) zeigen bereits, daß es mit ganz einfachen Mitteln möglich ist 
den Aussagenkalkül der mathematischen Logik abzubilden und nachzuvollziehen, also, wenn 
eine etwas unscharfe Ausdrucksweise erlaubt ist, ungefähr den Bereich der aristotelischen Logik. 
Vorarbeiten, die über diesen Kalkül hinausführen, sind an verschiedenen Stellen geleistetworden®). 
Es ist sicher, daß man mit den heutigen Mitteln weit über den Bereich des Aussagenkalküls 
hinausgehen kann, und daß — mit gewissen Einschränkungen — auch der Prädikatenkalkül 
und damit der weitere Bereich der mathematischen Logik entsprechenden Geräten zugänglich 
ist. Soweit sich also Teilgebiete der Mathematik als logisch-deduktive Formalismen aulfassen 
lassen, deren Teilaussagen durch kodifizierte Regeln aus den Grundannahmen folgen, werden 
sie auch der Abbildung und dem Vollzug durch Maschinen zugänglich sein?°). 

‚Diese Bemerkung macht deutlich, daß die Ziffernrechenmaschinen — trotz der Möglichkeit, 
auch allgemeine Probleme mit ihnen anzugreifen — nur eine spezielle Form logisch-mathe- 
matischer Maschinen darstellen. Sie verdanken ihre bevorzugte Entwicklung neben der praktischen 
Bedeutung einiger der numerischen Analysis zugänglichen Probleme vor allem der Tatsache, 
daß die modernen technischen Mittel eine ins einzelne gehende und sich selbst automatisch 
modifizierende Steuerung des Arbeitsprogramms mittels umfangreicher Speicher für Befehls- 
folgen und Zahlen bei sehr hoher Arbeitsgeschwindigkeit ermöglicht haben. Es ist aber sicher, 
daß die gleichen Mittel der maschinellen Bearbeitung auch ganz andersartigen Aufgaben (z.B. 
analytischer Art) zugute kommen würden, deren Lösung durch die große Zahl der eingehenden 
Bestimmungsstücke oder der notwendigen Operationen erschwert’ wird. 

Innerhalb der praktischen Analysis ergibt sich im Zusammenhang mit der Entwicklung 
programmgesteuerter Geräte die Aufgabe, die Rechenmethoden den neuen Möglichkeiten an- 
zupassen und sie in dem für jede Aufgabe aufzustellenden Rechenprogramm zu realisieren. Es 
ist leicht, Beispiele dafür zu geben, daß bisher übliche Näherungsverfahren nicht mehr zweck- 
mäßig sind, da sie die Vorteile der Programmsteuerung nicht auszunutzen gestatten, während 
andere Verfahren, deren praktische Durchführung bisher etwa an dem Umfang der Rechenarbeit 
scheiterte, jetzt neue Möglichkeiten eröffnen®). Die Entwicklung der mathematischen Verfahren 
und die technische Entwicklung müssen daher Hand in Hand miteinander vor sich gehen. 

Die in den nächsten Arbeiten dieser Folge beschriebenen Geräte, die von der Arbeitsgruppe 
numerische Rechenmaschinen im Institut für Instrumentenkunde in der Max-Planck-Gesell- 
schaft in enger Zusammenarbeit mit der astrophysikalischen Abteilung des Max-Planck-Instituts 
für Physik entwickelt worden bzw. noch in Entwicklung begriffen sind, gehören in den Bereich 
der programmgesteuerten elektronischen Rechenmaschinen. Als erstes wurde nach längeren 
Vorarbeiten, die sich z. T. auch auf Geräte ganz andern Typs erstrecken, um 1949/1950 mit dem 
Bau einer größeren mittelschnellen Maschine begonnen, welche als Speicher eine Magnettrommel 
besitzt, die schon früher entwickelt worden war”). Die Trommel ist eingerichtet zur Aufnahme 
von 2048 Dualzahlen zu je 32 Stellen (äquivalent 10stelligen Dezimalzahlen). Ein beliebiger 
Teil der Zahlenspeicher kann auch zur Aufnahme der Rechenkommandos dienen. Die Maschine 
kann die Rechenkommandos auf Grund vorangehender Rechenergebnisse selbständig abändern, 
so daß sie die Flexibilität der großen modernen Rechenautomaten haben dürfte. Die Fertig- 
stellung dieser Maschine, deren Bau ziemlich weit vorgeschritten ist, ist für 1953 vorgesehen. 

Während der Planungsarbeiten an der eben erwähnten Maschine ergab sich die Frage nach 
der Nützlichkeit einer kleineren programmgesteuerten Maschine. Der Abstand zwischen der 
Leistungsfähigkeit größerer programmgesteuerter Maschinen und der einer normalen Tisch- 
Rechenmaschine ist so groß, daß es sinnvoll schien zu untersuchen, wo die untere Grenze des 
elektronischen Aufwandes liegt, mit dem bereits eine größenordnungsmäßige Überlegenheit 
erzielt wird. Ebenso sollten in noch höherem Maße der Betriebssicherheit dienende Rücksichten 

2) N. Wiener, Cyberneties, Control and Communication in the animal and in the machine. Cambridge, 
Mass., 1949. The Human Use of Human Beings, London 1950. 

3) Durch K. Zuse (1944 unveröffentlicht, vgl. folgende Anmerkung) sowie D.M. McCallum und 
J.B.Smith, Electronic Engineering 1951 p. 126 u. 458. 

4) 8.z.B.K.Zuse, Archiv für Mathematik 1, 441 (1948). Herrn Zuse danken wir für die Möglichkeit 
in unveröffentlichte Arbeiten über diesen Gegenstand Einsicht zu nehmen. 

5) Eine Bearbeitung des hierdurch angedeuteten Problemkreises (speziell der den Maschinen gesetzten 
Grenzen) auf der Grundlage der Theorie der „berechenbaren“ Funktionen (im Sinne von A.M. Turing, 
Proc. London Math. Soc. 42, 230 (1937), s. auch R.Peter, Rekursive Funktionen. Budapest 1951) ist in 
Angriff genommen. 

®) Im Bereich der Quantenmechanik der Moleküle hat A. v. Mohrenstein eine solche Untersuchung hier 
durchgeführt. 


2 ?) H. Billing, Numerische Rechenmaschine mit Magnetophonspeicher, Z. angew. Math. Mech. 29 


(1949) 8. 1. 
4 


r 
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genommen werden, und die Programmfertigung sollte möglichst einfach durchführbar sein. Es 
zeigte sich, daß mit einem Aufwand von knapp 500 Elektronenröhren und einer kleinen Magnet- 
trommel für etwa 10°—10° zehnstellige Dezimalzahlen allgemein gesprochen etwa die Leistungs- 
fähigkeit mittelgroßer Relais-Maschinen, wie sie jetzt noch an mehreren Stellen in Gebrauch 
oder im Bau sind, erreichbar sein dürfte. Allerdings mußte hierbei auf die Möglichkeit, die 
Befehle während des Ablaufs der Rechnung automatisch zu ändern, verzichtet werden. 

Die nun fertiggestellte Maschine (s. die folgende Arbeit) ist demgemäß kein verkleinertes 
Abbild der noch im Bau befindlichen größeren, sondern eine eigene Entwicklung auf Grund der 
eben dargestellten Gesichtspunkte. Beide Maschinen arbeiten im reinen Dualsystem mit auto- 
matischer Übersetzung aus dem bzw. ins Dezimalsystem bei der Eingabe bzw. Ausgabe. Die 
arithmetischen Grundoperationen werden bei beiden Maschinen in Serienform ausgeführt. Das 
bedeutet im wesentlichen, daß bei Additionen die einzelnen Ziffern der Summanden zeitlich 
nacheinander von der niedrigsten» Stelle beginnend addiert werden. Bei. Parallelmaschinen 
werden alle Ziffern des Addenden gleichzeitig addiert und bei der Addition anfallende Überträge 
nachträglich behandelt. 


II. Die Göttinger bandgesteuerte Rechenmaschine 6 1 
Von A. Billing, W. Hopmann, A. Schlüter 


(Arbeitsgruppe numerische Rechenmaschinen und Max-Planck-Institut für Physik, Göttingen.) 


Übersicht über Aufbau und Leistung der Maschine 

Die wesentlichen Einheiten der bandgesteuerten Rechenmaschine sind (Bild 1): 

1. Der statische Speicher. Er dient zur Aufnahme von numerischen Eingangswerten und 
Zwischenresultaten und besteht aus einer rotierenden Trommel, deren Oberfläche mit einer 
magnetisierbaren Schicht bedeckt ist. 

2. Das eigentliche Rechenwerk mit dem Akkumulator und den beiden Hilfsregistern Md 
und Mr. Im Akkumulator fallen die Resultate der Rechenoperationen an, während die beiden 


Hilfsregister zur Aufnahme der bei der jewei- 
N ligen Rechenoperation benötigten Operatoren 
Verteiler 


Befehlsw, 


— bei der Multiplikation z. B. Multiplikand 
und Multiplikator — dienen. 

3. Der Verteiler. Dieser dient einerseits 
als Verbindungsglied zwischen statischem Spei- 
cher und Rechenwerk. Andererseits wird in 
ihm zur Aus- und Eingabe die Umwandlung von 
Dezimal- in Dualzahlen und umgekehrt durch- 
geführt. 

4. Das Ein- und Ausgabewerk. Es dient 

Bild 1. Aufbau der Maschine zum Eingeben von Zählen und Rechenbefehlen 

in die Maschine und zum Drucken der Resul- 

tate. Es besteht aus einer Schreibmaschine, deren Tasten zur Eingabe mit Kontakten ver- 

bunden sind und zur Ausgabe von Elektromagneten betätigt werden können, und aus einem 
Satz von Lochbandsendern zur automatischen Eingabe. 

5. Das Befehlswerk. Dieses ist das eigentliche Steuerorgan. Es bringt mit Hilfe eines 
Schaltwellengenerators die vom Eingabewerk einlaufenden Rechenbefehle zur Ausführung. 

6. Das Synchronisationswerk. Es sorgt für die richtige zeitliche Abstimmung der in den 
einzelnen Maschinenelementen auszuführenden Operationen aufeinander. Es ist starr mit der 
Speichertrommel gekuppelt, so daß Drehzahlschwankungen der Trommel sich nicht störend 
auswirken. 


Die Rechengenauigkeit der Maschine beträgt 32 Dualziffern, das entspricht fast 10 Dezimal- 
ziffern. Das Rechenwerk rechnet mit festem Komma, welches hinter die erste Dezimalziffer 
gelegt wurde. Jedoch darf die erste Dezimalziffer nicht größer als 7 sein. Der Rechenbereich 
liegt also zwischen +7,999999999. Bei Überschreiten des Rechenbereichs ertönt eine Warnung 
bzw. beim Arbeiten mit Programmsteuerung wird der Ablauf der Rechnung unterbrochen. 
Sämtliche Rechenoperationen werden in Serienform durchgeführt. Dies bringt gegenüber 
Parallelmaschinen neben einer erheblichen Verringerung des technischen Aufwandes den Vorteil, 
daß bei Ausfall eines einzigen Schaltelementes die Rechnung praktisch immer völlig verdorben 
wird, so daß der Ausfall sofort auffällt. 

Das Rechenwerk der Maschine enthält spezielle Einheiten für Addition, Subtraktion, 
Multiplikation, Division und das Ziehen der Quadratwurzel. Bei der Multiplikation wird zu- 


nächst das 20stellige Produkt errechnet, aber anschließend automatisch auf 10 Stellen abgerundet. 
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Auf Vermeidung systematischer Abrundungsfehler wurde großer Wert gelegt. Einzelheiten über 
die Rechenkommandos, wie über die Dauer ihrer Durchführung werden im Abschnitt ‚‚Befehls- 
plan‘ besprochen. 

Der statische Speicher der Maschine enthält Platz zur Aufnahme von 26 Zahlen zu je 
32 Dualziffern. Jeder dieser 26 Zahlenspeicher ist im Laufe der Rechnung auf Kommando hin 
nach spätestens 1/,, Sekunde aufrufbar zur Entnahme bzw. Eingabe. Bei der Eingabe in einen 
bereits von früher her gefüllten Speicher wird der vorherige Inhalt automatisch gelöscht. Der 
Speicherinhalt bleibt auch bei Abschalten der Maschine erhalten. Als zusätzliche Zahlenspeicher 
können die ohnehin zur Programmsteuerung vorgesehenen Lochbänder verwendet werden, da 
es möglich ist, Zwischenresultate in ein Lochband zu lochen und sie auch wieder vom Lochband 
automatisch abzulesen. Das Ablesen vom Lochband ist aber nur in der Reihenfolge möglich, 
in der die Zahlen gelocht sind. Wenn man die einzugebenden Zahlen auf verschiedenen der 
insgesamt 4 Lochbänder unterbringt, bleibt noch die Wahlmöglichkeit, welches der Lochbänder 
man jeweils ablesen will. 

Die Eingabe von Zahlen bzw. von Rechenkommandos ist auf zwei verschiedene Weisen 
möglich. 

1. Eingabe vom Tastenfeld einer Schreibmaschine aus: Zur Eingabe von Zahlen und von 
Rechenkommandos, wie zum Ausschreiben von Resultaten wird eine handelsübliche Schreib- 
maschine mit elektrischem Antrieb (Mercedes Electra) verwendet. Unter den Tasten dieser 
Maschine sind zusätzlich Kontakte angebracht. Nach Anschlagen der gewünschten Taste wird 
einerseits wie bei jeder Schreibmaschine das angeschlagene Zeichen auf das eingespannte Pro- 
tokollpapier geschrieben. Andererseits wird durch Schließen des mit der Taste verbundenen 
Kontaktes die zugehörige Ziffer bzw. das zugehörige Rechenkommando der Rechenmaschine 
zugeleitet und löst dort durch Betätigen eines Relais die gewünschte Operation aus. 

Dieses Verfahren bietet den Vorteil, daß beim Rechnen vom Tastenfeld aus das zugehörige 
Rechenprotokoll automatisch mitgeschrieben wird. Um Rechenfehler auszuschließen, die durch 
zu schnelle Aufeinanderfolge von zwei Tastenanschlägen auftreten könnten, wird das gesamte 
Tastenfeld während der Ausführung eines Befehles blockiert. Dies wirkt sich praktisch nur bei 
den Befehlen Multiplikation, Division und Wurzelziehen aus, da sämtliche anderen Befehle in 
weniger als Y/,, Sekunde ausgeführt werden. 

Die Maschine rechnet stets mit zehnziffriger Genauigkeit. Um die Eingabe von Zahlen zu 
erleichtern, werden für Rechnungen mit kleinerer Genauigkeit die ‚„‚nicht bedeutsamen“ Nullen 
am Ende einer Dezimalzahl automatisch angehängt. Dies wird dadurch erreicht, daß das Vor- 
zeichensymbol ‚‚positiv oder negativ“ (p oder n) am Ende der Zahl eingegeben wird. Mit der 
Eingabe des Vorzeichensymbols ist das Anhängen der nicht bedeutsamen Nullen automatisch 
gekoppelt. Anschlagen der Tastenfolge 783p bedeutet mithin Eingabe der Zahl +7,830000000. 
Um das Rechenprotokoll übersichtlich zu gestalten, wird nach Eingabe jeder Zahl (und zwar 
nach dem Vorzeichen) wie nach Eingabe logisch zusammenhängender Befehlsgruppen (siehe 
Tabelle 1) automatisch ein Zwischenraum im Protokoll gelassen. 

2. Eingabe vom Befehlsband aus: Der Einsatz elektronischer Rechenmaschinen zur Durch- 
führung längerer Rechnungen lohnt sich nur, wenn sich der Rechengang so gestalten läßt, daß 
eine im Vergleich zur ganzen Rechnung kurze Folge von Rechenbefehlen immer wieder auf neue 
Eingangswerte bzw. auf bei der Rechnung angefallene Zwischenresultate angewandt wird. Ein 
Musterbeispiel für derartige Rechnungen ist die Iteration. 

Zur Eingabe sich wiederholender Befehlsfolgen werden Lochbänder benützt, die an den 
‚Enden zu einem endlosen Band zusammengeklebt werden können. Zur Vereinfachung der Kon- 
struktion wurde das handelsübliche in Fernschreibern verwendete Lochband eingesetzt, welches 
in jeder Zeile bis zu 5 Löchern enthalten kann. Dies ergab den Vorteil, daß für die Herstellung 
wie für das Abtasten der Lochbänder handelsübliche Geräte der Firma Lorenz nach unbedeu- 
tenden Abänderungen unmittelbar verwendet werden konnten. Die Anzahl von 5 Lochmöglich- 
keiten in einer Zeile gestattet die Darstellung von 2? — 32 verschiedenen Zeichen. Eine hinter 
dem Abtaster des Lochstreifensenders geschaltete Relaispyramide mit 5 Eingängen und 32 Aus- 
gängen bewirkt, daß entsprechend der jeweiligen Lochkombination in der abgetasteten Zeile 
der zugehörige Ausgang der Pyramide unter Strom gesetzt wird. Dieser Strom betätigt das 
gleiche Relais, das bei Arbeiten von der Tastatur aus sonst bei Drücken der zu dem betreffenden 
Befehl gehörenden Schreibmaschinentaste erregt wird. 

Der Lochstreifenabtaster tastet normalerweise 7 Zeilen pro Sekunde ab. Die meisten 
Operationen können von der Rechenmaschine in der Zwischenzeit ausgeführt werden. Bei 
Operationen, die längere Zeit benötigen, wartet der Lochstreifenabtaster automatisch bis zur 
Beendigung der Ausführung. 

Bei größeren Rechnungen ist ein einzelnes Lochband nicht ausreichend, vor allem, wenn 
ein Lochband noch als zusätzlicher Zahlenspeicher benutzt werden soll. Einbau von insgesamt 
4* 
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4 Eingabelochbändern ist vorgesehen. Bei endlos zusammengeklebten Befehlsbändern;kann man 
das Band durch einen eingelochten Bandbefehl am Ende jedes Umjaufs anhalten lassen. Der 
Rechner kann dann auf Grund der während des Befehlszyklus gedruckten Resultats entscheiden, 
ob der Zyklus noch einmal wiederholt werden oder zu einem anderen Befehlsband übergegangen 
werden soll. Auf Einbau eines konditionellen Entscheidungsbefehles wurde verzichtet. Derartige 
konditionelle Befehle werden in den großen Rechenautomaten verwendet, um entsprechend dem 
im vorangehenden Rechengang errechneten Vorzeichen einer Zahl einen Befehlszyklus zu wieder- 
holen oder im Hauptprogramm fortzufahren. Bei der relativ niedrigen Rechengeschwindigkeit 
dieser Göttinger Maschine erschien es nicht lohnend, dem Rechner diese Entscheidung mittels 
eines Automaten abzunehmen. 2 Glimmlampen zeigen dem Rechner laufend das Vorzeichen der 
‘ im Akkumulator befindlichen Zahl an, so daß es zum Treffen der Entscheidung nicht nötig ist, 
den Akkumulatorinhalt für diesen Zweck auszudrucken. Der Übergang zu einem anderen Band 
kann auch automatisch mittels eines in das erste Band gelochten Befehls bewirkt werden. Dieses 
ist z. B. wertvoll, wenn man einen speziellen Befehlsgang mehrere Male an verschiedenen Stellen 
in den Hauptbefehlsgang einbauen will. 


Die Herstellung der Lochbänder kann von der Schreibmaschinentastatur der Mercedes- 
Electra aus erfolgen. Zu diesem Zweck sind die 32 von der Schreibmaschine her erregbaren 
Eingaberelais über eine Gleichrichterpyramide mit 32 Eingängen und 5 Ausgängen derartig mit 
den 5 Eingängen eines Bandlochers verbunden, daß bei Niederdrücken einer Schreibmaschinen- 
taste unmittelbar die entsprechende Lochkombination gelocht wird. Auf ganz ähnliche Weise 
können auch errechnete Resultate — zusätzlich zum Ausdrucken — in das Band gelocht werden 
Jede Zeile im Lochband entspricht dabei einer Dezimalziffer. 


Diese Herstellungsart eines Befehlsbandes hat den Vorteil, daß 


1. das Protokoll für das Befehlsband automatisch mitgeschrieben wird und 


2. von der Rechenmaschine während der Lochbandherstellung die einzelnen befohlenen 
Operationen gleichzeitig ausgeführt werden, so daß man bei Beendigung der Lochbandherstellung 
einen vollen Befehlszyklus mitgerechnet hat. Das kann für die Kontrolle der Richtigkeit des 
Befehlsprogrammes sehr nützlich sein. 


Ein Nachteil dieser Herstellungsart besteht darin, daß die Rechenmaschine für die Dauer 
der Lochbandherstellung nicht für andere Zwecke benutzt werden kann. Durch Verwendung 
eines zusätzlichen Handlochers, der unabhängig von der Rechenmaschine zur Lochbandher- 
stellung benützt werden kann, ließe sich dieser Nachteil leicht vermeiden. 


Zyklische Vertauschung von Speicherplätzen: Um die Zerlegung eines Rechenganges in 
eine sich vielfach wiederholende gleichbleibende Befehlsfolge zu erleichtern sind 16 der statischen 
Speicher in 4 Gruppen zu je 4 Speichern aufgeteilt. Diese Speicher sind bezeichnet mit «a0, al, 
a2, a3, 60---b3; c0---c3; dO---d3. Durch einen gesonderten Befehl: „‚Zyklische Vertauschung‘“ 
kann der Inhalt eines jeden Speichers einer Gruppe in den Speicher mit der nächst höheren 
Nummer derselben Gruppe verschoben werden, wobei die Nummer 3 +1 identisch 0 ist. Es 
wird also bei dieser zyklischen Vertauschung folgende Verschiebung der Speicherinhalte bewerk- 
stelligt (a0) > al <aly>a2 <a2y > a3 <a3> >a0. <a0) bedeutet Inhalt des Speichers «a0 
vor der. Vertauschung. 


Eine naheliegende Anwendung ergibt z. B. die punktweise Integration einer gewöhnlichen 
Differentialgleichung. Die Lösung sei bekannt für die Punkte, ©,_ı, %,. Die zugehörigen Lösungs- 
werte seien in die Speicher b1, bO gegeben. Daszur Lösung der Differentialgleichung angewandte 
numerische Verfahren ergibt das Rechenprogramm, um aus den Lösungen für die Abszissen 
%n—1, %, die Lösung für die Abszisse @,.;ı zu berechnen. Diese Lösung wird nach 53 ein- 
geschrieben. Wenn jetzt der Befehl ‚‚zyklische Vertauschung“ folgt, stehen hinterher die Lösungs- 
werte für 24-1, %, @n+ı in den Speichern 52, bl, 50. Man kann jetzt genau den gleichen 
Befehlsgang ein zweites Mal anwenden, um die Lösung für @,;2 zu berechnen usw. Im Anhang 
ist ein derartiges Beispiel genauer ausgeführt. 


Zum Drucken von Resultaten wird die Schreibmaschine verwendet. Unter den Tasten 
für die Ziffern 0 bis 9 wie für die Vorzeichensymbole p und n sind Elektromagnete angebracht, 
welche von der Rechenmaschine her zum Ausdrucken der in Dezimalform umgewandelten 
Resultate erregt werden. Nach jeder ausgedruckten Zahl wird automatisch ein Zwischenraum 
geschrieben. Zum Ausschreiben von Tabellen ist ein Zahlenzählwerk eingebaut, welches nach 
einer einstellbaren Anzahl von in einer Zeile ausgeschriebenen Zahlen (bis zu 20) den Wagen- 
rücklauf der Schreibmaschine betätigt. Die Anzahl der pro Zahl auszuschreibenden Ziffern ist 
von Hand aus willkürlich einstellbar (1—10 Ziffern). Dabei muß jedoch berücksichtigt werden, 
daß die Zahl beim Ausschreiben nicht abgerundet, sondern einfach abgebrochen wird, so daß 


man zweckmäßigerweise die Anzahl der auszuschreibenden Ziffern nicht zu sparsam bemißt, 
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Dies ist die einzige Stelle in der Rechenmaschine, wo auf Abrundung verzichtet wird. Da einmal 
ausgeschriebene Zahlen jedoch kaum vielfach in einem Befehlszyklus wieder verwendet werden, 
dürfte der hierdurch entstehende systematische Fehler keine Rolle spielen. 


Hauptteile der Maschine 


Der statische Speicher dient zur Aufnahme von 26 in Dualform umgewandelten Zahlen 
zu je 32 Dualziffern.- Die Anzahl von 26 Speicherplätzen ist ein wenigknapp. Eine Steigerung 
der Speicherkapazität auf 100 oder 1000 Speicherplätze würde bei einem Neubau der Maschine 
ohne merklich erhöhten Aufwand möglich sein. Der Speicher ist technisch realisiert durch eine 
mit 50 U/sek rotierende zylindrische Trommel, die mit einer magnetisierbaren Schicht bedeckt 
ist. Ihre Oberfläche rotiert an Elektromagneten vorbei, die zum Einschreiben der Zahl wie zum 
Ablesen benützt werden. Ein zu einem einzelnen Schreibmagneten gehörender um den Trommel- 
umfang herumlaufender Streifen der Trommeloberfläche wird eine Spur genannt. Jede Spur 
enthält Platz zur Aufnahme von 4 Zahlen zu je 32 Dualziffern, einer zusätzlichen Vorzeichenziffer 
und 5 Leerziffern. Letztere werden als Pufferraum zwischen den Speichern bzw. zum Erkennen 
einer Rechenbereichsüberschreitung benötigt. Beim Ablesen der Zahlen aus den statischen 
Speichern wie in allen internen Verbindungsleitungen der Rechenmaschine werden die Zahlen 
dargestellt durch Impulszüge, bei welchen die Einsen durch Impulse und die Nullen durch Fehlen 
von Impulsen wiedergegeben werden. Der niedrigste Dualziffernwert der Position 2-2 erscheint 
zuerst. Es folgen die Positionen 2728 bis 2+? und schließlich in der Position 2+5 die Vorzeichen- 
zilfer. Einzelheiten über den Speicher sind im Abschnitt ‚Speichertrommel‘“ zu finden. 

Der Akkumulator ist das eigentliche Rechenelement der Maschine. Er erhält ein 
Addierwerk, welches ermöglicht, daß eine dem Akkumulator zugeführte Zahl zum vorherigen 
Akkumulatorinhalt hinzuaddiert bzw. -subtrahiert wird. Da bei der Multiplikation Zahlen von 
doppelter Länge anfallen, ist der Akkumulator zur Aufnahme einer 72stelligen Dualzahl befähigt. 
Der Akkumulator enthält — wie auch die drei Register Verteiler, Md und Mr — je einen dy- 
namischen Speicher. Ein dynamischer Speicher”) enthält (vgl. Bild 2) einen Schreibmagnet S, 
der die zu speichernde Zahl in eine Spur der Magnettrommel eingibt. Um einen gewissen Winkel 
versetzt — beim Akkumulator eine halbe Trommelumdrehung — befindet sich auf derselben 
Trommelspur ein Lesemagnet Z, der nach einer in die- 


sem Fall zwei Zahlenlängen entsprechenden Zeit den Schreibm, 
eingeschriebenen Dualzahlenzug wieder abliest. Der E 
abgelesene Impulszug wird verstärkt, wobei unvermeid- DO Löschm. 
liche kleine Formänderungen der Impulsgestalt auto- 

matisch regeneriert werden, und wird schließlich dem an 
Schreibkopf wieder zugeführt. Auf dem Rückweg vom 12% 


Lese- zum Schreibkopf werden die auf die Trommel 
aufgeschriebenen Impulse dieser Spur durch das Gleich- Bild 2. Dynamischer Speicher 

feld eines Löschmagneten gelöscht. Der eingegebene 

Impulszug läuft damit in dem aus Schreibkopf — Magnetschicht — Lesekopf und Verstärker 
gebildeten Kreis bis zu seiner willkürlichen Löschung durch Unterbrechen dieses Kreises un- 
verändert um. Beim Akkumulator ist vor dem Schreibkopf ein Addiergerät eingefügt. Über den 
zweiten Eingang des Addiergerätes kann zu der im Akkumulator umlaufenden Zahl eine andere 
vom statischen Speicher her als Ziffernzug einlaufende Zahl addiert bzw. subtrahiert werden. 
Die Subtraktion wird entsprechend durch Addition der Komplementzahl vorgenommen. 

Der Verteiler: Bei der Addition vom statischen Speicher zum Akkumulator müssen 
gleichwertige Ziffernpositionen der beiden zu addierenden Impulszüge zu gleicher Zeit an den 
beiden Eingängen des Addiergerätes zusammentreffen. Da der Ziffernzug im Akkumulator 
mit der Periode von zwei Zahlenlängen umläuft, die 4 in einer Spur des statischen Speichers 
untergebrachten Zahlen aber nur eine Dauer von einer Zahlenlänge haben, ergibt sich die Not- 
wendigkeit, aus dem 1. und 3. Quadranten des statischen Speichers entnommene Ziffernzüge 
vor Addition zum Akkumulator um eine volle Zahlenlänge zu verzögern. Diesem Zweck dient 
der Verteiler (Bild 2), und zwar in ihm der von L1 nach 8 geschlossene Kreis von der Umlaufs- 
dauer einer Zahlenlänge. Auch bei der Rückführung einer Zahl vom Akkumulator zum Zahlen- 
speicher wie beieiner Überführung vom Zahlenspeicher zum Md-bzw. Mr-Register ergibt sich 
die gleiche Notwendigkeit der Verschiebung um eine Zahlenlänge. Daher werden alle vom sta- 
tischen Speicher entnommenen oder in ihn hineingeschriebenen Zahlen zunächst durch den 
Verteiler gelenkt. Für den Befehl ‚‚zyliksche Vertauschung‘ wird die notwendige Verschiebung 
des Inhalts einer gesamten Spur um eine Zahlenlänge ebenfalls mittels des Verteilers vorge- 
nommen. 

Das Md- und Mr-Register werden zur Multiplikation, Division und zum Quadrat- 
wurzelziehen benötigt. Beide sind dynamische Register von 2 Zahlenlängen Umlaufsdauer. 
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Jedoch kann die Umlaufsdauer durch eine einfache elektronische Schaltung auf Befehl hin in 
Mr um eine Ziffernlänge verlängert und in Md um eine Ziffernlänge verkürzt werden. Das be- 
deutet, daß bei Betätigung der elektronischen Schaltung bei jedem Umlauf der Inhalt von Mr 
bzw. Md mit 2*+! bzw. 2! multipliziert wird. Zur Multiplikation nimmt Mr den Multiplikator 
und Md den Multiplikanden auf. Die Multiplikation geschieht dann auf die gewöhnliche Weise 
durch Addition der Partialprodukte zum Akkumulatorinhalt bei Verschiebung des Multipli- 
kanden um jeweils eine Ziffernlänge pro Addition. 

Die Division wird nach dem von v. Neumann?) beschriebenen Verfahren durch- 
geführt, welches für maschinelle Division im Dualsystem besonders geeignet ist. Der Dividend 
wird dabei zunächst in den Akkumulator und der Divisor nach Md verbracht. Der errechnete 
Quotient wird Ziffer für Ziffer in Mr aufgebaut und am Schluß der Operation nach Löschung 
des Restes des Dividenden zum Akkumulator überführt. 

Das zum Wurzelziehen verwendete Verfahren ist von Zuse°) angegeben worden. 
Der Radikand steht dabei zunächst im Akkumulator. Die Wurzel selbst wird aus technischen 
Gründen in Md und Mr gleichzeitig aufgebaut und ersetzt zum Schluß im Akkumulator den 
Radikanden. s 

Ein- und Ausgabewerk stellen die Verbindung zwischen Rechner und Maschine 
dar. Das Eingabewerk dient der Eingabe von Zahlen und Kommandos von der Schreibmaschinen- 
tastatur oder den Befehlsbändern her, während das Ausgabewerk lediglich Resultate druckt bzw. 
ins Lochband locht. Bei Zahleneingabe und -Ausgabe wird die notwendige Umwandlung vom 
Dezimal- ins Dualsystem und umgekehrt in dem langen Zyklus des Verteilers durchgeführt 
(Weg L2 nach $ Bild 2). Näheres hierüber findet sich im Abschnitt Ein- und Ausgabe. 


Das Befehlswerk bildet das technische Bindeglied, welches die vom Eingabewerk 
her eingegebenen Kommandos zur Ausführung bringt. Es enthält die von den Schreibmaschinen- 
tasten bzw. Lochbandabtastern her erregten Relais und den Schaltwellengenerator, der die zum 
jeweiligen Befehl zugehörigen Schaltwellen herstellt, welche zu den richtigen Zeiten die zu dem 
betreffenden Befehl gehörenden Verbindungswege zwischen den Recheneinheiten öffnen bzw. die 
Recheneinheiten steuern. Die Relais legen also nur den auszuführenden Befehl fest, die eigent- 
liche Operation wird rein elektronisch gesteuert. 

Das Synchronisationswerk ist als Hilfseinheit für die statischen und dyna- 
mischen Speicher, wie für das Befehlswerk anzusprechen. Es besteht aus der Uhr und einigen 
Zählwerken. Die Uhr ist ein fest mit der Speichertrommel verbundenes eisernes Zahnrad mit 
144 Zähnen, welches in gleichen Abständen pro Trommelumdrehung 144 Impulse abgibt. In 
den Speichern werden die Uhrimpulse verwendet, um die an den dortigen Lesemagneten an- 
fallenden gestaltmäßig und zeitlich schwach verformten Ziffernimpulse zu scharfen Impulsen 
zu regenerieren, die mit den zugehörigen Uhrimpulsen koinzidieren. Die Zählwerke des Syn- 
chronisationswerkes dienen dem Zweck, die 144 Impulse pro Trommelumdrehung und zusätzlich 
die seit Beginn eines Befehles abgelaufene Anzahl von Trommelumdrehungen abzuzählen, so daß 
es möglich ist, die Schaltwellen des Befehlswerkes bei Abstimmung auf die jeweilige Trommel- 
stellung zu den richtigen Zeiten beginnen und enden zu lassen. 


Befehlsplan 


Insgesamt werden 31 Tasten der Schreibmaschine und die Zwischenraumtaste zur 
Steuerung der Maschine verwendet. Jeder Taste entspricht eine bestimmte Lochkombination 
auf dem Befehlsband. Die Tasten sind bezeichnet mit den 10 Ziffern 0, ---, 9, den beiden Vor- 
zeichen (die mit 9 bzw. n bezeichnet sind, um Verwechslungen mit dem Additions- bzw. Sub- 
traktionszeichen zu vermeiden), den Buchstaben a, b, c, d, und den Zeichen + — mx x:DL 
]1/h>BFI. Ein Befehl besteht im allgemeinen aus mehreren Tastenanschlägen. Die in Frage 
kommenden Kombinationen von Tastenanschlägen sind in dem Befehlsplan (Tabelle 1) an- 
gegeben. ‚Einige andere nicht aufgeführte Kombinationen sind auch noch erlaubt, bieten aber 
keine wesentlichen neuen Möglichkeiten. 


Eine Zahl Z (|Z| < 8) wird in die Maschine eingegeben,.indem zunächst die Stelle vor dem 
Komma (auch wenn sie gleich Null ist), dann höchstens 9 weitere Stellen durch Anschlag der 
entsprechenden Zifferntasten und zum Schluß das Vorzeichen (p oder n) einzeln eingegeben 
werden. Die Zahl wird dabei ziffernweise im Verteiler als Dualzahl aufgebaut. Wenn weniger 
als 10 Ziffern eingegeben werden, werden bei Eingabe des Vorzeichens die fehlenden Ziffern 
als Nullen automatisch ergänzt. 


°) Preliminary Discussions of the logical Design of an Electronie Computing Instrument by A. W. Burks 
Goldstine,v. Neumann, Institute for Advanced Study Princeton 1947. 

°) Rutishauser, Speiser, Stiefel, Programmgesteuerte digitale Rechengeräte. Z.angew. 
Math. Phys. 1 (1950), S. 356. 
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Tabelle 1, Befehlsplan für die bandgesteuerte Rechenmaschine @l. 


Angegeben ist der Inhalt der Register A (Akkumulator), M (Multiplikandenregister) und z (Zahlspeicher) 

ausgedrückt durch den Inhalt vor Ausführung des Befehls (in <)). — = keine Änderung; 0 = Löschung. 

Speicherplätze z= 0,..., 9 und (zyklische Speicher) «0,...,a3; b0,..., b3; c0,..., e3; d0,..., d3. 

Z bedeutet eine höchstens 10-stellige Dezimalzahl, die ziffernweise eingegeben wird, beginnend mit der Stelle 
vor dem Komma (<7) und endend mit der Eingabe des Vorzeichens (p oder n). 


Operationen mit eingegebenen Zahlen Operationen mit gespeicherten Zahlen 
Nr. Befehl A M 2 Nr. Befehl A M 2 
aA CAS zZ a A 1b 2 A + tn 
23 De SAD 2 — —_ 2b —z (A> — <> — — 
dal Zm — Z _ 3bl mz — RM — 
4a ZIR <A + Z.<My> 0 — 4b x2 <Ay + (<)-<M>) 0 —_ 
5a Zx <A — Z»<MS 0) _ 5b 2 (Ay — (2) -<MS). 0 — 
6a Z: ANZ 0 — 6b 2 LAY: \<2) 0 _ 

Drucke die Zahl Z Drucke die Zahl 2) 

9a ZL _ — _ 9b Lz _ —_ u 

Drucke und loche die Zahl Z Drucke und loche die Zahl <z) 

Andere Operationen 

108] 0 LAS == 13 & zyklische Vertauschung <a0y—al 
11 Y VcA> En = b, <, d entsprechend <aly—>a2 
12 h (43:2 Yu sV | | <a2y—a3 
a <a3ay a0 


14 Ba Start Lochband «a 
Bb, Be, Bd entsprechend 
15 F _Lochband halt 
16 I _ Überlochen einer Fehllochung 
und Löschung eines Vorbefehls 


1) hiernach soll Befehl 4) oder 5) folgen. 


Mit der jetzt im Verteiler befindlichen Zahl Z kann nun direkt gerechnet werden. Nach 
dem Anschlagen der + oder — Taste (Befehl 1a bzw. 2a) wird Z zum Inhalt des Akkumulators 
(der im Befehlsplan durch A angedeutet ist) hinzugezählt bzw. von ihm subtrahiert. Durch 
Anschlag der Taste m (Befehl 3a) wird nach der Eingabe einer Zahl diese ins Multiplikanden- 
register gebracht und steht dort zur Multiplikation bereit. Die Multiplikation einer eingegebenen 
Zahl mit einer Zahl, die im Md-Register steht, wird durch den Befehl 4a oder 5a bewirkt; beim 
Befehl 4a wird das Produkt mit 64stelliger Genauigkeit zum Akkumulatorinhalt addiert, 
beim Befehl 5a dagegen subtrahiert. Division des Akkumulatorinhaltes durch eine eingegebene 
Zahl wird dadurch bewirkt (Befehl 6a), daß man nach der Eingabe der Zahl die Taste : anschlägt. 

Man kann die eingegebene Zahl aus dem Verteiler in einen beliebigen der 26 Speicher 
schaffen (Befehl 7a), indem man nach der Eingabe der Zahl die Taste > und anschließend die 
gewünschte Speichernummer anschlägt. Soll die Zahl in einen der ‚gewöhnlichen‘ Speicher 
gebracht werden, schlägt man die betreffende der 10 Zifferntasten an; soll sie in einen zyklischen 
Speicher kommen, wählt man zunächst durch Anschlag eines der Buchstaben a, ---, d die be- 
treffende Spur und dann durch eine der Zifferntasten 0, ---,3 den entsprechenden Sektor. 
Das Symbol z im Befehlsplan bedeutet also entweder eine der Ziffern 0, ---, 9 oder eine dieser 
16 Kombinationen. Bei der Ausführung dieses Befehls wird die zuvor im Speicher z befindliche 
Zahl gelöscht, die übrigen Register der Maschine bleiben unberührt. 

Die Befehle 8a) und 9a) besitzen nur geringes praktisches Interesse. Durch Befehl 8a 
wird die Zahl wieder in eine Dezimalzahl zurückverwandelt und ausgedruckt. Sie ist damit 
aus der Maschine verschwunden und steht zu weiteren Rechnungen nicht mehr zur Verfügung. 
Bei Befehl 9a wird sie außerdem ins Ergebnisband gelocht. 

Gemeinsam ist allen Befehlen, die mit gerade eingegebenen und dann im Verteiler stehenden 
Zahlen sofort weiter operieren, daß auf die Eingabe unmittelbar das Operationssymbol folgt. 
Der Befehl] wird dann sofort ausgeführt und dabei der Verteiler gelöscht. 

Wird dagegen eines der in den Befehlen 1a:---9a auftretenden Operationssymbole ohne 
vorherige Eingabe einer Zahl (also bei leerem Verteiler) gegeben, so wirkt es nur als Vorbefehl 
und die betreffende Operation wird erst ausgeführt, wenn anschließend eine Speichernummer 
(0, +, 9 bzw. a0, --, d3) angeschlagen wird. Die Befehle (1b, --, 9b) werden dann ganz analog 
mit der im Speicher z befindlichen Zahl ausgeführt; es sind nur folgende Unterschiede zu be- 
merken: Befehl 7b schafft bei leerem Verteiler die im Akkumulator befindliche Zahl nach Ab- 
rundung auf 32 Dualstellen in den Speicher 2 und löscht den Akkumulator dabei, während die 


- Befehle 8b und 9b die im Speicher z befindliche Zahl ausdrucken (bzw. zusätzlich lochen). 
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Die Befehle 10, 11, 12 wirken direkt auf den Akkumulatorinhalt, bedürfen daher nur eines 
Tastenanschlages und werden sofort ausgeführt. 

Befehl 10 dient mehr der Bequemlichkeit. Falls der Multiplikand bereits im Akkumulator 
steht, erlaubt er die Überführung des Multiplikanden vom Akkumulator zum Multiplikanden- 
register. Dies ist z. B. bei der Bildung mehrfacher Produkte nützlich. Multiplikation der Inhalte 
der drei Speicher 0, 1, 2 wird z. B. mittels der Tastenanschläge m0x 1] x 2 durchgeführt. 


Da nach den Befehlen 3a, 3b und 10 das Multiplikandenregister besetzt ist und das Vor- 
zeichen des Multiplikanden bereits im Akkumulator abgesetzt wurde, um dort im Verein mit 
den Vorzeichen von Multiplikator und Befehl 4 oder 5 die tatsächliche Ausführung der Mul- 
tiplikation zu steuern, darf zwischen Besetzen des Multiplikandenregisters und Ausführung der 
Multiplikation keine Addition, Subtraktion, Division und Radizierung liegen. Letztere beiden 
verbieten sich ohnehin, da sie das Multiplikandenregister benötigen. D.h. praktisch, daß die 
Multiplikation möglichst hintereinander weg ausgeführt werden soll. 


Befehl 11 zieht die Quadratwurzel aus dem Akkumulatorinhalt und überführt das Resultat 
in den Akkumulator. : 

Befehl 12 halbiert den Inhalt des Akkumulators. 

Befehl 13 führt die bereits besprochene zylklische Vertauschung der 4 in einer zyklischen 
Speichergruppe befindlichen Speicherinhalte durch. Er ist einzeln anwendbar auf die 4 zy- 
klischen Speichergruppen a bed. 

Die Befehle 14, 15, 16 dienen der Steuerung der Maschine von einem der 4 vorgesehenen € 
Lochstreifen aus. 


Befehl 14 erlaubt das gewünschte Lochband in Betrieb zu setzen oder von einem Band 
auf ein anderes überzugehen. Wenn der Befehl 14 von einem Lochband aus gegeben wird, so 
läuft das neue Lochband an und das alte wird automatisch stillgesetzt. Auf diese Weise kann 
man eine willkürliche Reihenfolge, in welcher die 4 Lochbänder abwechselnd die Komman- 
dierung übernehmen sollen, auf den Lochbändern verschlüsseln. 


Befehl 15 bringt das Lochband zum Stillstand, ohne ein neues in Betrieb zu setzen. Dieser 
Befehl wird durch eine Lochbandzeile gegeben, die in allen 5 Positionen ungelocht ist. Der 
Anfang eines Lochbandes ist normalerweise nicht gelocht. Um bei Inbetriebsetzung des Loch- 
bandes das Band über diesen ungelochten Teil hinweglaufen lassen zu können, ist der Stopp- 
befehl 15 technisch so durchgeführt, daß er nur wirksam wird, wenn nach dem Startbefehl 14 
von dem betreffenden Band vorher mindestens eine mit Löchern versehene Zeile abgetastet ist. 


Bei Eingabe von Zahlen von einem Lochband aus ist es meist erwünscht, nach Ablesen 
einer einzelnen Zahl automatisch auf ein anderes Band überzugehen, welches kommandiert, 
wie die abgelesene Zahl im weiteren Rechengang verwendet werden soll. Das Lochband, welches 
die einzugebenden Zahlen enthält, wird im allgemeinen von einer vorhergehenden Rechnung 
aus automatisch gelocht worden sein. Dieses hat zur Folge, daß das Zahlenlochband von’vorn- 
herein zwischen den einzelnen Zahlen den erwünschten Befehl 14 bzw. 15 nicht enthält. Am 
Ende jeder von der Rechenmaschine her ausgelochten Zah] befindet sich aber ohnehin das Vor- 
zeichensymbol p bzw. n. Zum Bandwechsel oder Bandstop nach Eingabe einer Zahl vom ‚di 
Zahlenband her sind daher mit jedem der Bandabtaster, von Hand aus setzbare Schalter ver- c 
bunden, welche bewirken, daß nach Abtasten eines in dem betreffenden Band enthaltenen 
Vorzeichensymbols automatisch je nach Schalterstellung auf ein gewünschtes anderes’ Band 
übergegangen wird oder Stopbefehl 15 ausgeführt wird oder die Bandabtastung unverändert 
weiterläuft. 


Befehl 16 wird verwendet, wenn beim Lochen des Bandes ein Irrtum unterlaufen ist. 
In diesem Fall werden nachträglich alle 5 Lochpositionen der falsch gelochten Zeile mit Löchern 
versehen. Diese Lochverteilung — sie entspricht Befehl 16 — löst in der Rechenmaschine keine 
Wirkung aus, die Maschine geht anschließend zum Abtasten der nächsten Bandzeile über. 


B 

Fast alle im Befehlsplan angegebenen Rechenoperationen benötigen eine Ausführungsdauer Re: 
von weniger als Y/, Sekunde. Bei Steuerung von einem Lochband aus beträgt die Abtast- 3 
geschwindigkeit 7 Zeichen pro Sekunde, so daß bei all diesen Operationen die Rechengeschwin- 4 
digkeit durch die Abtastgeschwindigkeit des Lochbandes begrenzt ist. Eine Ausnaıhmeblden 
lediglich die Befehle 4, 5, 6, 11. Bei Befehl 4, 5, 6 muß zur Ausführung der Multiplikation bzw. Be: 
Division nach dem Anschlagen der Zifferntaste z zusätzlich knapp /, Sekunde und nach Befehl 11 
zusätzlich knapp ?/, Sekunde bis zur nächsten Befehlsgabe gewartet werden. In diesem Fall 
wird das Lochband für die Dauer der Befehlsausführung angehalten. Die volle Multiplikation 
zweier Speicherinhalte — beispielsweise für die Befehlsfolge m0 x 1 — benötigt etwa %/, Sekunde 
+!/, Sekunde = 0,9 Sekunden. A: 


e) 
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Ein-und Ausgabe 

Bei der Zahleneingabe muß aus der Folge der vom Eingabewerk her eingeschriebenen 
Dezimalziffern die entsprechende Dualzahl aufgebaut werden. Dies geschieht im Verteiler, der 
hierzu auf den langen Zyklus von 2 Zahlenlängen geschaltet ist. Folgendes Verfahren wird an- 
gewendet: In 9 zusätzlichen statischen Speichern, die sich sämtlich im 1. bzw. 3. Quadranten 
der Speichertrommel befinden, sind als Umrechnungszahlen die Dualdarstellungen derf Zahlen 
1x 10710 pis 9 x 10710 auf 1x 2-% genau in Dualform ständig für die gesamte Lebensdauer der 
Rechenmaschine gespeichert. Durch Aufruf von einer Zifferntaste her wird die zugehörige 
Umrechnungszahl in den Verteiler eingegeben und der Verteilerinhalt anschließend dual mit 10 
multipliziert. Nach 10maliger Wiederholung dieses Vorganges ist der Aufbau der gewünschten 
Dualzahl beendet. Zur Durchführung der Dualmultiplikation mit 10 enthält der Verteiler eine 
einfache Zusatzeinrichtung, um entsprechend des Dualbildes von 10 (d.i. ZO LO) den Verteiler- 
inhalt mit 2! bzw. 2? zu multiplizieren (das geschieht durch eine Verschiebung der Zahl um 1 
bzw. 3 Stellen), und ein Addiergerät, um die beiden auf diese Weise gebildeten Partialprodukte 
zu addieren (Bild 2). Das Vorzeichensymbol wird nach der letzten Ziffer der Zahl eingegeben 
und löst einen Mechanismus aus, der einerseits zum Anhängen nicht bedeutsamer Nullen evtl. 
noch ausstehende Zehnermultiplikationen ausführt und andererseits bei negativem Vorzeichen 
einen zusätzlichen Impuls in die Vorzeichenstelle des Verteilers eingibt. 

Bei der Ausgabe erfolgt die Rückumwandlung in eine Dezimalzahl ebenfalls im Ver- 
teiler und zwar auf die übliche Weise durch Übertragung der unmittelbar vor! dem’ Komma 
stehenden aus 4 Dualziffern bestehenden Gruppe in die entsprechende Dezimalziffer, darauf- 
folgende Löschung der vor dem Komma stehenden Dualziffern und anschließende Multiplikation 
der restlichen Dualzahl mit 10. Die ausgegebenen Dezimalziffern können gedruckt oder auch 
zusätzlich in ein Band gelocht werden. Es ist Vorsorge getroffen, daß das Vorzeichensymbol 
am Ende der Zahl ausgeschrieben wird, so daß die Darstellung der Dezimalzahlen bei Ein- und 
Ausgabe identisch ist. Das angewendete Ausgabeverfahren arbeitet nur richtig, wenn die erste 
ausgegebene Ziffer die unmittelbar vor dem Komma stehende Dezimale darstellt. Dies ist der 
wesentliche Grund, weshalb der Rechenbereich der Maschine auf das Gebiet |Z] < 8eingeschränkt 
wurde. 


Abrundung 


Nach einer Multiplikation ist das Produkt im Akkumulator auf 58 gültige Dualstellen 
hinter dem Komma errechnet. Vor Ausgabe in einen statischen Speicher muß eine Abrundung 
auf 29 Dualstellen hinter dem Komma stattfinden. Dieses wird bewirkt durch Addition einer 1 
in die Position 273° und nachträgliches Fortlassen des bei der Position 2-30 beginnenden Restes 
der Zahl. Dieses Verfahren ist richtig, wenn der Erwartungswert des Restes 1 x273° beträgt, 
was nach Multiplikation der Fall ist. Nach Anwendung des Halbierungsbefehles auf den 
Akkumulatorinhalt gilt jedoch obige Voraussetzung über den Erwartungswert des Restes nicht, 
falls der Akkumulator vor der Halbierung bereits eine auf 29 Stellen hinter dem Komma ab- 
gerundete Zahl enthielt. In diesem Fall weiß man, daß nach der Halbierung von Position 27°! 
beginnend der Rest lauter Nullen enthält. Man darf daher die Zahl jetzt nieht nach dem nor- 
malen Verfahren abrunden. Ein systematischer Abrundungsfehler wird in der Rechenmaschine 
auch in diesem Fall vermieden, indem die Maschine vor der Abrundung prüft, ob der Rest 
genau 0,5 x 22° beträgt. Ist dies der Fall, so wird zur Abrundung die Position 2=?° unabhängig 
von ihrem vorherigen Wert zu einer 1 gemacht und der Rest fortgelassen. Dieses Verfahren 
bewirkt, daß der Rest von genau 0,5 x 272° mit gleicher Wahrscheinlichkeit nach oben oder 
unten abgerundet wird. Auch nach mehrfacher Anwendung des Halbierungsbefehles vermeidet 
dieses erweiterte Abrundungsverfahren einen systematischen Fehler. 


Die Speichertrommel 

Die Speichertrommel (Bild 3) besteht aus einem vollen Bondurzylinder von L=17cm 
Länge und D=8,3cm Durchmesser und wird mittels direkter Kupplung mit der Welle eines 
0,3kW Drehstrommotores mit 50 U/sek um die eingezeichnete Achse gedreht. Die Trommel- 
achse liegt horizontal. Da der Abstand zwischen der Trommeloberfläche und den Lesemagneten 
nur Y,,mm beträgt, muß die Trommel sauber gearbeitet und ausgewuchtet sein. Die magne- 
tisierbare Oberfläche ist hergestellt durch Aufspritzen einer zur Magnetophonbandherstellung 
entwickelten Magnetit-Dispersion (y Fe, O,), die von der Badischen Anilin- und Sodafabrik 
Ludwigshafen, unter der Typenbezeichnung LCH hergestellt wird. Unter allen von uns aus- 


_probierten Dispersionen eignete sich diese Type besonders gut, da einerseits bei ihrer Verwendung 
_ die an den Leseköpfen anfallenden Leseimpulse besonders groß waren (0,2 V) und andererseits 


die Spannungsamplitude des durch die Körnigkeit der Dispersion hervorgerufenen unvermeid- 


_ lichen ren um den Faktor 100 unter der Höhe der Leseimpulse liegt. Die Dicke der 
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aufgebrachten Magnetitschicht beträgt etwa 0,03 mm. Auf der Trommel beträgt = tangentiale 
Abstand zwischen 2 benachbarten Ziffernpositionen, die sogenannte Ziffernbreite " — Se = 2mm, 
die Breite einer Spur einschließlich Sicherheitsabstand zwischen benachbarten Spuren 10 mm. 
Die Trommel] rotiert in einem Eisenkäfig, der gleichzeitig zur Halterung der Magnetköpfe dient. 
Am linken Ende der Trommel befindet sich das Zahnrad für die Uhr. Es folgen die 3 Spuren 
für die dynamischen Register, Mr und Md, Akkumulator mit jeweils einem Schreib- und einem 
Lesemagneten. Die Repräsentation der Dualziffern O0 und 1 erfolgt in den Speichern durch 
tangentiales Magnetisieren senkrecht zur Achse bis zur Sättigung in der einen oder anderen 
Richtung (+Richtung oder — Richtung). In allen Verbindungsleitungen der Rechenmaschine 


Mr Md Akk VW. 


Uhr Mr Md. Akk. MV. Y RER DER; 


Bild 3. Die Magnettrommel 


werden Einsen durch Impulse und Nullen durch Fehlen der Impulse dargestellt. Um in den 
dynamischen Speichern Einsen in Nullen umschreiben zu können, werden die 4 dynamischen 
Spuren durch dauernd eingeschaltete Löschmagnete nach Passieren der Leseköpfe in det Minus- 
richtung magnetisiert (das bedeutet Löschung des Inhaltes) und die Einsen hinterher von den 
Schreibköpfen neu eingeschrieben. Es folgen 5 weitere Spuren, die in den 1. und 3. Quadranten 
die testen Umrechnungszahlen 1 x 10-10 his 9 x 10-10 und in den nullten und zweiten Qua- 
dranten die statischen Speicher 0 bis 9 enthalten. In den statischen Speichern erfolgt die 
magnetische Ziffernrepräsentation auf eine etwas andere Art dadurch, daß zur Darstellung der 
Ziffer 0 die erste Hälfte der Ziffernbreite in +Richtung und die zweite Hälfte in —Richtung 
magnetisiert wird. Bei der Ziffer 1 sind die Magnetisierungsrichtungen in beiden Hälften um- 
gekehrt. Dieses Verfahren bietet den Vorteil, daß auf eine vollständige Löschung des Speicher- 
inhaltes vor dem Einschreiben einer neuen Zahl verzichtet werden kann, da Einsen über Nullen 
und Nullen über Einsen unmittelbar übergeschrieben werden können. Ein zusätzlicher Lösch- 
magnet ist hier nicht notwendig. Stattdessen benötigt man ein einfaches elektronisches Gerät, 
welches die einlaufende Impulsdarstellung der Zahl in eine elektrische Welle umwandelt, die 
zur Darstellung von Nullen in der ersten Ziffernhälfte positive und in der zweiten Ziffernhälfte 
negative Spannungen hat. Die entgegengesetzte Wellendarstellung wird zur Darstellung der 
Einsen gebildet. Die beim Ablesen vom Speicher notwendige Rückwandlung in einen Impulszug 
ist bei beiden Arten magnetischer Ziffernrepräsentation gleich einfach. 


Die letzten 4 Spuren ganz rechts auf der Trommel nehmen die 4 zyklischen Gruppen a, 
b, c, d zu je 4 statischen Speichern auf. Die Ziffernrepräsentation ist die gleiche wie bei den 
übrigen statischen Speichern. 


Die Schreibköpfe, welche bei den statischen Speichern auch zum Lesen verwendet werden, 
ähneln den vom Magnetophon her bekannten Schreibmagneten. Sie erhalten einen lamellierten 
ringförmigen Eisenkern, der an der der Trommel zugewandten Seite einen 0,05 mm breiten 
Schlitz hat. Die dem Kern umgebende Magnetspule hat 250 Windungen. Die Eingangsimpedanz 
der Köpfe beträgt etwa 100 Ohm. Aus Abschirmungsgründen ist jeder Magnetkopf in einem 
Kasten aus Imm starkem Kupferblech eingebaut, aus welchem vorne mit dem Kasten ab- 
schließend die beiden Magnetpole heraussehen. 


Zuverlässigkeit 

Bei der großen Anzahl von Schaltelementen werden sich Ausfälle der Maschine nicht ganz 
vermeiden lassen. Soweit diese Ausfälle durch allmähliche Alterung von Widerständen oder 
Röhren auftreten, kann man diese gealterten Elemente vor dem Ausfall auffinden. Dies ge- 
schieht anläßlich von in regelmäßigen Zeitabständen durchgeführten Prüfungen, bei denen die 
beiden Hauptbetriebsspannungen der Maschine von -+280 V und —180 V um -+20%, geändert 
werden. Bei intakten Schaltelementen bleibt die Maschine rechenfähig, während Alterung 
von Elementen sich durch eine Verringerung des Änderungsbereiches bemerkbar macht. Die 
Wartung der Maschine ist relativ einfach. Da auf kompakte Bauweise verzichtet wurde, sind 
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sämtliche Schaltelemente zur Behebung von Störungen leicht zugänglich. Durch einen mit der 
Maschine vertrauten Techniker lassen sich Störungen im allgemeinen recht schnell beheben. 
Nachtrag während der Korrektur: 


Die Maschine ist im Sommer 1952 fertiggestellt und hat während der ersten Monate 
ihrer Erprobung so zufriedenstellend gearbeitet, daß sie zur praktischen Arbeitsaufnahme im 
September 1952 ihren endgültigen Aufstellungsort im Max-Planck- Institut für Physik, Göttingen, 
einnehmen konnte. 


Rechenbeispiele 
Zur Illustration der Anwendung der Rechenbefehle geben wir nun noch zwei Beispiele. 
1.) Es soll der Ausdruck 
(a —0,3 -4,5)b] :0,2 = c 
berechnet werden. Die Zahl a sei bereits im Speicher 3, die Zahl 5 im Speicher c2, während das 
Ergebnis c auf den Speicherplatz 4 gebracht werden soll. Hierfür schreiben wir kurz: 
dx): D’=X02%, c—4. 
Die Befehlsfolge, die entweder ins Lochband einzulochen, oder von Hand aus durch Anschlag 
der entsprechenden Tasten der Schreibmaschine zu geben wäre, lautet dann (mit den von der 
Maschine von allein geschriebenen Zwischenräumen) 
+3 089 m 459 x] X 2.029: —4. 
Das sind 20 Tastenanschläge bzw. Lochreihen. Die Maschine benötigt hierzu etwa 4 sec. 


2.) Es soll die lineare Differentialgleichung 


y" = fl® y); !(® y) = pe) 4 + g(®) 
schrittweise mit der Schrittweite A integriert werden. Hierzu benutzt man zweckmäßig folgendes 
Verfahren, das wir als das der ‚direkten Interpolation‘ bezeichnen wollen. Führt man für alle 
Punkte (d.h. für die Stellen & = x,, wobei &,+1ı = &, + h) die Bezeichnungen 


EN Zn EEE hr 3 
H-nra egal Bed wtn= ne, vw) 


ein, wobei %, die berechnete Näherung für yan der Stelle x, sei, so folgt nach diesem Verfahren 
der Wert von y an der Stelle v»—=n +1 aus den zuvor berechneten y=nundv=n—|1) 
nach der Formel 


[2 Yn — Yn—ı En In+iı =e 10 f, A 


sah Da 


In+ı = nr Yuzıt IMn+ı - 


Der pro Schritt erzeugte Fehler ist hierbei von der Ordnung 


ee 

240 da® 
Die Folge der 7, stehe auf Band 5; die Folge der 9, auf Band c. Die Besetzung der Speicher 
nach dem nten Schritt sei “ 


Fi Br <a0) VE ze <b0) 0,1 zn <0) 
Faı = (aly mı= (bl). 


Die dynamischen Register seien leer. Die Verschlüsselung eines Integrationsschrittes, der zu 
einem entsprechenden Zustand der Speicher mit (n+1) an der Stelle von n führt, sieht dann 
so aus: (In der zweiten Spalte ist der Inhalt des Accumulators nach Durchführung der jeweiligen 
Befehlsgruppe angegeben): 


Befehlsgruppe Inhalt des Acc Bemerkungen 


1.1p+ 1 Die Zahl 1 wird in den Acc. addiert. 
DABE 1 Die Zahl auf Band 5, die an der Reihe ist (hier ?„.+.1) wird 


in den Verteiler geschafft. Nach Ablesen des Vorzeichens 
von Pn+ı kehrt der Rechengang zum Band « zurück. 


Befehlsgruppe Inhalt des Acc Bemerkheen - EEE 
3:.—1 d P„+ı wird in den Speicher 1 gebracht und dort zur späteren E 


Bildung von f„,ı aufgehoben (s. 21.) 7 
4. —1 1 — Parı £ 
> 2 0 1 —Pn+ı wird in den Speicher 2 gebracht. Der Acc wird 
£ dabei automatisch gelöscht. (Vgl. 2.) Be. 
6. +a0 je u 
RO ef 0:1 Division durch 0,1, da eine einfache pn mit E10 | e 


(> 8!) nicht möglich ist. “ 
8 tal 10f, +fa—ı 4 
9, Be RE 9n+1 wird abgelesen und in den Verteiler gebracht. g+ 
10. —3 2 9n+1ı wird im Speicher 3 zur späteren Bildung von Tun auf- a 
gehoben (s. 20) | $ 
P 2 $3 
11. +3 10 tn + In+ı r 
12 ae, b 1 „ Yn—ı h 
3. +50 >. m Se e 
14. +50 „ +29, Die Bildung der eckigen Klammer ist damit abgeschlossen. €e 
195752 Yn+ı1 5.)!. Damit ist das neue % berechnet. Die folgenden Befehle Er 
dienen nur der Ausgabe dieses Ergebnisses und der Neu- 
u 3 ordnung der Speicherinhalte, damit der nächste Integrations- 
. ‚schritt nach demselben Schema durchgeführt werden kann. 
16. 5 Ya+ı Zyklische Vertauschung der Speicher der ‚Gruppe b, um 
Yn.ı auf den Platz bO zu bringen. 
17. — b0 0 
B. 18. Lb0 0 Yn+ı wird von der Schreibmaschine geschrieben und in er 
> Ergebnisband gebracht. 
19, a 0 Zyklische Vertauschung analog 16. : 
21. x. In+ı + (Pr+ı ' Yn+1) { ER 
e en = Jar a 
le 22. > a0 0 Damit ist: der Zyklus geschlossen. 27 


Insgesamt sind dies 54 Einzelbefehle, bzw. Lochreihen im Befehlsband a (= A cm Band- 
länge). Einschließlich der Ablesung der Bänder b und ce und der Ausgabe benötigt die Maschine 
etwa 12sec für einen solchen Integrationsschritt, für 50 Schritte also etwa 10 min. 


Als Beispiel wurde nach diesem Schema die Differentialgleichung y’ = yyaas —ıı mit der Re 
Maschine integriert, deren Lösung zu den gewählten Anfangsbedingungen (y = -Y&) bekannt ist, 


Br: Unsere Erfahrungen bei der Verschlüsselung einer großen Zahl verschiedener Probleme e Eu 
haben gezeigt, daß es möglich ist, in wenigen Tagen oder Wochen die Aufstellung der aa, on 
pläne für Probleme zu erlernen, wenn diese in einer Form vorliegen, daß sie direkt mit ı den üb- 
lichen Methoden der numerischen Analysis bearbeitet werden können. Es entfäll t hier ‚ja Er pen 
wesentliche Komplikation für die Befehlsgebung, die bei einer a ı Maschine 
bei der Aneinanderfügung verschiedener Be und allen mehr „logise en“ 


entstehen. 


2 


Eingegangen am 29. März 1952. 


_ tervall ©): - ©, die Momente M,, M,, M,,.. 
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KLEINE MITTEILUNGEN 


Zur Approximation einer Wahrscheinlich- 
Bei Yorteilung mit Hilfe ihrer niederen Mo- 
mente. 


Um eine Funktion unter Verwendung ihrer niederen 
Momente zu approximieren, bedient man sich zweck- 
mäßigerweise eines Systemsvon Näherungsfunktionen, 
die bereits irgendeine Ähnlichkeit mit der zu approxi- 
mierenden Funktion besitzen. Bei Annäherung einer 
Wahrscheinlichkeitsverteilung wird man z.B. daran 
denken, die Gaußverteilung in die Approximations- 
methode einzubauen, d.h. die Hermiteschen 
Polynome zu verwenden. Im Folgenden soll kurz ge- 
zeigt werden, daß man neben Forderungen über die 
Werte der niederen Momente auch noch andere For- 
derungen erfüllen kann. Man kann zusätzlich Forde- 
rungen über einzelne Funktionswerte erfüllen, und 
man kann außerdem erreichen, daß sich die ge- 
wünschte Näherungsfunktion in vorgegebenen Teil- 
bereichen möglichst wenig von Null unterscheidet. 

Erfüllung von Forderungen über 
die Werte der ersten m Momente. Es 
sei vorerst an einige aus dem Schrifttum bekannte 
Dinge erinnert. Eine Funktion w(x), die in einem In- 
. haben 
möge, soll in dem betreffenden Intervall durch die 
Näherungsfunktion 


W(2) = a, H(2)W, (2) + H(a) W;(e)) Aa 
+0, H(@) Wa(a) + a 


dermaßen angenähert werden, daß 


% 
fW(e)de = M,, 
@, 


T- 


% 
Kr Wis)aa=M,, ==. 
% 


Dabei sind die x; zunächst unbestimmte Koeffizienten, 
fernerist H(x)eine willkürlich wählbare Funktion, die 
zweckmäßigerweise irgendwelche Ähnlichkeiten mit 
der zu approximierenden Funktion hat, und außerdem 
sind die Funktionen W;(x) Polynome, die so einge- 
richtet werden, daß die Berechnung der Koeffizienten 
a;recht einfach vor sich geht. 


Analog der Herleitung der Legendreschen 
Polynome setzt man an: 
Wi; 
W,=z-+a, | 
w=ne4beltb, el 2) 
an | 


- und bestimmt die Größen a, b;, e;,  - : derart, daß 


F His) Walz) Wn(a) da = 0 für m=#n. 


Dies wird erreicht, wenn man a, b;, c;, * so bestimmt, 
daß 


FH) z Wı(z)da&=0. für m<n. 
% 


‘ Man gelangt, wenn die Momente von H(x) mit p,, Pı, 
Ps," bezeichnet werden, damit zu dem Gleichungs- 
system 


9 tam =, 

%+bmn +bm =; 

BB +bm +bmı=0; 

BB +4R+&aMmTt5Mm=0; 
mMtar tom +5 =0, 
Ps + CıPı + %aP3 + CP, = 0 


Mit 9, #0, S % == 0 und dem Nicht- Verschwinden 
aPı 

der übrigen in Frage kommenden Determinanten 

sind die Größen a, b;, e;, * - - bestimmbar. 

Unter Benutzung der Polynome W;(x) kommt man 
zu dem gewünschten Gleichungssystem, aus dem die 
Größen a; bequem erhältlich sind. Multipliziert man 
nämlich 


W=HW+&HW,+0,HW,+:-- 


mit W, bzw. W, bzw. W, usw. und integriert von ®, 
bis ©,, dann entsteht das Gleichungssystem 


%, 
&SHWide +, /HW,W,de 


a +8/HW,W,de+:'=M, 
&JSHW,W,de +0, /H Widx 
+%8/HW,W,de+---=M, +taM,, 
of HW,W,de +0,/H W,W,dx 
+0, /HW2d&e + -=M,+b,M,+b,M, 


%, 
in dem alle Integrale / H W„W,„ds für m ++ n ver- 


2 
schwinden. Die &; ergeben sich einfach zu 


M, M,+aM 
ee, 
SH W}dx /HW}dx 
en Rı 


._M+hM+b,M, 
RE WM, ... 


% 
SH Wida 
T, 


Wie bei Verwendung der Legendreschen Poly- 
nome berechnen sich die Koeffizienten unabhängig 
voneinander. Aus der Näherungsfunktion, die sich 
unter Zugrundelegung der ersten m Momente ergibt, 
gewinnt man die Näherungsfunktion, die sicb unter 
Zugrundelegung der ersten m + 1 Momente ergibt, 
durch Addieren eines dem (m + 1)-ten Moment ent- 
sprechenden Gliedes. Die den ersten m Momenten 
entsprechenden Glieder bleiben hierbei unbeeinflußt. 


Forderungen über einzelne Funk- 
tionswerte. Hat maneine Näherungsfunktion W 
unter Zugrundelegung der ersten m Momente gefunden 
und wünscht man, daß die Näherungsfunktion nicht 
nur diese m Momente aufweist, sondern auch durch 
den Punkt (7, Y%) geht, so braucht man nur das 
(m + 1)-te Moment als eine vorerst unbestimmte 
Konstante aufzufassen und diese Konstante dann der- 
art zu bestimmen, daß die Näherungsfunktion an der 
Stelle # den Wert y hat. Das ist gleichbedeutend da- 
mit, daß man xm-+1 so bestimmt, daß die Näherungs- 
funktion durch den Punkt (7, %) geht. Wie man leicht 
einsieht, wird diese Forderung erfüllt, wenn 


. (48). 


Mit einem Querstrich über einer Größe ist angedeutet, 
daß die betreffende Größe an der Stelle x zu nehmen 
ist. Man sieht, daß als Näherungsfunktion hier also 
auch Formel (1) mit den durch (2) gegebenen Poly- 
nomen W; zu verwendenist, wobei aber nur die Koeffi- 
zienten &g, &, &g, *'", &m mit Hilfe von (3) zu berechnen 
sind, während &m +1 nach (4a) zu berechnen ist. Wir 
schreiben den zu benutzenden Näherungsausdruck V 
folgendermaßen: 

V=W 4 KnPe HWm+ı . (4b), 
H Win+ıi 
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wobei W den durch (1), (2) und (3) gegebenen Nähe- 
rungsausdruck bedeutet, für den nur die Werte der 
ersten m Momente vorgeschrieben sind. 

Neben dem soeben behandelten Fall betrachten 
wir auch die Möglichkeit, daß die Näherungsfunk- 
tion V, die unter Zugrundelegung der ersten m Mo- 
mente gefunden ist und durch den Punkt (7, 7) geht, 
außerdem auch noch durch einen Punkt (7, Y) geht. 
Bei entsprechender Behandlung dieser Aufgabe fassen 
wir nunmehr auch das (m + 2)-te Moment und hier- 
mit &m+2 als zunächst unbestimmte Konstante auf 
und bestimmen &m-+2 derart, daß die gestellten For- 
derungen erfüllt werden. Esist unschwer zu erkennen, 
daß man das Gewünschte erreicht, wenn man zur 
Näherungsfunktion (4b) ein weiteres Glied addiert und 
somit den folgenden Näherungsausdruck benutzt: 
ER 

H Wm+1 


F - W)H Wm+ı — (FG — WHWm+1$ (5). 
H Wn+ı H(Wm+ı Wm+2 — Wm+1Wm+2 
-H (Wn+1 Wm+2 — Wm+2 Wm+1) 


Mit einem bzw. zwei Querstrichen über einer Größe 
ist angedeutet, daß die betreffende Größe an der 
Stelle # bzw. @ zu nehmen ist. 

Für den Fall, daß die Näherungsfunktion außer den 
vorstehend erfüllten Forderungen auch noch durch 


= 


einen Punkt (z, 7) oder noch weitere Punkte gehen 
soll, wären entsprechende Formeln abzuleiten. Im 
allgemeinen ist dies jedoch für die praktisch auf- 
tretenden Bedürfnisse nicht von Bedeutung. 

Forderungen über bestmögliches 
Verschwinden der Näherungsfunk- 
tion ineinem Teilbereich. Praktisch 
recht belangreichist die Frage, auf welche Weise man 
eine zusätzliche Forderung über bestmögliches Ver- 
schwinden der Näherungsfunktion in einem vorge- 
gebenen Teilbereich, der auch mehrteilig sein darf, zu- 
wege bringen kann. Eine solche Fragestellung taucht 
auf, sobald man für H (x) irgendeine sich über einen 
unendlichen Bereich erstreckende Funktion wählt, 
währenddessen die zu approximierende Wahrschein- 
lichkeitsverteilung nur in einem endlichen Bereich von 
Argumentwerten existiert. 

Wir nehmen uns zunächst den Fall vor, daß die 
Werte der ersten m Momente vorgegeben seien und 
daneben gefordert wird, daß sich die Näherungsfunk- 
tion in einem Bereich z, --- x, möglichst wenig von 
Null unterscheidet, wobei „möglichst wenig unter- 
scheiden‘ im Sinne der Methode der kleinsten Qua- 
drate verstanden sein möge. Das Gewünschte wird 
hier offensichtlich erreicht, sofern wir unter Zu- 
grundelegung der ersten m Momente vorerst eine 
Näherungsfunktion W gemäß (1), (2) und (3) ausfindig 
machen und dann weitere Glieder gemäß Formel (1) 
addieren, in denen zwar die durch (2) gegebenen Poly- 
nome Wm+1, Wm+2, Wm+3, **: verwendet werden, 
aber die Koeffizienten &n-+1, X&m-+2,&m+3,'** derart 
bestimmt werden, daß 


© 
S(W +am+1H Wm+i 
©s + o&m+2H Wm+2 +: ')?dz = Min. 


Die Koeffizienten am +1, &m +2, &m-+3, * "sind aus nach- 
stehendem Gleichungssystem erhältlich: 


17 
&m+ı/HWin+ıda 
A + o&m+2/H? Wm+1Wm+2de + = 


— /WHWm+1da, 
(6a). 


> %m+1/H? Wm+1 Wm+2 da 


+0m+2/H?Wm+2de + :- 18 
—/WHWm+2dx 
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Die gesuchte Näherungsfunktion läßt sicb schreiben: 
U=-W +&m+ıH Wn+i +&m+2 HWm+2 + a (6b), 


wobei die Koeffizienten am+1, &m-+2,'* durch (6a) 
bestimmt sind. 


Als nächsten Fall fassen wir die Möglichkeit ins 
Auge, daß die Werte der ersten m Momente gegeben 
sind, daß ferner die Näherungsfunktion durch einen 
Punkt (zZ, Y) gehen soll und daß sich die Näherungs- 
funktion außerdem in einem Bereich 2; -- - x, möglichst 
wenig von Null unterscheiden möge, Um in diesem 
Falleine Lösung zu finden, sehen wir uns Formel (4b) 
nochmals an und überlegen uns, daß wir auch noch 
zahllose andere Formeln bilden könnten, die dasgleiche 
wie Formel (4b) leisten. Durch lineare Kombination 
solcher Ausdrücke vermögen wir leicht eine Formel zu 
bilden, die einerseits dasselbe wie Formel (4b) leistet 
und andererseits beliebig viele unbestimmte Kon- 
stante linear enthält. Eine so beschaffene Formel ist 


U=-W+5BB; 
ci 
7 W (Ta). 


== —= HAWm+i und ZPßi —=] 
Tı 


H Wın+i 


Über die ß,; verfügen wir nun derart, daß wir fordern, 
% 
es möge f U?dz= Min. sein, und gelangen zu dem 


Ts; . 
folgenden Gleichungssystem, aus dem die ß; berechen- 
bar sind: 


BFkBeS B,) d« 
u +Pß3S(Bz— Bı)(Bz — B)de+ = 
— /(W+ B,)(B,— B,)de, 
BaS (Ba — Bi) (Bz — B,)de (7b). 
+ Bf (Br BP de + = 
— S(W + B,) (Ba — B,) da 


. 


Bun. de, sap ui 


Ppıtßet+ß +" =1 : 


Da alle B; hier den Faktor y = enthalten, lassen 


sich in den vorstehenden Gleichungen entsprechende 
Kürzungen durch diesen Faktor vornehmen, 


Als weiteren Fall betrachten wir schließlich die 
Möglichkeit, daß die Werte der ersten m Momente vor- 
gegeben sind, daß ferner die Näherungsfunktion durch 
die Punkte (2,9%) und (7, y) gehen soll und daß sich 
die Näherungsfunktion außerdem in einem Bereich 
&; *"" &, möglichst wenig von Null unterscheiden möge. 
Um hier eine entsprechende Lösung wie im vorstehen- 
den Fall abzuleiten, ziehen wir Formel (5) heran und 
erkennen leicht, daß die folgende Formel dasselbe wie 
(5) leistet und dabei beliebig viele Konstante enthält: 


VS WER BEE (88) 
mit 
Be I En 
H Wnti Be 
+ _G-MEWnH—- G-WAWm _ 
H Wm+iH (Wm+i Am+i+1 — Wm+i Wm+i+1) 
HE (Wm+i Wm+i+1 — Wm+i Wm+i+1) 
und 


ZB; =1. 
€ 


% 
Wir bestimmen die ß; wieder derart,daßf U?de=Min. 


% 
und gelangen zu dem folgenden Gleichungssystem, 


aus dem die ß; berechenbar sind: 


n 


EN 
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BaS (Ba — Bi)? de 
= +ß3/(Bs —B)(B, —B)de+:-- 
—/(W + B)(Bs — B,)de , 
B2S (Bz — B) (B, — B,) de 
TRflB = Byraa tr =) 
— /(W + B,) (B, — B,) de 


l 


(8b). 


Pıt+B+Bs+ -=1 

Hinweis. Zum Schluß sei auf ein System von 
Funktionssystemen hingewiesen, dessen Verwendung 
bei Approximation von Wahrscheinlichkeitsvertei- 
lungen bisweilen nützlich ist. 

Wählen wir 4 (x) = »4e-= mit x, = 0 und ler 
so vermögen wir für jedes positive ganzzahlige A (sowie 
). = 0).die Polynome W7(x) der Formel (1) allgemein 
anzugeben mit 


Wr = ak — (k +) (1) a1 | 
ERLARFA— (br 


—(k+A) (+2 —1)(k+2— 2) (k) ars. 


+ Ye k+m!, % 

wobei R=0% 1, 2,3, **. Beider Bildung der zuge- 
hörigen Koeffizienten a7; beachten wir, daß hier 

fH Wide =kt.(k +9)! 

{) 
Für} =1z. B. sieht die Reihe der ersten Näherungs- 
glieder folgendermaßen aus: 
M, &% M,— 2M, 


ee ran a N) 
M,—-6M,+6M, 
VERSEETET 62494], 


wobei die Momente M7. selbstverständlich auf den 
Nullpunkt der x-Skala bezogen sind. 

Durch geeignete Wahl von A können wir bereits der 
Funktion H(x) irgendeine Ähnlichkeit mit der zu 
approximierenden Funktion aufprägen, indem wir 
ein erstes Hauptcharakteristikum berücksichtigen. 
Für H(x) = »4e-2 gilt im Bereich 0 :--- oo, daß das 
auf den Erwartungswert bezogene Streuungsquadrat 
0: —=i +1, die Schiefe — is: und der Exzeß 

„A+3 +1 
=.) REIER 
wachsendem A gegen die entsprechenden Werte der 
Gaußverteilung. 

Göttingen. 


Schiefe und Exzeß streben bei unbegrenzt 


Otfrid M. J. Mittmann. 


Numerische Integration 
impliziter Differentialgleichungen 


Implizite Differentialgleichungen treten in den An- 
wendungen nur verhältnismäßig selten auf und sind 
daher in der praktischen Mathematik kaum betrachtet 
worden. So wird z. B.in dem umfassenden Werk von 
Collatz[1]lediglich in einer kurzen Bemerkung im 
Anhang ein Verfahren zur Lösung der Gleichung 
F(x, y, y') = 0 angegeben, das überdies an gewisse ein- 
schränkende Vorsausetzungen gebunden ist. Im 
folgenden werden zunächst einige rechnerische Ver- 
einfachungen für das von Collatz vorgeschlagene 
Verfahren beschrieben und dann weitere Methoden 


angegeben, die von den erwähnten einschränkenden 


Voraussetzungen frei sind und i.a. auch weniger 
Rechenarbeit erfordern. 
1. Nach dem Vorschlag von Collatz wird die 


_ implizite Differentialgleichung 


Besyy0n re.) 


mit der Anfangsbedingung 


y(%) = Yo, y’(z)=yY (Filz, Yo, Yo) =(0) (2) 
durch die Substitution y’ =! in das explizite System 
da dy > 
d£ —=f(7,9y,t), 7 HT REN (9) 
mit 
F; 

EN 
S®, 9, 0) P,+ıR, (4) 


übergeführt, wobei die Lösungsfunktionen «(t), y(t) 
den Anfangsbedingungen 
z(y)=%, Yo) = % 
genügen. Das Verfahren ist nur durchführbar, wenn 
für die (als zweimal differenzierbar angenommene) 
Lösung y (x) von (1), (2) gilt i 
y’#+0, Fıty'Fy#d. 
Während man nun gewöhnlich zur Lösung eines 
Systems von zwei Differentialgleichungen nach einem 
der bekannten Differenzenverfahren zwei Differenzen- 
schemata benötigt, kann man für das spezielle System 
(3) auch mit einem Schema auskommen, Das kann 
man auf verschiedene Weisen erreichen, z. B. folgen- 
dermaßen: 
Setzt man zur Abkürzung f* =/f, so gilt mit der 
Schrittweite Ah 


v = mVifk HihVimf—, 


so daß man mit den üblichen Koeffizienten 


ADS, 


1 
uw +1)... w+i ! 5 
u(w +1) z er) du  .(5) 


U 
0 
für das Extrapolationsverfahren die Formeln 


je 
RN ED VA ee ee (0) 
i=0 


E 
Yx+1 — 9, = Rh ER avi: 
er 
r 
—= tg (&% +1 — 27) + RM? I ia; Vin fı-ı (7) 
i=1 


bekommt. Analog erhält man für das Interpolations- 
verfahren mit den Koeffizienten 
0 


;. RT RR 


du . (8) 


i! 
1 
das Gleichungssystem 


& 
dr 2, = RED RU TR a ee re LO) 
0 


{ 
Y41ı 7 9 = barı (&rrı — @e) + Ben ib, Vi1 fe 
| (10) 


das normalerweise durch Iteration gelöst werden 


kann, wobei man zweckmäßig die Näherungswerte 

1) +1 R y DjB, lo). 
a» YL) aus geschätzten Werten x)’, y; nach 
den Formeln 


#-+1) Ian () (v) 
2,41 241 = har er.ı » %+1> (+1) 


an) =12.. (1) 
+1) WIESN +1) (®) 
ir larhı — 241] » 
20,12...) 
berechnet. 


2. Denkt man sich in (1) die Lösung y (x) eingesetzt 
und differenziert man nach x, so ergibt sich nach 
Division durch F,’ die mit (1) äquivalente Differen- 
tialgleichung zweiter Ordnung 

IE Fy 2 y' Fy 
1 Fyr 


N nn Zu Sa a an Zr 2 EEE 
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mit der Anfangsbedingung (2). Die hierfür notwendige 
Voraussetzung Fy’ = 0 muß man sowieso machen, da- 
mit (1) eine eindeutige Lösung besitzt. Allgemeiner 
kann man auf diese Weise die implizite Differential- 
gleichung n-ter Ordnung 


EUER REM) —.V, erde) 
in die explizite Gleichung (n + 1)-ter Ordnung 


FytyFty" Fr... +ym Fe 
FR) (15) 
verwandeln, die nach jedem der bekannten Verfahren 
für Differentialgleichungen höherer Ordnung (vgl.[1]) 
gelöst werden kann. Noch allgemeiner kann man 
durch Differenzieren nach « ein beliebiges implizites 
System von Differentialgleichungen n-ter Ordnung 
überführen in ein explizites System von ebensoviel 
Gleichungen (n + 1)-ter Ordnung; man muß dabei 
nur voraussetzen, daßin den differenzierten Gleichun- 
gen die Koeffizienten der (n + 1)-ten Ableitungen 
eine nicht singuläre Matrix bilden, so daß man nach 
diesen Ableitungen auflösen kann. 
3. Man kann die Gleichung (1) aber auch inte- 
grieren, ohne sie überhaupt explizit zu machen. So 
wird in [2] gezeigt: Setzt man 


ynHl) = 


‚ ’ Fe, 9 y') 
,9y,Yy)=YV -— u ——..(16 
ee Fy' (%o, Yo, Yo) ee 
und bildet — ausgehend von einer weithin willkür- 


lichen Funktion y; (2) — die beiden Funktionenfolgen 


Yr-+1(2) =/f(e, y (2%) ’ y(2)) ’ 
20.12... ae 


so konvergiert die Folge y, (x) unter sehr allgemeinen, 
in der Praxis wohl stets erfüllten Voraussetzungen 
gegen die Lösung y (x) von (1), (2). Diese Methode ist 
eine Kombination des Verfahrens von Newton zur 
Berechnung von Nullstellen und des Verfahrens von 
Picard-Lindelöf zur Lösung expliziter Differential- 
gleichungen und läßt sich genau wie dieses (,,Ver- 
fahren der wiederholten Quadratur‘‘) numerisch ver- 
wirklichen, wenn man die in (17) auftretenden Inte- 
grale mittels numerischer Quadraturformeln berech- 
net. Hinsichtlich der Rechenökonomie ist diese 
Methode allerdings denselben Einwänden ausgesetzt 
wie das Verfahren der wiederholten Quadratur. Nume- 
risch zweekmäßiger ist es, das Newtonsche Null- 
stellen verfahren mit einem der Differenzenverfahren 
für Differentialgleichungen zu kombinieren: 


“ Hat man sich durch eine Anlaufsrechnung (etwa 
nach 2.) gewisse Anfangswerte y, = y(a%), y, = y(2,) 
beschafft, so berechnet man y, ,, mittels der Koeffi- 
zienten (5) aus 


er ; 
Yı = Y, hd zZ Viye 18) 
BE 
und Yır1 — etwa mittels des Newtonschen Ver- 
fahrens — als Nullstelle der Gleichung 
Flaryı» 941» M4)=0 - - (19). 


Während diese Methode offenbar dem Extrapolations- 
verfahren entspricht, führt die Übertragung des Inter- 
polationsverfahrens auf das Gleichungssystem 


T: 
Yrı yeah 2 BiVYyyrı, di wiein (8). . (20) 


F(a.11> Ya+ı> Y.41) = 0 er (21). 


Um daraus y,,], Y;.,.1 zu berechnen, beschafft man 


sich am besten durch Schätzung von Vr y, ,, und an- 
schließende Ergänzung des Differenzenschemas Nähe- 


rungswerte Pi vn G=0,1,...,r) und daraus 


durch Einsetzen in (20) eine Näherung en Durch 
Berechnung der Nullstelle von 


F(ay41» Yarı> Yarı) = 


_ erhält man ven und daraus y{_ nach (20). So 


k+1 

fortfahrend kann (20), (21) durch Iteration gelöst wer- 
den, wobei man in den weiteren Iterationsschritten 
(20) durch 


(+1) LAN ae ’$+1) ’) 
Y:41 — Y4ı = Ra (v Br Yerı) > 


la) 2 sen 22) 
ersetzt. Auch wird esi. a. zweckmäßig sein, für DEE 
nicht die genaue Nullstelle der Gleichung 


(ar, Yan» n)=0  -. (3) 


zu nehmen, sondern sich mit der ersten Newtonschen 
1 
7 P—l) . 
Verbesserung des Wertes y,” |, ' zu begnügen. 


Zu einer Fehlerabschätzung kann man wie folgt 
gelangen. Bezeichnet man mit y die exakte Lösung 
von (1), (2), mit y wie bisher die Näherung und setzt 


,=Yy 4» u u 7 k=0,1,2,..., 
so folgt aus den Gleichungen 
0 =F(a,, Y.» 2) =F(a,, Y,> y;) . . (24) 
und dem Mittelwertsatz der Differentialrechnung die 


Beziehung 
* 
&% er ee Ek > 
Fy kr 


wobei der Stern andeuten soll, daß als Argument ein 
Zwischenwert der Argumente in (24) zu nehmen ist. 
Schreibt man nun z. B. die Formel (20) des Inter- 
polationsverfahrens mit 


r fi : 
RL () bi 


ı=0o 


NTSEESE (OB) 


in der Form 
Ya+ı 79T n EB: Yarızaa 
und analog für die exakte Lösung 
Yrı-d, mh Ep Yerıot % 


mit dem ‚ Quadraturfehler“ g7, so folgt durch Sub- 
traktion 


r ‚ 
SE+1 2ER BR E41 F IE 


und wegen (25) 


&41-g 7 IE: 
RE Nie 


Diese Fehlergleichung kann genau wie die entspre- 
chende Gleichung (2,3) in [3] weiterbehandelt werden 
und führt zu denselben Fehlerabschätzungen wie in 
[3], wenn man die dortige Lipschitzkonstante K jetzt 
durch 


r * 
u dee), 


e=0 \Fyr/a+1—e 


Ei 


y' 


K= max 


a 


a 


ersetzt, wobeisich das Maximum auf einen (konvexen) 


Bereich im (x, y, y’)-Raum bezieht, der die angenä- 
herte und die exakte Lösung enthält. — Ganz ent- 
sprechend findet man übrigens durch Anwendung des 
Mittelwertsatzes auf die för 

Konvergenz des Iterationsprozesses (22), (23) die Be- 
dingung Kha, <1 hinreichend ist, wenn man das 


Maximum in (27) in einem les, die Näherungs- = 
werte enthaltenden Bereich nimmt. ro Zu 
3 Rn T : ’ 7 
u 3 x \ u 
ET 


eichungen (23), daß für die 


DV Er, 


2, 


iz 


Anw(z, y,2) = 


un = E 
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Auch. die in diesem Abschnitt betrachteten _Ver- 
fahren lassen sich leicht auf Differentialgleichungen 
höherer Ordnung und Systeme solcher Gleichungen 
verallgemeinern. 


4. Verschiedene der geschilderten Methoden wurden 
an der Diffentialgleichung 


zu? — Yy’tea+td)y — y=0 


BR > 
el NA 1) = |/ — 
=), vwn=y- 

(mit der exakten Lösung y = V x? + 5) erprobt. 


Insbesondere zeigte sich, daß das Interpolationsver- 
fahren nach 2. und nach 3. etwa den gleichen Zeitauf- 
wand erfordert. I. a. wird die Wahl des Verfahrens 
von der Form der Differentialgleichung abhängen. Je 


i 3 Fr E 
nachdem, ob die Funktion F,ryr, oder Fr 


leichter zu berechnen ist, wird man ein Verfahren nach 
1. bzw. 2.oder nach 3.nehmen. Bei der Wahl zwischen 
1. und 2. wird man i. a. wohl 2. bevorzugen, da der 
Rechenaufwand für beide Verfahren stets ungefähr 
gleich ist und man sich bei 2. die Argumente x,., bei l. 
nur die Parameterwerte i7; vorgeben kann. 
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Vereinfachte Herleitung der Weissingerschen 
Zirkulationsgleichung für den schiebenden 
Tragflügel. 


I. Einleitung. 


Die Schwierigkeiten, die sich bei der Aufstellung 
einer der Prandtlschen Zirkulationsgleichung ent- 
sprechenden Bestimmungsgleichung für schräg nen: 
strömte Tragflügel ergeben, sind bekanntlich in der 
Singularität begründet, die der Ausdruck für die indu- 
zierte Geschwindigkeit am Ort der tragenden Linie be- 
sitzt. Während Weissinger!)für denschiebenden 
Flügel diese Singularität berücksichtigt, indem er eine 
passende Anfahrströmung heranzieht, soll sie hier zur 
Verbesserung der Übersicht explizit ermittelt und 
durch eine geeignete ebene Schicht dem Flügelschnitt 
zugeordneter Wirbel kompensiert werden. 


II. Bereehnung der induzierten Ge- 
schwindigkeit desschrägen Wirbel- 
bandes. 


Die Aufwärtsgeschwindigkeit w in einem Punkte 
(2, Y, z), herrührend von einem in x-Richtung ge- 
legenen Wirbelband von der Breite b, welches sich von 
der tragenden Linie x(y) ab bis ins positiv Unendliche 
erstreckt, ist gegeben durch 


+b/2 
j N 
IMülerer 
x (1 52 ee) d 
Men) — 2)? +n — y® +2 


3) J.Weissinger, Derschiebende Tragflügel bei gesunder 
Strömung. Jahrbuch der deutschen Luftfahrtforschung 1940. 


- —bl2 


v 


Kleine Mitteilungen 65 


wenn I/”(n) die Dichte der abschwinimenden Wirbel 

bedeutet, mit der das Band längs der Breitenkoordi- 
: b b 

nate n belegt ist. — 7 ne cY 


Liegt der Aufpunkt oberhalb der tragenden Linie, 
so gilt speziell für den schiebenden geraden Flügel 
wegen & =ntgö, 2 = ytgö 

+b/2 
, MEN 
inw— I’(n) — — = ' 
“ (n—-y’+? 
—b/2 
Vo-wi+t)+2 
Spalten wir im Folgenden von allen mit z — O0 endlich 
bleibenden Ausdrücken ihren Grenzwert ab, so haben 
wir zunächst mit 


Sue 
f I’(n) dn 
inw= I X 
1 ar (2), 
—»]2 
— 4n(w ai w,) — sin6- 
+5/2 
2 n — y)?d 
T'(n) BEN. Lo, 
AR ((n— y)°+22) Y(n — y)? +2? cos? 6 


wobei 0(z) eine Größe kennzeichnen soll, die mit z 
gegen Null geht. Weiterhin gilt dann auch 


+b/2 


— 4r(w— w,) =sind | TmM-T'y) dn 
m — | 


—b/2 


+b/2 
Ir (mans j (n — y)’ dn 


((n-y)?+2?) VY(n-w)? + 2?cos?ö 


—b/2 
+0(z). 


Setzen wir zur Abkürzung 
+bj2 
— 4rw, =sinö | At) dn (3), 
—b/2 


so finden wir für den Rest nach elementarer Inte- 
gration 


—4n(w —w— w) 
= I*(y)sinö |In(n —y+Y(n —y)®+ 2? cos? 6) 
+b/2 


1, Vn—w?+zcos2ö+(n— y)sind 


en ein = 
2sinö Yn — y)? + 2? c0s?d — (n — y)sinö 


n 
- [2 


Der erste Teil in der Klammer ergibt 


Ma — 
Be — ") au 2 v)+ x? cos? ö 


In Sr 


b b F 2 2 
-\a t3)+y)\5 +y) +2*eostö 


„er GEN + 


2? cos? ö 


+ in (&+,)+ VE +) +#oso] 


2? cos? ö 
mit z— O0 im Zähler 


er) 


2? cos? ö 


In 
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Der zweite Teilliefert den Grenzwert 


1 1+sinö 
sind 1-—sind' 


Insgesamt ergibt sich demnach 


—4r(w— w-— w)=I'(y)sinö- 


er. 


sin Ö 


1-+sinö 
1 — sind 


)+oi. 


2? cos? ö 
Setzen wir noch 


— 4 w, = I’(y)sinö- 


Ir (a & ) ; 


A l1-+ sin ) (4), 


(y)cos®®ö sind 1-—sinö, 
und 
: 4 )z] 
— Arw, = —2T’(y)sinöln——, I(y)vgl.(7) (5), 
3 (Y) iy) (y)vgl.(7) (5) 
so wird mit den Abkürzungen (1) bis (5) 
vw ru Tw rw 0) 2... 66). 


III. Bestimmung der zugehörigen 
Zirkulation im Flügelschnitt. 

Die vom Wirbelband an der tragenden Linie indu- 
zierte Geschwindigkeit ist in den Bestandteilen w,, 
w, und w,endlich, während w, mit z— Ologarithmisch 
über alle Grenzen wächst. Diese Singularität läßt sich 
auf folgende Weise berücksichtigen. Denken wir uns 
den schiefen tragenden Wirbelfaden, wie in Wirklich- 
keit, über die tragende Fläche ausgebreitet und deren 
Wirkung an der Stelle yfürz — 0 wiein der Prandtl- 
schen Theorie durch die Wirkung des ebenen Flügel- 
schnitts ersetzt, so läßt sich eine Verteilung ebener 
gebundener Wirbel über die Tiefe /(y) angeben, die 
diese Singularität gerade kompensiert und eine be- 
stimmte Zirkulation I, zur Folge hat. 

Da die endlich bleibenden Geschwindigkeiten die 
Bestimmung der Zirkulation im Prandtlschen Sinne 
erlauben, gilt dann die Zirkulationsgleichung 

w+ w + w 
Fin al) Plan + pe)tne m 
wenn /(y) die Flügeltiefe, V die Anströmungsgeschwin- 
digkeit und x(y) den Anstellwinkel an der Stelle y be- 
deutet, senkrecht zur tragenden Linie gemessen. 

Zur Bestimmung der kompensierenden ebenen Ver- 
teilung über die Sehne des Flügelschnitts von der 
Länge !(y) gehen wir aus von der Formel 


” [2 
2nwe -- (Oo rat 
ö 


die die induzierte Aufwärtsgeschwindigkeit in einem 
Punkte oberhalb der Hinterkante liefert. Durch Teil- 
integration findet man, wenn y endlich bleibt, 


I _ E12 _L za\d 
2nw= (7 In Me 3 
4 0 


(8), 


Setzen wir 
+) =0Cp 
wobei die Variable p eingeführt sein soll durch 


E= (1 - cosg), 


normale) 
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so ergibt sich 
nr 


a 0a, Im yaresor + dp 
0 


und mit z > O im zweiten Bestandteil 


4 lz| 


IIWe = Car = 2 Me (10). 


Die Geschwindigkeit w, kann also durch w, aufge- 
hoben werden, wenn gilt 


Kae ? 
[63 re (y) sin ö. 
Die zugehörige Zirkulation 7, ergibt sich dann zu 
Ky) 
a F d a ir I(y) = 6) 11 
BA IE I IE Ed BEA U) Bar na (ıl). 
0 


Damit erhalten wir schließlich die folgende Bestim- 
mungsgleichung für die Zirkulation 


+5/2 
1 | I’(n)dn 


Urs) 


ne 64 
'(y) sin ö| 
RT iu 


rd witz]. 


12(y) cos? ö 


sindoö 1-—sinö 


die die Weissingersche Zirkulationsgleichung ohne Be- 
rücksichtigung der Querströmung darstellt. 
Hannover. K. Jaeckel. 


Der wirksame Austrittsquerschnitt 
bei Zyklonen. N 
Um für die Auslegung von Zyklonen einen Anhalt 
über die Größe des beim Austritt sich bildenden Kerns 
zu gewinnen, soll die gesamte Strömung als isentro- 
pische Expansion vom Eintritt bis zum Austritt auf- 
gefaßt werden, so daß nach Cr oc co eine wirbelfreie 
Strömung vorliegt; aus dieser Annahme folgen für die 
azimutale Komponente wy und für die axiale Kom- % 
ponente w, der Strömungsgeschwindigkeit w im Aus- @ 
trittsrohr die Beziehungen 
WZErWVcX = Ur... 
= Ws rkgen...oe (2) 
(r: Abstand von der Rohrachse, = <X(w, df) 
df: Flächenelement eines achsensenkrechten Schnitte® 
des Austrittsrohres); hierbei ist vorausgesetzt, daß im 
Austrittsquerschnitt in allen Meridianschnitten gleiche 
Strömungsverhältnisse vorliegen und weiterhin die 
Änderung der radialen Geschwindigkeitskomponente 
w, in Achsenrichtung vernachlässigbar sei. Ferner 
werde angenommen, daß infolge Druckausgleiches im 
Kern der Druck p.„ des Außenraumes herrsche. Die 
Konstanten C und K sind durch die Bedingung des 
Unveränderlichkeit des Massendurchsatzes und des 
Dralls bestimmt. Da durch ein Flächenelement 
df=r:dr:dp pro Zeiteinheit die Masse (n: Flächen- 


owondf=gKrdrdp 


strömt (o: Dichte, @: Azimut), so liefert der Konti- 
nuitätssatz die Gleichung 


. (3) 


(I+eE)m=2rKfordr 
TK 
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mo: die durch die Eintrittsöffnung des Zyklons in der 
Zeiteinheit strömende Masse, € m,: die bei der Ver- 
brennung im Zyklon in der Zeiteinheit erzeugte Gas- 
masse, "x: Kernhalbmesser, «: Austrittsrohrhalb- 
messer; da ferner der resultierende Vektor 
D=/S T7xdmü® 

des Drehimpulses, welcher durch die Austrittsöffnung 
in der Zeiteinheit transportiert wird, in Achsenrich- 
tung weist, brauchen bei der Integration nur die 
axialen Komponenten (7; x w;), des Vektorproduktes 
berücksichtigt zu werden. Somit wird 


. 4 
D=2nK/o(rwcosy)rdr 
"K 


. (4) 


(R,: Halbmesser bei Eintritt, Zyklon, w,: Eintritts- 
geschwindigkeit), so daß C = wu wird. Unter Be- 


a 
=2nKCford=m BoW 
"K 


rücksichtigung des Energiesatzes 


”»—1 


2 e cent 
w° = Wmax ( — (2) > 
\Poo 


(wmax: Endgeschwindigkeit beim Ausströmen ins 
Vakuum vom Zustand 2%, 090, der im wesentlichen 
mit dem Eintrittszustand p,, o, identisch ist: 


1 % . . A 9 
SL ET WE ET a el Fe 


T, 
i: Enthalpie, A: mechanisches Wärmeäquivalent, AT: 
bei Feuerungszyklonen die Temperaturänderung bei 
konstantem Druck) und der Änderung der Dichte 
bzw. des Druckes bei isentropischer Expansion kann 
wegen 


. (5) 


we = Kr Ct 2.2... (6) 
die Dichte als Funktion von r ausgedrückt werden. 
Wird X und C nach (3) und (4) in (6) fürr =rx, 
w= wr unter Beachtung von (5) eingesetzt, so er- 
geben elementare Umformungen für den Eintritt- 
querschnitt #5 wegen 

m = FEW 
folgenden Ausdruck 


FF 1 [6 w ax? a —1 
Ba EEG, Dies Seth teen) PT) 
„—1 
mit $ 
“ als 
WE, &,) = ( 4) du (8) 
& 
! 
A4= ERS . (9) 
Yı-=(@) 4 
y Poo 
= LEERE IE 
a=(lrE) Ro: Wo - & ( ) 


Die Größe «& stellt die Unbekannte dar, welche durch 
numerische Rechnung aus (7) zu bestimmen ist. — 
Für x = 1,4 (Luft) ist I elementar auswertbar: 


m,9-Yı-Z2hı-2)+e+8) 
1 
1-ı-4 
5 ® 
ST N ern 
IH Via 
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Nach Ermittlung von « ergibt sich der Kernhalbmesser 
zu 


FRE: 1 Wy a 
Ro FE) wmas FEN), 
1—.o? 5 4 
Poo 
Die aus (7) und (11) folgende Beziehung 
”—1 
Po | * : 
1 << (?e) (12) 
Pa 
zeigt, daß für die numerische Rechnung nur ein 
kleiner Zahlenbereich zu untersuchen ist. — Für den 


Durchsatz durch den Zyklon ist der Umstand von Be- 
deutung, daß mit wachsendem «a die untere Grenze in 
(8) stark steigt und zusätzlich die obere Grenze 
schwach fällt, so daß der Wert des Integrals abnimmt; 
für, = a wird nach (9) und (10) bei 


”—1l 
GR 
Pa 
”—l 


1+ 5 l Ro wo )( =) 2 
+ € awmazx Pa 
der Durchsatz zu Null; a, liegt unterhalb der oberen 
Grenze von (12). 

Mit steigendem r.:@ kann der Kernquerschnitt in- 
folge der bisher vernachlässigten inneren Reibung des 
Mediums einen zunehmenden Anteil des aufgezwunge- 
nen Dralls m, R, w, übernehmen. Falls in der äußeren 
Zone des Kerns infolge Impulsaustausches eine axiale 
Geschwindigkeit auftritt und somit aus Gründen der 
Kontinuität längs der Rohrachse eine entgegengesetzte 
Strömung aus dem Außenraum in den Kern vorhan- 
den ist, kann unter der Annahme eines für turbulente 
Rohrströmung üblichen Potenzgesetzes die axiale 
Komponente im Kern iin der Form 


Y 1/n 
1i= we -c (1- ) ) 
TK: 


geschrieben werden; ferner kann unter der Voraus- 
setzung konstanter Winkelgeschwindigkeit im Kern 


(13) 


. (14) 


Y 


WW — 
ar PK TK 


YES IB) 


gesetzt werden. Soweit diese Ansätze zur Darstellung 
der Kernströmung anwendbar sind, ergibt sich für den 
durch den Kern in der Zeiteinheit transportierten 
Drall D; nach kurzer Zwischenrechnung unter Beach- 
tung der für w,, = w;(r,) und wo, = wy(r,) gelten- 
den Beziehungen (1) und (2) 


De EEE mas 
h er Fr Wo & 0 


Die Konstanten c,,, sind hierbei durch die Euler- 
schen Integrale erster Gattung 


= 1222141) 


1 rk max Vo? — 1.2) 


»(4 1+2) 
(2 er 2) 
n 


gegeben; im Falle des 1/7-Gesetzes vonv. Kärmän 
wird 

6, = 0,0392. 
Bei relativ kleinem Kernquerschnitt ist somit der nach 
(16) berechnete Wert des durch den Kern übertrage- 


nen Drehimpulses vernachlässigbar. 
Berlin-Lichtenberg. W. Bader. 
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BUCHBESPRECHUNGEN 


M. Eßlinger, Statische Berechnung 
von Kesselböden. VIII+ 100 S. m. 21 Abb. 
Berlin-Göttingen-Heidelberg 1952. Springer-Verlag. 
Preis kart. 10,50 DM. 

Verf. hat sich die dankbare Aufgabe gestellt, die 
Berechnungsverfahren für Kesselböden einer ein- 
gehenden Untersuchung zu unterziehen. Hierbei ist 
einerseits eine kritische Sichtung des bisherigen 
Schrifttums vorgenommen worden, andererseits 
werden verschiedene Ergänzungen gebracht, die sich 
auf die Art der Wölbung des Bodens und seine Be- 
festigung, sowie auf Durchbrüche beziehen. Der 
1. Abschnitt bringt einen Überblick über die Rech- 
nungsarten; der 2. Abschnitt enthält die Herleitung 
von Einzelformeln für Boden, Krempe und Mantel 
(Tabellen und Kurven); der 3. Abschnitt enthält 
Anwendungsformeln für verschiedene Kesselformen 
auf Grund der Überlegungen des 1. und 2. Ab- 
schnittes. Der 4. Abschnitt bringt Zahlenbeispiele, 
welche im 5. Abschnitt z.T. mit Messungen von 
Siebel verglichen werden. Im 6. Abschnitt wird 
der Spannungsverlauf für die durchgerechneten Bei- 
spiele diskutiert. Das Buch schließt eine Lücke inner- 
halb der Konstruktions- und Berechnungsbücher 
und wird deshalb in der Praxis gute Aufnahme finden 


Dresden. H.Neuber. 


Dr. W. Krull (o. Prof. a. d. Univ. Bonn), Ele- 
mentare u. klassische Algebra vom 
modernen Standpunkt, Bd.I, 2. Aufl. 
(Sammlung Göschen, Band 930.) 136 S. Berlin 1952. 
Verlag Walter de Gruyter u. Co. Preis brosch 2,40 DM. 

Es handelt sich um den bemerkenswerten Versuch, 
mit den Begriffen der modernen Algebra dadurch 
vertraut zu machen, daß sie in die Darstellung der 
elementaren und klassischen Algebra eingeführt 
werden. Dem formalen Rechnen wird von Anfang 
an der Begriff des Körpers zugrunde gelegt. Das dem 
Bändchen gesetzte Ziel ist namentlich in den ersten 
Abschnitten recht gut erreicht. _ 

Ob freilich auch der Anfänger ein durchweg 
klares Bild erhält, ist bei der gedrängten Darstellung 
fraglich. Der knappe Raum für den sehr reichen 
Inhalt läßt manchen Ansatz willkürlich erscheinen, 
Aus diesem Grunde kann auch der Gruppenbegriff 
in seiner fundamentalen Bedeutung nicht so heraus- 
treten, wie es wünschenswert wäre. e 

Die Anlehnung des die höhere Gleichungstheorie 
behandelnden Abschnittes an die geschichtliche Ent- 
wicklung ist ein immanenter Ausgleich des bezeich- 
neten unvermeidlichen Mangels. Dem Leser wird 
dadurch der Weg gewiesen, mittels der einschlägigen 
Literatur etwa verbliebene Lücken in der begriff- 
lichen Klarheit zu schließen. Dazu könnte übrigens 
auch der vorgesehene zweite, Aufgaben enthaltende 
Band die erforderlichen Anleitungen geben. Erst 
nach seinem Erscheinen ist ein abschließendes Ur- 
teil möglich. 

Dresden, Draeger. 


H. M. Hansen (Prof. a. d. Univ. von Michigan) 
und P. F. Chenea (assoc. Prof. a. d. Uriv. von Michi- 
gan, Mechanics of Vibration. 4178. m. 
345 Abb. New York und London 1952, J. Wiley and 
Sors, Crapmen a. Hall. Preis 8$. 

‘Das Buch entstand aus einem schwingungstech- 
nischen Kursus an der technischen Abtcilung der 
Universität Michigan. Die theoretischen Grund- 
lagen sind nur in knapper, auf die unmittelbare An- 
wendung hinzielender Form gebracht, um auch Le- 


sern mit geringer mathematischer Vorbildung ent- 


gegenzukommen. Dennoch gehen einzelne Teile des 
Buches zu abstrakteren Formulierungen über, welche 


auch für wissenschaftlich geschulte Leser in anre- 
gender Weise abgeleitet werden. Nach einführenden 
Bemerkungen über Definitionen, Schwingungsarten, 
Dimensionen und Maßeinheiten, sowie kinematischen 
Überlegungen bei einfachen harmonischen Bewe- 
gungen folgen 3 Hauptteile. Der 1. Teil behandelt 
Schwingungen mit einem - Freiheitsgrad: Freie 
Schwingungen ohne Dämpfung (Ableitung und Lö- 
sung der Bewegungsgleichungen für Translation und 
Rotation, zusammengesetzte Federn, Anwendungen, 
Energiemethode); erzwungene Schwingungen ohne 
Dämpfung (Natur der erzwungenen Schwingung, 
Ableitung und Lösung der Bewegungsgleichung, Re- 
sonanz, Meßinstrumente, erzwungene Pendelschwin- 
gung, kritische Drehzahlen); erzwungene Schwin- 
gungen mit Dämpfung (Natur der Dämpfung, Sy- 
steme mit linearer Dämpfung, Ableitung und Lö- 
sung der Bewegungsgleichung, Resonanz, Energie- 
bilanz). Der 2. Teil behandelt Systeme mit meh- 
reren Freiheitsgraden: Klassische Methode (Art 
der Freiheitsgrade, Ableitung und Lösung der Be- 
wegungsgleichungen für freie Schwingungen, er- 
zwungene Schwingungen ohne Dämpfung bei 2 Frei- 
heitsgraden, desgleichen mit Dämpfung, Orthogo- 
nalität der Hauptschwingungsarten, Schwingungen 
vieler Freiheitsgrade, Schwingungen von Systemen 
mit Zwangsbedingungen und die Lagrange’schen 
Gleichungen, potentielle Energie und Dissipations- 
funktion, Anwendung auf ein Zentrifugalpendel, An- 
wendung auf ein dissipatives System); Bewegungs- 
methode und Einführung komplexer Variablen (der 
Vektor und seine komplexe Darstellung, Impedanz, 
Anwendungen, verzweigte Systeme, Übertragbarkeit, 
Systeme mit ebener Bewegung, Schwingungen einer 
federnd abgestützten Plattform); Lösung der all- 
gemeinen Frequenzgleichung (Allgemeines, Verfahren 
von Holzer, Verfahren von Graeffe). Im 
3, Teil werden spezielle Probleme behandelt; Sy- 
steme mit kontinuierlicher Massenverteilung (Ein- 
führung, Ableitung und Lösung der Bewegungs- 
gleichungen, Schwingungen einer Saite, Längs- 
schwingungen eines Stabes, _Drehschwingungen, 
Biegeschwingungen, Anwendungen der Energie- 
methode); nicht-harmonische Schwingungen (Natur 
der nicht-harmonischen Bewegung, freie Schwin- 
gungen mit Dämpfung, Wirkung eines zeitlich un- 


stetigen Kraftverlaufes) ; nicht-Jineare Schwingungen 


(Natur der nicht-linearen Systeme, ungedämpfte 


freie Schwingungen bei nicht-linearer Federkraft, 


Bewegung eines Pendels bei großen Amplituden, 
stückweise lineare Federkraft, nicht-Jineare Dämp- 
fung). Es folgen Ausführungen über die Durch- 
rechnung von Aufgaben, die mit den im Textteil 
aufgeführten Problemen in engem Zusammenhang 
stehen. Das Buch berührt weitgehend die für die 
Praxis maßgeblichen Problemstellungen und bringt 
eine große Zahl mit klaren Skizzen versehener Bei- 
spiele, die sich an die verschiedensten Gebiete der 
technischen Praxis anlehnen und zum besseren Ver- 
ständnis, sowie zum Selbststudium des Lesers bei- 
tragen; daher kann ein guter Erfolg erwartet werden, 


Dresden. H. Neuber. 


L. Locher-Ernst, Einführungindiefreie 
Geometrie ebener Kurven. (Elemente 
der Mathematik vom }öFeren Standpunkt aus. Bd. TI) 
88S. mit 167 Abb. Verlag BirkLäuser. Basel 1952. 
Preis brosch. 12.50 Sch w. Fr. 

Die genannte Bandreihe des Bareler Verlags Birk- 
häuser will einige der schönsten Gebiete der Mathe- 
matik in leicht faßlichen und im Umfang mäßigen 

erblicken sowohl Fachmathematikern als auch 


. Freunden der Mathematik darbieten. Der vorliegende _ 
erste Band dieser Reihe bringt eine bisher fehlende _ 
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elementare Einführung in einen modernen und in- 
teressanten Zweig der Geometrie. Einer seiner be- 
deutendsten Vertreter, J. Hjelmslev, nannte 
ihn „natürliche Geometrie“. Die Entwicklung der 
ersten Grundlagen dieser rein anschaulichen Geo- 
metrie geht in der Hauptsache auf A. Kneser und 
C. Juel zurück. 

Locher-Ernsts Einführung in dieses Gebiet 
setzt keine besonderen Kenntnisse voraus. Ihr Ziel 
"besteht darin, die den Verlauf ebener Kurven regeln- 
den Gesetze kennen zu lernen und eine Übersicht der 
einfachsten Kurvenformen zu gewinnen. Die pro- 
jektiv-topologischen Eigenschaften der ebenen Kurven 
werden aus dem unmittelbar Evidenten in begriff- 
licher Strenge entwickelt. Als frei bezeichnet der 
Verfasser diesen Zweig der Geometrie, weil seine 
Untersuchungen frei von einem Koordinatensystem 
geführt werden. Die Grundlage seiner Entwicklungen 
ist der Begriff der ‚elementaren stetigen Bewegung“. 
An die Spitze sind einige einfachste, anschaulich 
selbstverständliche Aussagen axiomatischen Charak- 
ters gestellt, darunter ein Stetigkeitsaxiom. Aus 
diesen lassen sich weitere einfache Sätze projektiver 
Natur beweisen. Es wird hierauf definiert, was man 
unter den Begriffen ‚elementarer Bogen‘, ‚elemen- 
tare Kurve‘ und ‚einfacher Bogen‘ zu verstehen hat. 
Der letztere wird, was die bei ihm möglichen Singu- 
laritäten und seine allgemeine Struktur anbetrifft, 
genauer untersucht, worauf die weiteren Begriffe 
„C-Bogen“, „Eilinie“ und ‚‚Spiralbogen‘“ entwickelt 
werden können, 

Bei der Ableitung allgemeiner Sätze über elemen- 
tare Kurven erweisen sich manche Eigenschaften, die 
man für wesentlich algebraischer Natur anzusehen 
geneigt wäre, als viel allgemeinerer Natur. Der Auf- 
bau der freien Geometrie wird bis zu der Konstruktion 
der Gebilde dritter Ordnung bzw. Klasse mit Hilfe 
von Eilinien durchgeführt. Diese Konstruktion ist 
wie manches andere in dem Werk ein Beitrag des Ver- 
fassers zu dem neuen Zweig der Geometrie. Der Vor- 
trag ist bei aller Exaktheit leicht verständlich, wozu 
die zahlreichen Zeichnungen wesentlich beitragen. 
Die beigefügten Aufgaben sind geeignet, den Leser zu 
selbständiger Beschäftigung anzuregen. | 


Freiberg (Sa.) = W. Schmid. 


F. Riegels. Die Strömung um schlanke, 
fastdrehsymmetrische Körper. (Mit- 
teilungen aus dem Max-Planck-Institut für Strö- 
mungsforschung N. 5.) 72S. mit 12 Abb. u. 4 Tab, 
Göttingen 1952. Selbstverlag Max-Planck-Institut 
für Strömungsforschung. Preis brosch. 7,— DM. 

Mit der Arbeit von Fuhrmann ‚Theoretische 
experimentelle Untersuchungen an Ballonmodellen‘“ 
begann im Jahre 1910 eine Reihe von Untersuchungen, 

' welche die Strömung um rotationssymmetrische Kör- 
per mit der Methode der Achsbelegung durch Quellen 
und Senken bei rein axialer Anströmung, bzw. durch 
Dipole bei Queranströmung behandeln. 

Die Berechnung der Geschwindigkeitsverteilungum 
Rotatienskörper beliebiger Gestalt gelingt in der 
vorliegenden Arbeit verhältnismäßig bequem durch 
Belegung der Körperoberfläche mit Quell- 
ringen, deren Intensitätsverlauf längs der Körper- 
oberfläche aus einer Integralgleichung iterativ be- 
stimmt wird, und zwar entsprechen Quellringe mit 
konstanter Intensität über den Umfang der axialen 
Anströmung, solche mit einer bestimmten periodischen 
der Queranströmung. Begünstigt wird diese Behand- 
lung durch bereits vorliegende Tabellen des Ver- 
fassers, welche die Berechnung der Geschwindigkeits- 
komponenten der genannten Quellringe stark verein- 
fachen. Die Methode wird zunächst auf das Rotations- 
und das dreiachsige Ellipsoid angewandt, für die be- 
reits die exakte Lösung vorliegt, und ergibt für alle 

_  Anströmungsrichtungen bereits in der ersten Nähe- 

rung befriedigende Übereinstimmung mit der exakten 


Lösung. Weiterhin wird noch ein Körper behandelt, 
dessen Querschnitte dreieckig mäßig deformiert, in 
Längsrichtung elliptisch veränderlich sind. Hierbei 
ergibt sich, daß die Berechnung für beliebig vorge- 
gebene Formen, soweit sie im Rahmen der theoreti- 
schen Voraussetzungen bleiben, keine zusätzliche Er- 
schwerung bedeutet. Auch gerade hierdurch ist die 
Arbeit als ein erheblicher Fortschritt in dem anfäng- 
lich genannten Aufgabengebiet zu bezeichnen. 


Weilburg. In S:chiller., 


K.Hain, Angewandte Getriebelehre.‘ 
408 S. mit 491 Abb. Hannover-Darmstadt 1952, 
Hermann Schroedel Verlag KG. Preis brosch. 21,50 DM, 
geb. 24,— DM. 

Das Buch gibt die zur Untersuchung vorhandener 
oder zum Entwurf neuer Getriebe notwendigen, be- 
kannten und durch eigene Untersuchungen des Ver- 
fassers ergänzten Verfahren und Methoden an. Es 
beschränkt sich dabei bewußt auf periodische ebene 
Getriebe, so daß z. B. die Rädertriebe nicht behandelt 
werden, dafür aber neben den Kurventrieben be- 
sonders ausführlich die Vierzylinderkette (Gelenk- 
viereck) und die n-Zylinderketten (n= 6, Sund mehr). 
Die Darstellung bevorzugt die geometrischen Methoden 
und sieht oft von den Herleitungen benutzter For- 
meln oder Beziehungen ab, was für den Lernenden 
eine gewisse Erschwerung bedeutet. Aber es wird vor 
allem Wert darauf gelegt, daß der Praktiker die Ver- 
fahren findet, die notwendig sind und die er auf Grund 
der angegebenen Beschreibungen auch anwenden 
kann. Manche Verfahren werden nur kurz, doch unter 
Angabe entsprechender Literatur, gestreift. Das aus- 
führliche (doch nicht vollständige) Literaturverzeich- 
nis von 269 Nummern zeugt von dem Fleiß und von 
der Mühe, die der Verfasser in sein Buch gelegt hat, 
inein Buch, das für den Getriebepraktiker von großem 
Wert sein wird. Ausstattung und Druck des Buches 
sind gut. 


Aachen. W.Meyer zur Capellen. 

Dr. phil. H. v. Sanden (o. Prof. a. d. TH Hannover), 
Praktische Mathematik. (Teubners Ma- 
thematische Leitfäden, Bd. 44). Zweite, erweiterte 
Aufl. 120 S. mit 25 Abb. Leipzig 1951. B. G. Teub- 
ner Verlagsgesellschaft. Preis kart. 3,80 DM. 

Die erste Auflage dieses Buches wurde bereits in 
dieser Zeitschrift, 29 (1949), S. 191, besprochen. In 
der vorliegenden Neuauflage wurden zwei Para- 
graphen, über das Iterationsverfahren und über die 
Ausgleichung mit Gewichten, hinzugefügt und der 
Abschnitt über harmonische Analyse und Synthese 
(unter direktem Bezug auf die Tafeln von Zipperer) 
erweitert. Außerdem wurden einige Beispiele neu 
aufgenommen. Bei der Durchsicht ist man immer 
wieder überrascht von der Fülle der praktischen 
Rechenerfahrung, die Verf. in das Buch hineingear- 
beitet hat und die in anderen Büchern meist nicht in 
dieser Art vermittelt wird. Dem Leser wird nicht nur 
Aufklärung über sein Problem gegeben, sondern er 
wird darüber hinaus oft noch auf eine Zahl weiterer 
Gesichtspunkte hingewiesen, die ihm sonst vielleicht 
entgangen wären. Nicht nur der Student, sondern 
gerade auch der schon längere Zeit in der Praxis 
stehende Ingenieur, der während seiner Ausbildung 
nichts über praktische Mathematik gehört hat, wird 
das Buch mit großem Gewinn für seine Arbeit be- 
nutzen. 


Dresden. G.Opitz. 


Dr.-Ing. W. Meyer zur Capellen, Instrumen- 
telle Mathematikfür den Ingenieur. 
(Fachbücher für Ingenieure.) 3838. m. 190 Abb. 
Essen 1952. Verlag W. Giradet. Preis geb. 27,80 DM. 

Das vorliegende Werk will den Ingenieur mit 
den Anwendungsmöglichkeiten der handelsüblichen 
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mathematischen Instrumente vertraut machen. Des- 
wegen ist das Buch nicht nach den verschiedenen 
Arten der Instrumente gegliedert, sondern es sind 
nacheinander ihre Verwendungsmöglichkeiten in den 
folgenden vier Gebieten, Arithmetik und Algebra 
(92 S.), Geometrie und Trigonometrie (13 S.), Infi- 
nitesimalrechnung (215 S.) und Harmonische Analyse 
und verwandte Probleme (54 S.) behandelt. Natür- 
lich wird auch der Aufbau und die Theorie der In- 
strumente wie die Fehlermöglichkeiten bei ihrer An- 
wendung erläutert. Das Hauptgewicht liegt aber auf 
der Beschreibung ihres Gebrauchs bei den verschieden- 
sten Problemen, die allen möglichen Gebieten ins- 
besondere der Ingenieurwissenschaften entnommen 
sind. Dabeiist in manchen Fällen gezeigt, wie dieselbe 
Aufgabe mit verschiedenen Instrumenten behandelt 
werden kann. Die Lösungen sind im allgemeinen bis 
in alle Einzelheiten durchgeführt. Ein Schrifttums- 
verzeichnis, das in der Hauptsache die Anwendungen 
berücksichtigt und ein alphabetisches Sachverzeichnis 
schließen das Buch. Offenbar hat der Verfasser eine 
sehr große Arbeit an das Buch gewendet. Schon allein 
die Sammlung der zahlreichen schönen Beispiele und 
ihre Bearbeitung wird viel Zeit gekostet haben. So ist 
aber auch ein Buch entstanden, aus dem sich der 
Ingenieur in vielen Fällen Rat holen kann. 


Dresden. Willers. 


Technische Hausmitteilungen des Nordwestdeutschen 
Rundfunks. Sonderheft:Unterlagenfür UKW- 
Netzplanungen. 44S. mit 39 Abb. Hamburg 
1952. H.H. Nölke Verlag. Preis brosch. 6,50 DM. 

Das Heft enthält drei Aufsätze. In den beiden 
ersten wird über amerikanische Arbeiten und über 
deren Anwendung bei den westdeutschen Planungen 
auf dem UKW-Rundfunk- und Fernsehgebiet be- 
richtet. Die dritte führt den mathematisch wenig 
geschulten Leser in die wahrscheinlichkeitstheore- 
tischen Grundlagen des in den ersten Arbeiten be- 
handelten Stoffes ein. 


Dresden. Willers. 


6. Tintner (Prof. am Iowa State College of Agri- 
culture and Mechanie Arts), Econometrics. 
XIII + 370 S. New York 1952. John Wiley & Sons. 
Preis 5,75$. 

Das Buch führt den Leser in einen Zweig der Wirt- 
schaftswissenschaft ein, der in Deutschland erst in 
jüngster Zeit Beachtung gefunden hat. Die Ökono- 
metrie nimmt eine Mittelstellung zwischen statisti- 
scher und mathematischer Ökonomie ein und sucht 
durch Anwendung moderner statistischer Methoden 
auf'theoretische Modelle quantitative Resultate zu 
erzielen und die Richtigkeit ökonomischer Theorien 
nachzuweisen. Verf. betont den rein feststellenden 
Charakter dieser Vorgehensweise; er ist sich der 
Grenzen (Widerspruch zwischen Modell und Wirk- 
lichkeit) und Mängel (Fehlen der für die statistischen 
Schlüsse nötigen Voraussetzungen) der Methode voll 
bewußt und schildert die dadurch auftretenden 
Schwierigkeiten in lebendiger Weise. Es spricht mehr 
für den Autor als für den dargestellten Gegenstand, 
wenn man-am Schluß fast jedes der zahlreichen Bei- 
spiele Wendungen wie ‚‚das Ergebnis muß mit ge- 
bührender Vorsicht gedeutet werden‘ oder ‚diese Er- 
gebnisse sind wiederum von zweifelhafter Gültigkeit“ 
findet; die kritische Einstellung des Verf. tritt z. B. 
auch darin zutage, daß er bezweifelt, daß eine rein 
zahlenmäßige statistische Erhebung möglich ist ohne 
eine zugrunde liegende ökonomische Konzeption oder 
Theorie (— die freilich nicht formuliert zu sein 
braucht —), oder daß die Ökonometrie wesentlich zur 
Erforschung von Entwicklungsvorgängen beitragen 
kann. Trotzdem steht für den Verf. die praktische Be- 
deutung der Ökonometrie — von der „laissez-faire 
economy“ bis zur zentral geplanten Wirtschaft — 
außer Frage, 


Den Leser dieser Zeitschrift wird vor allem inter- 
essieren, welche mathematischen Hilfsmittel der 
Ökonometrie als Werkzeug dienen: 

Es handelt sich einmal im wesentlichen um die 
Herausarbeitung linearer Regressionen zwischen 
mehreren Veränderlichen und den Nachweis der Signi- 
fikanz der Regressionskoeffizienten. Die abgeleiteten 
Beziehungen werden als Angebots-, Nachfrage-, Pro- 
duktionsfunktionen usw. gedeutet, wobei die Deu- 
tung eine weitere Problematik mit sich bringt. Breiten 
Raum nehmen in diesem Zusammenhang die ‚‚ge- 
wogene“ Regression und die „Identifikation“ von 
Gleichungen ein, zwei Probleme, die speziell für die 
Ökonometrie von Bedeutung sind, da dort die Mög- 
lichkeit zu Experimenten fehlt. Außerdem werden 
noch die Diskriminanz-Analyse und verwandte Fragen 
behandelt. 

Ein zweites großes Gebiet ist die Untersuchung von 
Zeitreihen, also Trendberechnung (wobei die logisti- 
sche Kurve immer noch eine Rolle zu spielen scheint), 
Abspaltung von oszillatorischen und periodischen Be- 
standteilen, Autokorrelation, Korrelogrammanalyse 
usw. Besondere Bedeutung mißt Verf. den stochasti- 
schen Differenzengleichungen bei, mit deren Hilfe er 
sämtliche bei Zeitreihen auftretenden Phänomene ein- 
heitlich zu erfassen hofft. 

In einem Anhang werden die Elemente der Ma- 
trizen- und Determinantentheorie und entsprechende 
Rechenmethoden (lineare Gleichungssysteme, Eigen- 
werte) gebracht. An Voraussetzungen wird vom Leser 
die Kenntnis ökonomischer Theorien und der Grund- 
tatsachen der Statistik verlangt. 

Das sonst sehr breit geschriebene Buch verzichtet 
auf Beweise und gibt dafür jeweils die Originalliteratur 
an. Das Namenverzeichnis geht über neun Spalten und 
legt Zeugnis ab von dem großen Fleiß, mit dem das 
umfangreiche Material zusammengetragen wurde. So 
ist ein Buch entstanden, das einmal einen Überblick 
über das bisher Erreichte gibt, zum anderen aber viele 
Fragestellungen aufzeigt, die bisher noch unerledigt 
sind — vor allem bei den Prüfmethoden für kurze 
Zeitreihen —, und damit für die künftige Forschung 
richtungweisend ist, wobei die zu erwartenden Ergeb- 
nisse natürlich nicht allein der Ökonometrie zugute 
kommen werden. 


Dresden. G. Opitz. 


Dr.L.Vietoris (o. Prof. a. d. Univ. Innsbruck) u. 
Dr. 6.Lochs, Vorlesungen über Diffe- 
rential- undIntegralrecehnung. 4158. 
mit 159 Abb. Innsbruck 1951. Universitätsverlag 
Wagner. Preis geb. S138,—, DM 30,—. 

Bei der Abfassung des Werkes, das aus Vorlesungen 
an der Technischen Hochschule Wien und der Uni- 
versität Innsbruck entstanden ist, haben sich die 
Verff. von folgenden Grundsätzen leiten lassen: 
‚„‚1. Der Leser soll mit nichts Unnötigem belastet 
werden. 2. Auf Klarheit der Begriffe ist das größte 
Gewicht zulegen. 3. Verständlichkeit ist wichtiger als 
Kürze. 4. Beweise sind die beste Einübung der Be- 
griffe. 5. Unrichtige Gedankengänge können nicht 
verstanden werden; daher darf nirgends versucht 
werden, dem Leser durch Aufnahme von, Unrichtig- 
keiten das Verständnis scheinbar zu erleichtern.‘“ — 
Man darf sagen, daß die konsequente Einhaltung 


dieser Grundsätze ein Lehrbuch geschaffen hat, das - 


in seiner geschickten Stoffauswahl und seinem 
methodischen Aufbau dem Studenten wertvolle Hilfe 
beim Durcharbeiten seiner Vorlesungen geben wird 
und auch zum Selbststudium vortrefflich geeignet ist. 

Der Inhalt entspricht etwa dem Stoff, der an einer 
TH.in den Vorlesungen: Höhere Mathematik I und II 
vorgetragen zu werden pflegt, vertieft durch eine 
gründliche Einführung der analytischen Grundbegriffe 
im Reellen und im Komplexen; auf eine ausführliche 
Erörterung des Zahlbegriffes, insbes. der Irrational- 
zahlen wird verzichtet. Im Einzelnen sei als bemerkens- 
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wert hervorgehoben, daß der Differentialquotient und 
das Integral nahezu parallel entwickelt werden, was 
in mancher Hinsicht den Aufbau der Analysis er- 
leichtert. Es sei z.B. auf den natürlichen Loga- 
rithmus, der als bestimmtes Integral eingeführt wird, 
und (als dessen . Umkehrung) ‚die Exponential- 
funktion, ferner auf die Herleitung der Ta ylorschen 
Formel hingewiesen. 

Inhaltsübersicht: Grundbegriffe und vor- 
bereitende Darlegungen. — Die Ableitung oder der 
Differentialquotient. — Wichtige Begriffe und Hilfs- 
sätze. — Das Integral. — Geometrische Anwendungen. 
Winkelfunktionen. — Der Taylorsche Satz. — Un- 
bestimmte Formen. — Numerische Methoden. — 
Einiges aus der Algebra. — Integration der rationalen 
und verwandter Funktionen. — Analysis im Kom- 
plexen. — Das Wallis$che Produkt und die Stir- 
lingsche Formel. — Einiges aus der analytischen 
Geometrie des Raumes (u. a. Determinanten und Vek- 
toren). — Raumkurven. — Funktionen mehrerer Ver- 
änderlicher. — Implizite Funktionen. — Mehrere 
Funktionen mehrerer Veränderlicher. 

Dresden. Maruhn, 

Review of Elecetronie Digital Computers, 
Joint AIEE — IRE Computer Conference. 114 8. 


Eingegangene Bücher 


[at 


m. 100 Abb. New York 1952. American Institute 
of Electrical Engineers. Preis 3,50 $. 

Das Heft berichtet über eine vom American Insti- 
tute of Electrical Engineers und dem Institute of 
Radio Engineers im Dezember 51 gemeinsam in Phila- 
delphia abgehaltene Tagung über Rechenautomaten, 
die den Zweck hatte, Richtlinien für ihre Weiterent- 
wicklung zu gewinnen. In 19 Vortragsberichten wird 
eingehend die Erfahrung dargelegt, die bei Entwurf, 
Einrichtung und Gebrauch von 10 großen elektroni- 
schen Maschinen verschiedener Bauart gewonnen 
wurde, teils handelt es sich dabei um Parallel- teils 
um Serienmaschinen. Als Speicher wurden in der 
Hauptsache Ultraschallwellen-, Magnettrommel- oder 
Kathodenstrahlröhrenspeicher verwendet, wobei letz- 
tere sowohl das Strich-Punkt-System, wie das System 
der Fokussierung-Defokussierung benutzten. Dahei 
kamen sowohl die Erbauer wie die Benutzer zu Wort. 
In dem letzten Aufsatz, der die gegenwärtige und 
künftige Entwicklung der Rechenautomaten behan- 
delt, wird die Ansicht ausgesprochen, daß man erst 
am Anfang der Entwicklung stehe. Das Heft gibt eine 
ausgezeichnete, mit Bildern und schematischen Zeich- 
nungen reich ausgestattete Übersicht über den Stand 
der Entwicklung Ende des Jahres 1951. 


Dresden. Willers. 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 


EINGEGANGENE BÜCHER 


Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten). 


L. Bieberbach (ehem. o. Prof. d. Univers. Berlin), 
Theorie der geometrischen, Kon- 
struktionen (Lehrbücher und Monographien 
aus dem Gebiete der exakten Wissenschaften. 
Mathem. Reihe Band 13.) 102 S.m. 102 Abb. Basel 
1952. Verlag Birkhäuser. Preis brosch. 15,60 Fr. 
geb. 18,70 Fr. 


T.vanLangendonck, Funcoes ortogonais 
na resolucäo de problemas da tteoria 
da elasticidade, Tomo|]l, Generali- 
dades eTorcao.69S. mit 28 Abb. Sao Paulo 
1952. 

Dr. Bruno Eck, Ventilatoren. Entwur 
und Betrieb der Schleuder- und 
Schraubengebläse. Zweite verbesserte und 
erweiterte Auflage. XI + 304 S. mit 344 Abb. Ber- 
lin/Göttingen/Heidelberg 1952. Springer-Verlag. 
Preis geb. 36,— DM. 


Dr. Hans Reichardt, Gesetzmäßigkeiten 
der freien Turbulenz. VDI-Forschungsheft 
414. Zweite Auflage. 30.8. mit 18 Abb. Düsseldorf 
1951. Deutscher Ingenieur-Verlag G.m.b.H. Preis 
12,50 DM. 


Dr.R.M.E.Badrawy,EbenePlattengitter 
beilüberschallgeschwindigkeit. (Mit- 
teilungen aus dem Institut für Aerodynamik a. d. 
Eidgen. Techn. Hochschule in Zürich. Herausgegeben 
von Prof. Dr. J. Ackeret. No.19) 908. mit 
71 Abb. Zürich 1952. Verlag Leemann. 


0. Haupt (o. Prof. a. d. Univ. Erlangen), Ein- 
führungin die Algebra. TeilI. (Mathe- 
matik und ihre Anwendungen in Physik und Technik, 


Reihe A, Band 5, Teil 1). Zweite Auflage. XVIII + 
3708. Leipzig 1952. Akademische Verlagsgesell- 
schaft Geest & Portig K. G. Preis geb. 23,— DM. 


- Henry de Wolf Smyth, Atom-Energieund 

ihreVerwertungimKriege, übersetztund 
erläutert von Prof. Dr. Friedrich Dessauer. 
3528. Basel 1947. Ernst Reinhardt Verlag. Preis 
kart. 11,15 DM., geb. 14,50 DM. 


Hermes-Scholz, Mathematische Logik. 
(Enzyklopädie der mathematischen Wissenschaften 
Band 1,, Heft 1, Teil I.) 2. völlig neubearbeitete Aufl, 
8258. Leipzig 1952. B.G. Teubner Verlagsgesell- 
schaft. Preis geh. 8,20 DM. 


Dr. K.H. Weise (o. Prof. a. d. Univ. Kiel), Vom 
mathematischen Denken. (Veröffent- 
lichungen der Schleswig-Holsteinschen Universitäts- 
gesellsch. Neue Folge Nr. 2.) 30 S. Kiel 1952. Ver- 
lag Ferd. Hirt. Preis geh. 1,20 DM. 


M.A.N. Forschungsheft 1952 1. Halbjahr, heraus- 
geg. von der Maschinenfabrik Augsburg-Nürnberg. 
112 S. mit vielen Abbildungen. Augsburg-Nürnberg 
1952. Preis geh. 10,— DM. 


Dr. Fr. v. Krbek (Prof. a. d. Univ. Greifswald), 
Eingefangenes Unendlich. Bekenntnis 
zur Geschichte der Mathematik. VI + 3318. mit 
128 Abb. Leipzig 1952. Akademische Verlagsgesell- 
schaft Geest & Portig K. G. Preis geb. 22,— DM. 


Dr. W. Lietzmann (Prof. a. d. Univ. Göttingen), 
Anschauliche Einführung in die 
mehrdimensionale Geometrie. 2208. 
mit 157 Abb. und 1 Tafel. München 1952. Verlag 
R. Oldenbourg. Preis brosch. 19,50 DM. 
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NACHRICHTEN 


Paul Felix Nemenyi 
June 5, 1895 (Budapest) — March 1, 1952-(Washing- 
ton, D. C.). 

Dr.Ing. habil. Nemenyi was Head of the 
Theoretical Mechanies Section at the U.S. Naval 
Research Laboratory. During a life of travel and 
change in much of Europe and America he had devoted 
his major effort to pure mechanies. Among his dis- 
coveriesin elasticity, publishedinthisZeitschrift, 
are a general duality theorem (9, 488; 10, 383) and a 
beautiful analogy with plane irrotational flow (12, 
364): the latter states that given any plane isometric 
net, there exists a five parameter familiy of stress 
systems for which it is the net of stress trajectories. 
He reduced tbe entire theory of extensional stresses in 
thin shells to the single equation uy/un + (n? — 1) 
f’If = 0 (Proc. Nat. Acad. Sei. U.S.A. 29, 159 (1943)), 
enabling the exact solution for many shell shapes 
subject to unsymmetrical load and support. A strong 
exponent of inverse and semi-inverse methods, in 
collaboration with R. Prim he discovered many 
exact rotational and three-dimensional solutions of 
the equations of gas dynamics. He was an expert in 
hydrology, fishways, and the morphology of rivers. 
His numerous expository works are of unusual value 
for their organization andscope (e. g.this Zeitschrift11, 
(1921) S.59 u.5.450 Wasserbauliche Strö- 
mungslehre (Barth, Leipzig 1933, Adv. Appl. 
Mech. 2, 123 (1951)). Toward the end of his life, many 
of his ideas were worked out by his pupils. 

Bloomington, Indiana. C. Truesdell, 


Tagung der Deutschen Mathematiker-Vereinigung 
in München vom 4. bis 8. September 1952, 


Die Tagung war mit rund 150 Teilnehmern, unter 
denen sich auch eine Reihe namhafter ausländischer 
Gäste befand, gut besucht. Es wurden rund 50 Vor- 
träge gehalten, darunter diefolgenden Hauptvorträge 
von dreißig bis sechsig Minuten Dauer: 

Hofmann- Ichenhausen: Mathematisches aus 
früheren Tagen der altbayerischen Hochschulen, 

Kneser- Tübingen: Amerikanische Arbeiten 
zur mathematischen Wirtschaftslehre und Soziologie. 

Heinhold- München: Die Behandlung mathe- 
matischer Probleme mit Hilfe moderner Rechengeräte. 

Witt - Hamburg: Der Begriff des Filters als ver- 
allgemeinerte Folge. 

Rinow- Greifswald: Bericht über die innere 
Flächentheorie von A. D. Alexandrov. 

Godement-Nancy: Aus der Arbeit des Bour- 
baki-Kreises. 

Die übrigen Vorträge behandelten Fragen aus fol- 
genden Gebieten: Algebra, Gruppentheorie, Zahlen- 


theorie, Mengenlehre, Limitierungstheorie, Kom- 
plexe Funktionentheorie, Gewöhnliche und partielle 
Differentialgleichungen mit Rand- und Eigenwert- 
problemen, Geometrie einschließlich Differential- 
geometrie und Topologie sowie andere Einzelfragen. 

Ein Nachmittag war der Erörterung von Fragen 
des Mathematischen Unterrichts gewidmet. | Dabei 
sprachen: 

Behnke- Münster: Zur Neugründung der 
IMUK. Die Verantwortlichkeit der DMV für die 
Gestaltung des mathematischen Schulunterrichts. 

Cramer-Kempfenhausen: Die Aufgabe der 
Mathematik im Gymnasialunterricht und die Idee 
des mathematisch-naturwissenschaftlichen r Gym- 
nasiums. 

Mönkemeyer- Braunschweig: Die Lage der 
Höheren Schule in Niedersachsen mit Ausblicken auf 
einige Nachbarländer. 

Heinrich- München: Die Lage des Mathe- 
matikunterrichts in Bayern. 

Lösch. 


Kolloquium über programmgesteuerte Rechengeräte- 


und Integrieranlagen in der Rheinisch - West- 
phälischen Technischen Hochschule in Aachen am 
22. Juli 1952. 


Diese Tagung, die erste ihrer Art in Deutschland, 
war von etwa 200 Personen besucht, darunter zahl- 
reichen Gästen aus Industrie und Technik. 

Dr. habil.H. Bückner, Minden berichtete über 
den von ihm konstruierten z.Zt. in Entwicklung be- 
findlichen „Integromat“. Prof. Dr. A. Walther, 
Darmstadt sprach über Rechentechnik in Darmstadt 
— vorhandene Geräte, Ergebnisse und Pläne. Dr. 
F.J.Weyl, Experte der US-Navy, hielt einen 
Vortrag über die Entwicklung der Rechenautomaten 
in den Vereinigten Staaten. Prof. Dr. L.Bier- 
mann, Göttingen schilderte die Entwicklung elek- 
tronischer Rechenautomaten in Göttingen, wo, ge- 
messen an den im Vergleich zu den USA zur Ver- 
fügung stehenden gerinsen Mitteln, Erstaunliches 
geleistet wurde. Dipl.-Ing. K. Zuse, Neunkirchen 
gab neben einer Schilderung der historischen Ent- 
wicklung eine Übersicht über die z. Z. bei der Firma 
Zuse K. G. laufenden Entwicklungsarbeiten. 


Aachen. H.Cremer. 


Die näckste Tagung der Gesellschaft für angew. 
Mathematik und Mechanikfindet in der dritten Woche 
nach Ostern, vom 21.—25. April 1953, in Aachen 
statt. (Am 21. April: Sitzung des Vorstandsrates 
und Begrüßungsabend.) — Die örtliche Tagungs- 
leitung hat Herr Prof. Dr. Schultz-Grunow 
übernommen. 


ZUSCHRIFTEN AN DEN HERAUSGEBER 


Zu Hsien Chih Liu: Über die Entstehung 
von Ringwellen an einer Flüssigkeitsober- 
fläche durch unter dieser gelegene, kugelige, 
periodische Quellsysteme. Z. angew. Math. Mech. 
32 (1952), S. 211 —226. 

Das vom Verfasser behandelte Problem ist schon 
1922 von Lamb gelöst worden (Proc. London Mathe- 
matical Society, Second Series, Vol. 21, 1923); 
kompliziertere Quell-Systeme hat Havelock unter- 
sucht (Phil. Mag. 1942). Nach Abfassung der vor- 


liegenden Arbeit (1945) sind dann zahlreiche weiter- 
gehende Untersuchungen gemacht worden. Gegen- 
über der Lambschen Darstellungsweise ist die vom 
Verfasser gegebene Ableitung und Art der Diskussion 
etwas kompliziert; die Ergebnisse der letzteren gehen 
jedoch hinsichtlich der Wellenformation wesentlich 
über das bisher Bekannte hinaus. Insofern hat die 


Arbeit ihren Wert. Es ist jedoch erstaunlich, daß 


dem Verfasser die Literatur unbekannt geblieben ist. 
Hamburg. 'G. Weinblum. 
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